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Recoil protons in a cloud chamber have been used to measure the energy spectrum of the 


neutrons from C"%+H?. The neutrons are found to be monochromatic when produced by 
2-Mev deuterons. An accurate Q value of —0.27+0.02 Mev is computed from the range of 


the recoil protons. This Q value gives a mass of 13.00988--0.00004 for N¥. 





INTRODUCTION 


Bare neutrons from carbon bombarded by 
deuterons have been used as a monochro- 
matic source of fast neutrons of variable energy. 
In the present experiment the energy spectrum 
of these neutrons has been measured to deter- 
mine whether or not the neutrons from C™ are 
in fact monochromatic. 


EXPERIMENTAL ARRANGEMENT 


A thin carbon target was prepared by de- 
positing camphor black on a silver disk. Its 
thickness was determined by weighing and by 
comparing the yield of neutrons from it with the 
yield from a thick graphite target, using the 
known excitation curve! for the reaction to effect 
the comparison. Its thickness was found to be 
75 kev for deuterons of two million volts. 

The target was bombarded with 2.00-Mev 
deuterons from The Rice Institute pressure 
electrostatic generator. The energy of the neu- 


* This work was done in June 1942 under OSRD con- 
tract. 

** Now at Illinois Institute of Technology. 

*** Now at University of Wisconsin. 
.' T. W. Bonner, E. Hudspeth, and W. E. Bennett, Phys. 
Rev. 58, 185 (1940). 
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trons in the forward direction was measured by 
stereoscopically photographing the recoil protons 
in a cloud chamber filled to half an atmosphere 
with methane and 95 percent ethyl alcohol. The 
cloud chamber was 25 cm in diameter, and its 
near side was 12 cm from the source of neutrons. 
Tracks in a strip 8 cm wide from front to rear 
of the chamber could be measured. 







RESULTS 


The lengths of tracks of recoil protons in the 
cloud chamber were measured when the direction 
of the track was within 15° of the line drawn 
from the source of neutrons. Tracks as short as 
one centimeter in length could be measured 
reliably. The lengths were converted to lengths 
in standard air using (a) the pressure of gas in 
the cloud chamber when expanded as measured 
on a Bourdon gauge, (b) the vapor pressure of 
95 percent alcohol and water as given in the 
International Critical Tables, and (c) the atomic 
stopping powers of the gases in the cloud cham- 
ber as given by Livingston and Bethe.? Account 
was taken of the fact that the methane is 


* M. S. Livingston and H. A. Bethe, Rev. Mod. Phys. 9, 
274 (1937). ' 
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Fic. 1. Geometry corrections: The dotted curve is the 
correction for the finite length of the cloud chamber while 
the solid curve is for the finite depth of the illuminated 
portion of the chamber. S is the observed range in the 
cloud chamber. 


actually natural gas and is 13 percent ethane. 
This was checked by determining the density of 
the gas from the methane cylinder. 

Track lengths were converted to neutron 
energies using the Cornell range-energy curve for 
protons. The numbers of tracks in each energy 
interval of 0.05 Mev were counted. These num- 
bers were corrected for the variation with energy 
of the cross section for neutron-proton scattering 
and also for geometrical factors which made more 
probabie the measurement of short tracks than 
of long ones. 

Longer recoil protons have smaller probability 
of being measured because they leave the cham- 
ber at the distant wall or at the top and bottom 
of the illuminated space. (If they left the side of 
the stereoscopic space they could still be meas- 
ured on one of the stereoscopic photographs.) 
Hence, two geometrical corrections were re- 
quired. The first, for the finite length of the 
chamber, is 


W1i=(1/Lo)—1/(Lo+A—S), (1) 


where W, is the probability of observing a track 
of length S in the rectangular illuminated area 
of a cloud chamber of diameter \. Lo is the 
distance from the source to the chamber and S 
is the range in the chamber of the recoil proton. 
The correction for the finite depth of the chamber 
is of the simple form 


Wa=1-—4S sin0o/3rd (2) 


if the source is infinitely distant from the cham- 


(orbitrory scale) 
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Fic. 2. The energy spectrum of neutrons from a thin 
target of carbon bombarded by 2.0-Mev deuterons, 4 
methane filled cloud chamber at 0° to the deuteron beam 
was used for energy measurement. 


ber. % is the maximum angle that the recoil 
proton can make with the neutron and still be 
acceptable for measurement; d is the depth of 
the illuminated portion of the chamber. How. 
ever, if the neutron source is not infinitely 
distant, use of the above expression discriminates 
against the longer tracks. The correct expression 
for Wa when the neutron source is at a finite 
distance from the chamber may be obtained in 
a manner similar to that used in deriving the 
above expression ;* however, the resulting expres- 
sion is then rather involved and was actually 
solved numerically for only five different track 
lengths. A smooth curve was drawn through 
these five values and the resulting correction 
curve (normalized to one for a 10 cm track) is 
plotted on Fig. 1. 

The final energy spectrum of neutrons from 
C”+D is shown in Fig. 2. For comparison, the 
dotted curve shows the width of a group of 
monochromatic neutrons, the straggling being 
caused by (a) the thickness of the target, (b) the 
fluctuations in voltage of the Van de Graaff 
machine, (c) the angular straggling of the neu- 
trons which could reach the cloud chamber, (d) 
the angular straggling of the protons, which were 
measured out to 15°, and (e) the range straggling. 
The asymmetry on the low energy side of the 
group is probably to be attributed to small angle 
elastic scattering of neutrons in the wall and 
roof of the cloud chamber on the side toward 


*H. T. Richards, Phys. Rev. 59, 796 (1941). The ex- 
pression derived in this reference is for a different geometry 
than used in the present cloud-chamber measurements. 
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the neutron source. The wall was one-quarter 
inch of iron and the roof was one inch of glass. 

The constant background at low energies 
probably represents scattering by the walls of 
the room rather than slow neutrons from the 
source. A blank run was taken with no carbon 
on the silver disk and the number of neutrons 
observed was negligible. Therefore, not more 
than a few percent of the neutrons observed 
could have been produced by the molecular beam 
in the analyzer. 

The graph shows that the mean energy of the 
neutron group is 1.60 Mev. This becomes 1.65 
Mev when corrected because recoil protons were 
measured out to 15°. The mean deuteron energy 
in the target was 1.96 Mev. This gives an energy 
evolution of @= —0.27+0.02 Mev in the reac- 
tion. 
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CONCLUSIONS 


The neutrons from carbon bombarded by deu- 
terons are monochromatic, except for the high 
energy neutrons from C", when produced at 
2 Mev; and are therefore likely to be mono- 
chromatic when produced at lower bombarding 
voltages. 

The Q of the reaction is —0.27+0.02 Mev. 
This value is believed to be more accurate than 
previous determinations.‘ This Q value gives a 
mass of 13.00988+-0.00004 for N™. 


4 J. D. Cockcroft and W. B. Lewis, Proc. Roy. Soc. 154, 
261 (1936). Q=—0.28 Mev (from theshold of reaction). 

T. W. Bonner, Phys. Rev. 53, 496 (1938). Q=—0.25 
+0.05 Mev (from recoil protons in cloud chamber). 

W. E. Bennett, T. W. Bonner, E. Hudspeth, H. T. 
Richards, and B. E. Watt, Phys. Rev. 59, 781 (1941). 
Q=-—0.19+0.05 Mev (from recoil protons in cloud 
chamber). 
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The sensitivity of the theoretical neutron-proton scattering cross section to possible varia- 


tions in the quantities defining the potential has been investigated in the energy range from 
0 to 6 Mev. It is found that in this energy range the variations in exchange or tensor character 
and range of the potential which are consistent with what is known about the interaction do 
not result in so large a modification of the predicted cross section as do variations in shape. 

A well resembling a Yukawa potential (A exp(—ar)/r) is found to yield a cross section at 
6 Mev which is 10 percent lower than that given by a square well fitted to the epithermal 
neutron cross section and the binding energy of the deuteron. The Yukawa well is in better 
agreement with experiment than is the square well. One of the constants entering into the 
determination of the potential is the epithermal neutron scattering cross section in hydrogen. 
We have determined this by extrapolating the data of Frisch with the aid of the theory, and 


obtained a value of 20.8 10-* cm? 


INTRODUCTION 


URRENT nuclear theories attempt to cor- 
relate all known facts about nuclear scatter- 

ing and nuclear binding by means of the assump- 
tion that these phenomena are caused by a 
potential energy of interaction between the 
fundamental particles. Unfortunately, the present 


data are too meager and too crude to permit the 


* Now at Princeton University. 
** Contract No. W-7405-Eng-48. 






determination of this potential in any but the 
most schematic way. As a result, it is desirable 
to find out how critically the predictions for 
experimental results depend on each of the many 
parameters needed to specify the potential, and 
thus to indicate the range and accuracy of 
experiments needed to determine any feature of 
the potential within the limits desired. 

In this paper, the neutron-proton scattering 
cross section is investigated in the energy range 
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from 0 to 6 Mev. The object of this investigation 
is to estimate the extent to which the predicted 
cross sections will vary when the potential is 
given the largest possible modifications which 
seem to be consistent with all data other than 
that of neutron-proton scattering. When this 
work was done, there was so little systematic 
information about the neutron-proton scattering 
cross section, that our emphasis was placed on 
estimating the uncertainties in the predicted 
cross section, on the basis of- data then in 
existence. Since then, fairly complete measure- 
ments in the energy range have been made 
available,? and our results can now be used not 
only to provide closer limits on the potential, 
but also to indicate the direction in which the 
theory may be modified, if necessary, to obtain 
better agreement with experiments. 

The four main independently specificable fea- 

tures of the potential are: (a) range, (b) fraction 
of tensor force, (c) exchange character, and (d) 
shape. The depths of singlet and triplet wells in 
neutron-proton interaction are most conveni- 
.ently determined in terms of the other four 
features with the aid of two reliably khown 
constants, the binding energy of the deuteron, 
and the epithermal neutron-proton scattering 
cross section. 

Probably the feature of the potential about 
which the least is known is the shape. In fact, 
it is found that a square well fits all data other 
than that of neutron-proton scattering above 3 
Mev about as well as any other shape. Since the 
modifications in cross section resulting from 
likely modifications in shape are not expected to 
be large, the analytic simplicity of the square 
well may be utilized advantageously in investi- 
gating the sensitivity of the cross section to 
features other than the shape. The study of 
effects of variations in shape will then be deferred 
until the -other three parameters have been 
taken up. 


I. VARIATION IN RANGE 


The range of nuclear forces is usually deter- 
mined by means of a detailed investigation of the 


1C. L. Bailey, William E. Bennet, T. Bergstrahl, Richard 
G. Nucholls, H. T. Richards, and J. H. Williams, Phys. 
Rev. 70, 118 (1946); Phys. Rev. 70, 583 (1946). : 

? Bretscher (to be published). 
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variation of cross section with energy. In this 
way, Breit, Thaxton, and Eisenbud* have been 
able to estimate the range of proton-proton forces 
with considerable accuracy, obtaining a value of 
(2.8+0.4) X10-" cm?, with a square well. In the 
absence of grounds to the contrary, the tentative 
assumption is usually made that the range of 
neutron-proton forces is the same as that of 
proton-proton forces. This assumption, at pres. 
ent, obtains some qualitative support from the 
general features of nuclear binding, but its main 
quantitative support is twofold. First, it jg 
found** that if the neutron-proton singlet well jg 
taken to be the same as that for proton-proton 
interaction, the epithermal neutron-proton scat- 
tering cross section and other low energy neutron: 
proton cross sections, can be fitted within experi- 
mental error. Second, calculations of Rarita and 
Schwinger® show that a square well of radius 
2.8X10-" cm? is consistent with the observed 
quadripole moment of the deuteron. That the 
range of neutron-proton forces is not much less 
than that of proton-proton forces is made clear 
by the calculations of Sen,*® which indicate that 
a well of radius 2.4X10-" cm is the narrowest 
that can yield simultaneously the correct quadru- 
pole moment and magnetic moment. 

In the interests of seeking the simplest hy- 
pothesis consistent with what is now known 
about nuclear interactions, it seems desirable to 
retain the assumption of equality of all nuclear 
forces, and to try to obtain agreement with 
experiment within the limits of this hypothesis. 
If this assumption is made, it is then of interest 
to investigate the variations in the neutron- 
proton cross section which may result from the 
uncertainty of the proton-proton forces, as ob- 
tained from proton-proton cross sections. In ac- 
cordance with proton-proton data, a range of 
2.8 X10-" cm was taken to be the most probable, 
with permissable variations as large as 0.410 
cm. Calculations showed that this uncertainty 
in the radius introduced an uncertainty of 5 
percent in the cross section at 6 Mev, with 


( ‘ me Thaxton, and Eisenbud, Phys. Rev. 55, 1018 
1 x 
‘L. E. Hoisington, S. S. Share, and G. Breit, Phys. Rev. 
56, 884 (1939). 

5 W. Rarita and J. Schwinger, Phys. Rev. 59, 436 (1941) 
and 59, 556 (1941). 

* Sen (unpublished). 
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Fic. 1. Upper curves are drawn with expanded ordinate 
in order to show the extent of uncertainty arising from 
lack of knowledge of the exchange character of nuclear 


forces. 


correspondingly smaller uncertainties at lower 
energy. The cross section is larger, the larger the 
nuclear radius. The reason for the comparatively 
small sensitivity to nuclear radius, exemplified 
by the fact that a 5 percent change in cross 
section accompanies a 30 percent reduction in 
the “area” of the well, is that the depth is 
always adjusted to give the correct deuteron 
binding energy and epithermal neutron-proton 
cross section. 


Il. TENSOR FORCES 


Tensor forces are introduced into nuclear 
interactions in order to account for the observed 
quadrupole moment of the deuteron. Rarita and 
Schwinger® have shown that a variation in the 
fraction of tensor force of about 75 is consistent 
with the probable error in the quadripole mo- 
ment. It is, therefore, necessary to determine 
whether this uncertainty in the tensor force 
affects the cross section. It may be expected 
that the effect of this uncertainty will be small, 
because at moderate energies, the tensor force 
behaves very much like an addition to the 
ordinary radially symmetrical force, except for 
the production of a small D wave, which is re- 
sponsible for the quadrupole moment. The ex- 
pectation is borne out by the calculations of 
Rarita and Schwinger® and Sen,* who show that 
for a radius of 2.8X10-" cm, the tensor force 
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results in a reduction of only 2 percent in the 
triplet cross section at 0 Mev and at 2.82 Mev. 
Sen also shows that these cross sections are not 
changed appreciably when the fraction of tensor 
force is varied from 0.5 to 1. 

The effects of tensor forces at higher energies 
can be estimated from the results of Rarita and 
Schwinger,® who give the 3.S:+3D, cross section 
at 15.3 Mev, when the fraction of tensor force is 
adjusted to obtain agreement with the quadrupole 
moment of the deuteron. We find that this cross 
section is only 1 percent greater than that pre- 
dicted by ordinary forces at this energy. Hence, 
it seems safe to conclude that, in the range from 
0 to 6 Mev, the uncertainty in the fraction of 
tensor force should produce a negligible uncer- 
tainty in the cross section. In fact, even the 
complete neglect of the effects of tensor forces 
should result in an error of less than 2 percent in 
the cross section in this energy range. We have, 
therefore, decided that, for the purpose of study- 
ing the sensitivity of cross section to other 
parameters, the use of tensor forces is not neces- 
sary. Their effect will then be taken into account 
approximately, however, by means of a uniform 
reduction of 2 percent in all 3.S cross sections. 


Ill. EXCHANGE FORCES 


Exchange forces are introduced into nuclear 
interactions for the purpose of explaining the 
saturation character of nuclear binding. A fairly 
wide range of combinations of exchange and 
ordinary forces is, however, consistent with the 
present experimental data. It is, therefore, neces- 
sary to find out whether the cross section depends 
appreciably on the precise exchange character 
of the forces in the range from 0 to 6 Mev. 

The combinations studied in this paper were, 
in the terminology of Rarita and Schwinger, the 
“charged” and the ‘‘symmetric” theories, which 
represent about the widest possible variation in 
the exchange character consistent with the 
existence of saturation. The results of the calcu- 
lations showed that the contribution of the ‘‘ P” 
wave to the total cross section was so small that 
the exchange character did not make an appreci- 
able difference in this cross section. Even the 
angular dependence varied by less than 1 percent 
(see Fig. 1). In order to understand the smallness 
of the angular dependence, it is helpful to analyze 
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the sources of the P scattering. On account of 
the spin dependent tensor forces which operate 
only in the triplet state,® the phases are different 
for 3Po, 3P:, and 3P:2. The calculations showed 
that although the 3P» phase was fairly large from 
4 to 6 Mev, the angular dependence arising from 
it was almost cancelled by that arising from the 
3P;, and 3P: phases, which are smaller, but 
which occur with greater statistical weight. The 
reason for this cancellation is that the tensor 
forces average out to zero in the first approxi- 
mation. The net result is that practically all the 
angular dependence arises from the singlet state, 
which occurs with only one-fourth of the total 
weight. Consequently, both the angular depend- 
ence of the cross section and its uncertainty 
arising from the uncertainty in the exchange 
character, are very small below 6 Mev. 


IV. SHAPE 


Up to this point, the sensitivity of the cross 
section to possible variations in the parameters 
specifying the potential has been studied with 
the aid of the simplifying assumption of a square 
well. It is now necessary to estimate the effects 
of possible derivations from the square shape. 
The principal deviations from the square form 
which are significant for moderate energy scat- 
tering are (a) a tail and (b) a concentration of 
force near the origin. Hoisington, Share, and 
Breit* have studied such variations in the case 
of proton-proton scattering, and have found that 
the best agreement with experiment is obtained 
by making these two modifications simultane- 
ously. For example, they find that the expo- 


nential well, which is characterized mainly by a, 


long range tail without a great concentration 
of force near the origin, predicts both a wrong 
energy dependence for “‘.S’’ scattering and exces- 
sive “P’’ anomaly. On the other hand, the 
Yukawa potential (A exp(—ar)/r), which em- 
bodies both the tail and the necessary concen- 
tration of potential, is found to fit the proton- 
proton data as well as the square potential, while 
it is also in very good agreement with the epi- 
thermal neutron-proton cross section. It is neces- 
sary, however, to choose a range corresponding 
to a meson mass of 300 electron masses. The 
significance of the disagreement of this value 
with that of 200 electron masses obtained from 
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cosmic-ray experiments is not clear at present, 
particularly because of the uncertain status of 
the details of all current meson theories. In any 
case, it seems that a Yukawa potential with this 
range represents about the widest variation jp 
shape from that of a square .well which is cop. 
sistent with the hypothesis of equality of alj 
nuclear forces, and it is, therefore, the mogt 
promising one for study. 

In order to simplify the calculations, and to 
separate clearly the effects of the tail from thoge 
of the concentration of force near the origin, it jg 
obviously permissible to represent these two main 
features by an L-shaped potential, consisting of 
a deep short range well followed by a shallow 
long range well. The proportions of the L should 
be chosen to give a general resemblance to the 
above Yukawa potential in the range of impor- 
tance, which is in.the region from 10-" cm to 
5X10-" cm. This method of fitting yields some- 
what ambiguous results, but it turns out that 
almost any proportions of the L which lie within 
the range of ambiguity yield practically the same 
cross sections, if the depths are adjusted to yield 
the’ deuteron binding energy and epithermal 
neutron-proton cross sections. The small adjust- 
ment of singlet depth needed to give the epi- 
thermal neutron-proton cross section is easily 
within the range of the ambiguity in fitting. 
Since the L well now resembles the Yukawa 
potential which also gives the correct proton- 
proton and epithermal neutron-proton scattering, 
and because this agreement is not critically. 
dependent on the precise proportions of the well, 
it may be concluded that the small ambiguity 
in shape has no important effect on the cross 
section. 

A further requirement which the L-shaped 
potential must satisfy is that it gives the correct 
quadrupole moment of the deuteron. According 
to the work of Sen, the radius of the concentra- 
tion of the force near the origin which we have 
chosen (2.2X10-" cm) is definitely below that 
needed without a tail to give simultaneously the 
correct quadrupole moment and magnetic mo- 
ment. The tail is, however, very effective in 
producing quadrupole moment, which is propor- 
tional to the mean value of r*. Because the 
magnetic moment is proportional to r, the tail is 
less effective in producing it Approximate calcu- © 
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Fic. 2. Comparison of cross section for neutrons-proton scattering for two different well shapes. 
Specification of wells.—1. Square well-range 2.8 X 10-" cm. Singlet depth: 11.85 Mev; Triplet 
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depth: 21.2 Mev. 2. Narrow square well with tail: a. Square well-range 2.2 X 10-" cm. Singlet 
depth: 18.39 Mev. Triplet depth: 30.33 Mev. b. Square tail-range 2.2-5.6X 10-" cm. Single 


depth: 0.3 Mev. Triplet depth: 0.5 Mev. $ Experimental points. $ Williams, et al. A Bretscher, 


et al. A— Aoki, Z—Zinn, Seeley and Cohen. 


lations show that the well chosen by us’ can 
yield the correct quadrupole moment and mag- 
netic moment of the deuteron within experi- 
mental error, and that the fraction of tensor force 
needed to obtain the agreement is about the 
same as that obtained by Rarita and Schwinger® 
for a square well of radius 2.810-" cm. 


V. RESULTS 


The cross sections were calculated in the 
usual way, which is too well known to require 
further discussion.’ Our results (Fig. 2) show 
that the L-shaped well yields a cross section 
which is at 6 Mev, 10 percent lower than that 
predicted by the square well, with correspond- 
ingly smaller differences at lower energies. Below 
1 Mev, the difference in cross sections is practi- 
cally undetectable. It is of interest to understand 
why the cross sections behave in this manner. 
The agreement at low energies is, of course, a 
consequence of the fitting of the binding energy 
of the deuteron and the epithermal neutron- 
proton scattering cross section. At higher ener- 


gies, the ‘‘tail,”” which is only 0.5 Mev deep, 


™Mott and Maasey, Theory of Atomic Collisions, 
Chapter II. 









begins to have very little effect on the wave 
function, and the main part of the scattering 
is then due to the deep narrow well at the 
center. 

Because the tail is very effective at low ener- 
gies, however, the concentration of potential near 
the origin is not actually deep enough to produce 
by itself the correct binding or epithermal neu- 
tron scattering cross section. If the depth were 
such as to give the correct deuteron binding and 
epithermal scattering, then the narrowing of the 
effective radius would result in a considerable 
decrease in the triplet cross section, and a slight 
decrease in the singlet cross section.* The de- 
crease of effective depth, however, decreases the 
slope of the wave function at the edge of the 
well,* and thus introduces a tendency to increase 
the cross section. In the singlet state, this increase 
is the main effect, but in the triplet state, it is 
considerably over-balanced by the decrease re- 
sulting from the narrowing of the well. Because 
the increase of singlet cross section is less than 
the decrease of triplet cross section, and because 


of the greater statistical weight of the triplet 


*H. A. Bethe and R. F. Bacher, Rev. Mod. Phys. 8, 83 
(1936). 
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state, the net effect is to decrease the cross 
section at high energies. 


VI. COMPARISON WITH EXPERIMENTS 


The best early experiments are those of Aoki,® 
and those of Zinn, Seely, and Cohen,!® which 
cover a range from 2.1 to 2.9 Mev. Later experi- 
ments’? done after this work was completed, 
cover a range from 0.35 Mev to 6 Mev. In 
general, the experimental cross sections are in 
much better agreement with the L well than 
with the square well, wherever the two predic- 
tions differ. Although experimental errors are 
great enough to permit the assignment of a 
square well of range 2.4 10~—" cm, which would 
give a 6 Mev, a cross section 5 percent below 
that given by the well of range 2.8X10-" cm, 
the mean of the data is in definitely better 
agreement with the predictions of the L-shaped 
well. Besides, the range of 2.4 10-" cm is prob- 
ably the extreme limit of error in the proton- 
proton range, and not a particularly likely value. 
Since the square shape has nothing to recommend 
it, except analytic simplicity, it seems preferable 
to assume tentatively that the well resembles a 
Yukawa potential. 

One of the most important parameters used in 
the choice of a potential is the epithermal 
neutron-proton cross section. A very accurate 
value of this quantity was obtained by extrapo- 


* H. Aoki, Proc. Phys. Math. Soc. Japan 21, 232 (1939). 
1 W. H. Zinn, S. Seely, and V. W. Cohen, Phys. Rev. 
56, 260 (1939). 
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lating to zero energy with the aid of the 

the neutron-proton cross sections of Frisch, 
which were made in the energy range from 35 key 
to 490 kev. Because the uncertainties in the 
potential do not appreciably affect the crogs 
section at these low energies, the extrapolation 
is quite reliable. This is verified in Fig. 3, which 
shows that the variation of cross section with 
energy in this range is in very good agreement 


with the theory. The extrapolated value of the » 


epithermal cross section is 20.8X10-* cm?, as 
compared with the previously accepted value of 
(20+2) X10-* cm?.” 

More work on this problem is desirable, par. 
ticularly with regard to obtaining a more syste. 
matic study of the sensitivity of cross section to 
shape. Now that a Yukawa potential seems to 
give better agreement with experiment than a 
square well, an effort should be made to fit a 
Yukawa potential directly to all known data, 
including proton-proton scattering, neutron-pro- 
ton scattering, binding energy of the deuteron, 
and quadripole and magnetic moments of the 
deuteron. The prediction of the cross sections at 
energies above 6 Mev, where the exchange and 
tensor character of the forces is important, is 
particularly significant. Some theoretical work 
has already been done on this problem." It is 
necessary at higher energies to consider the 
possibility that the fundamental interactions 
cannot be described in terms of a potential, 
which is a function only of the relative coordi- 
nates. 

We are greatly indebted to Professor Bethe 
and Mr. Sen for sending us some very helpful 
manuscripts. It is also a pleasure to acknowledge 
gratefully the help and encouragement received 
from Professor Oppenheimer and Professor Serber 
in the solution of the problem. This work was 
carried out under the auspices of the Manhattan 
District. 
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2% V. W. Cohen, H. H. Goldsmith, and J. — my 
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13 J. Leites Lopes, Phys. Rev. 70, 5 (1946). A list of 
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Yield of Photo-Neutrons from U*** Fission Products in Heavy Water* 


S. BERNSTEIN, W. M. Preston, G. WoLFE, AND R. E. SLATTERY 
Clinton Laboratories, P. O. Box 1991, Knoxville, Tennessee 


(Received January 20, 1947) 


The photo-disintegration of the deuteron has been used to study the hard y-rays emitted 
by fission products of U**. The neutrons created in the process were used as the indicator of 
the presence of hard y-rays. The fission products were placed at the center of a 10” radius 
sphere of heavy water. Conclusions about the periods and yields of the hard y-rays were 
made from the total number of photo-neutrons being captured in a large tank of oil surrounding 
the sphere of heavy water. Eight half-lives were found: 2.5 sec., 41 sec., 2.4 min., 7.7 min., 
27 min., 1.6 hr., 4.4 hr., and 53 hr. The shortest one and longest one of these are least reliable. 
Eighty-five percent of the photo-neutrons appear in the two shortest half-lives, the 2.5-sec. 
component being three times as intense as the 41 sec. component. The total saturated activity 
of the photo-neutrons for an infinite amount of heavy water was calculated from the 10” 
radius sphere measurements to be about 16.5 percent of the saturated delayed neutron activity. 
It is calculated that there must be of the order of one to two photons of energy above 2.2 Mev 
emitted per fission by fission products with half-lives greater than one second. 





I. INTRODUCTION Il. APPARATUS 


Eye threshold property of the photo-disin- 
tegration of the deuteron can be used as a 
crude y-ray spectrometer. The reaction 


A schematic diagram of the apparatus is 
shown in Fig. 1. It consisted of a ten-inch radius, 
$-inch wall, aluminum sphere immersed in a 
large tank of oil. A “rabbit’’ containing a sample 


y+D—n+p 


requires that the energy of the y-ray be at least 
2.17 Mev! (the binding energy of the deuteron). 
The use of this reaction immediately ‘‘resolves”’ 
y-rays above the threshold from those below, 
independently of the relative intensities of the 
two groups. It is convenient to use the neutrons 
created in the process as an indicator of the 
presence of hard y-rays, since they may be 
slowed down by a moderating medium to 
thermal energies where they are easily detected. 
The total’ relative number of photo-neutrons 
produced is also a fairly easily determined quan- 
tity. It is a measure of the relative intensity of 
the hard y-rays. We have studied the periods 
and yields of photo-neutrons produced in heavy 
water by the y-rays from U** fission products 
and from these yields have estimated the cor- 
responding yields of the high energy y-rays 
themselves. 


* This document is based on work performed at Clinton 
Laboratory operated by Monsanto Chemical Company 
under contract No. W-35-058-eng-71 for the Manhattan 
Project and the information covered therein will appear in 
Division IV of the Manhattan Project Technical. Series, 
as part of the contribution of Clinton Laboratories. 

(1938) T. Rogers, and M. M. Rogers, Phys. Rev. 55, 263 


of enriched uranium oxide could be transferred 
from the center of the Clinton Laboratories 
chain-reacting pile? to the center of the sphere 
in about 0.25 sec. by means of a pneumatic tube. 
The sphere could be filled or emptied of heavy 
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Fic. 1. Schematic diagram of apparatus. 


*H. D. Smyth, Atomic Energy for en Purposes 


(Princeton University Press, Princeton, 1945), Chapter 


VI; or Rev. Mod. Phys. 17, 392 (1945). 
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water in about one minute by means of a siphon 
arrangement. The number of photo-neutrons 
created in the sphere of heavy water by the 
fission product y-rays was determined by study- 
ing the neutron distribution in the oil by means 
of fission chamber detectors. The chamber could 
be placed at chosen reproducible radial positions 
in the tank by means of a track and carriage 
arrangement. Since some of the periods involved 
are fractions of a second, it was necessary to keep 
the initial counting rates quite high in order to 
obtain good statistics throughout the run. An 
available amplifier and scaling circuit found to 
have a resolving time of 20 microseconds was 
adapted for use in the experiment. Exposure 
times and sample sizes were adjusted so that 
corrections for missed counts reached 5 percent 
only at the start of one or two of the runs. Fission 
pulses from the amplifier were fed into a scale of 
256. The data were recorded by photographing 
the scaling circuit interpolation bulbs with an 
open shutter oscilloscope camera. The film moved 
past the neon bulbs in such a manner that each 
bulb appeared on the film as a dashed line for the 
period the bulb was lit. A time signal was im- 
pressed on the film every 0.1 second by means of 
a neon bulb and a revolving shutter driven by 
a synchronous motor. Another signal lamp con- 
nected through the electropneumatic control 
switch recorded the length of an exposure. 
Figure 2 shows a reproduction of a typical film. 




























Il. METHOD 


The method employed to measure the source 
activity is based on the following theory. Con- 
sider a point source of Q neutrons/sec., having 
an arbitrary energy distribution (above thermal 
energy), surrounded by an oil bath large enough 
to absorb all the neutrons. It is known that the 
absorption in the oil is caused almost entirely 
by capture of thermal neutrons by the protons. 
The number of captures of neutrons of energies 
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Fic. 2. Sample of 


film recorded 
data. 


higher than thermal is negligible. In the steady 
state the rate of production of neutrons mys 
equal the rate at which they are captured, or 
Q= Sf (nv) .0dv, where (nv); is the flux of thermal 
neutrons at any point in the oil, o, is the macro. 
scopic cross section for capture at thermal 
energy, and the integral is extended over the 
oil volume. Measurements with a fission chamber 
detector give a counting rate Y at any point ip 
the oil which is proportional to (mv):, provided 
correction is made for the epicadmium counting 
rate. 

For spherical symmetry, Q« { YR°dR. In our 
case, since there is negligible neutron absorption 
in the heavy water sphere, the integral is taken 
from Ry to ©, where Ro is the sphere radius. 

The ion chamber was fixed in a given radial 
position, the sample was sent into the pile fora 
fixed exposure time, then sent to the center of the 
sphere, and the resulting neutron counting rate 
recorded on the moving film as a function of 
time “‘t’’ after the end of the exposure. When the 
activity had decreased to an acceptable value, 
the ion chamber was moved to the next radial 
position and another exposure made with the 
same sample. A single series usually comprised 
exposures in four to eight positions. A complete 


experiment required three complete series: (@) 


with heavy water in the sphere, in which case 
both delayed and photo-neutrons are present; 
(b) with the sphere empty, so that the activity 
is that caused by delayed neutrons alone; 
and (c) with the sphere empty and the fission 
chamber surrounded by cadmium, in order to 
measure the delayed neutrons with epicadmium 
energies. With heavy water in the sphere the 
number of epicadmium neutrons was negligible. 

The constancy of the pile power level and the 
detecting instruments was checked by repeating 
the initial exposure at the end of each expert 
ment. Instrument sensitivity alone was checked 
by means of a Ra- « -Be source. Variations in the 
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Fic. 3. Activity 
x R? vs. R curves wh 
at different times 
t after the end of 
a 30-sec. exposure. 
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pile power level were found to be less than one 
percent over the period during which a set of 
data was taken, so that the use of an exposure 
monitor was not necessary. 

For each exposure, a plot was made from the 
film data of the counting rate in the particular 
radial position as a function of ¢. From the group 
of curves forming a series, other curves were con- 
structed of (counting rateXR?) vs. R, for suc- 
cessive fixed values of ¢. Here, R is the radial 
distance of the detector from the source at the 
center of the sphere. The source activity at any 
time is then proportional to the area under the 
corresponding curve. 


TABLE I.* Delayed neutrons from U**, r=mean life, 
7,=half-life, E=energy in Mev, A=saturated exposure 
yield relative to the yield of the 22-sec. period. 











% tT E(Mev) A 
55.6 80.2 0.25 0.153 
22.0 31.8 0.56 1.00 

4.51 6.50 0.43 1.28 

1.52 2.19 0.62 1.45 

0.43 0.62 0.42 0.51 
0.05 0.07 0.15 


Total 4.54 








* Data for the first five periods from unpublished work of D. J. 
Hughes, J. Dabbs, and A. Cahn; for the shortest period from D. Hall. 











Examples of these (counting rateXR?) vs. R 
curves are shown in Fig. 3. With heavy water in 
the sphere, the curves on a semi-log scale are 
linear with a relaxation length (inverse slope) 
which is constant and equal to 1.25 inches. This 
indicates that the neutrons have been well 
thermalized before they escape into the sur- 
rounding oil. When the sphere is empty, the 
semi-log plot of the neutron distribution shows a 
decided initial curvature caused by slowing down 
in the oil. 

In general, the desired source activity due to 
photo-neutrons alone, at a time /;, is the differ- 
ence between the areas under the corresponding 
distribution curves, with and without heavy 
water, after correction for epicadmium activity. 
By computing a number of points, a curve may 
be drawn of source photo-neutron activity as a 
function of time. . 

We assume that the data can be represented 
by a formula of the form (1), where A, is the 
photo-neutron activity of the source as a func- 
tion of time ¢ after the end of the exposure, 


A,(t)=La(1—exp—[7T/r7;]) exp—[t/r], (1) 


T is the exposure time, 7; the mean life, and a,;, the 
yield, is the activity of the ith radioactive com- 
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Fic. 4. 3.75-hour 
exposure with 
shows u 
data, whicharede. 
composed into 
curves 2, 3 and4 
of half-lives, 53 
hours, 4.4 hours, 
and 1.6 hours, 
respectively, 





4 7 4*1.65 bes. 

















No2, Pf *27.3min. 





Fic. 5. Thirty- 
minute exposure 
with DO. Curve 
No. 1 shows un- 
analyzed data. 
Curve No. 2 of 
half-life 27 min- 
ry oe 
y subtrac 
from curve 1 the 
activities caused 
by the 53-hour, 
4.4-hour, and 1.6 
hour periods. 








3 
TimE “t" IN HOURS 


ponent for T=, t=0. The longest period is 
determined first, at large values of ¢ when all 
shorter periods have died out, and its activity 
subtracted from the total. The shorter periods 
are determined, successively, by a repetition of 
the process. A similar formula for A,(t), the de- 


layed neutron activity, was fitted to the decay 
curve taken with the sphere empty. In this case, 
the yields and periods* were taken from Table I. 


* The yields and periods of the delayed neutrons have 

been studied by several groups of workers on the uranium 
roject. One of the earliest studies is that of A. H. 

f. S. Levinger, R. G. Wilkinson, E. P. Meiners, and M.B. 
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Fic. 6. Five-min- 
ute exposure. 
Curve No. 1 

shows unanalyzed 

data. Curve No. 2 
of half-life 7.7 © 
minutes is ob- 

tained from curve 8 
No. 1 by sub- a 

tracting the activ- ; 

ities caused by the 

4 longer periods. F 











No.2, The 7.7 min 



























Fic. 7. Three- 
minute exposure. 
Curve No. 1 is 
the total activity 
with DO. Curve 
No. 2 is obtained 
from curve No. 1 
by subtracting the 
elayed neutron 
activity. Curve 
No. 3, showing 2.4 
minute half-life, 
is obtained from 
curve No. 2 by 
subtracting the 
activities caused 
by the 5 longer 
periods. 








By extrapolation to t=0, T=, we can then 
express the yields a; of each photo-neutron period 
as fractions of the total delayed neutron yield 
Adti=0, T= @). 

Sampson, Phys. Rev. 70, 111 (1946). The values in Table I 


were taken from more recent unpublished work of D. % 
taghes, }. Dabbs, and A. Cahn, and also from the work 
- Hall. , 
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IV. EXPERIMENTAL RESULTS 





A. Periods and Uncorrected Yields 





The longest delayed neutron component has 
a half-life of 55 seconds, so for t >10 min. the 
delayed neutron activity is negligible in com- 
parison with that of the longer period photo- 
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Fig. 8, Thirty. 
€XPosure 


second 

ream in 
urve No, 

total acti 1 tthe 


is 
obtained from 
No. 3 by subtract. 
ing the activities 
due to the6 longer 
periods, 











neutrons. Under these conditions, decay curves ponents are 53 hours, 4.4 hours, 1.6 hours, 27 
of the gross activity with D,O in the sphere minutes, and 7.7 minutes. The longest period, 
represent photo-neutrons alone. Figures 4-6 show 74=53 hr., was actually followed somewhat 
the gross activity from ¢=15 min. to ¢=60 hr., beyond the 60-hr. limit of Fig. 4, but its half-life 
and the 5 component periods into which the data. is considerably less certain than those of the 
may be analyzed. The half-lives of these com- shorter periods. 


TABLE II. Photo-neutron yields and periods. 
74 = half-life, r= mean life. 
yi = uncorrected photo-neutron yield relative to that of the 22-sec. delayed neutron period. 
E,=v+-ray energy in Mev, from unpublished work of Hughes, Spatz, and Cahn. 


FF, = factor to correct yields measured with the 10 in. radius sphere of D,O to an infinite sphere. 
7 
@= correction for y-ray absorption in the rabbit wall, Al tube, etc. 

8=correction for multiple Compton scattering. 
A.=A Xa$/F,, yield corrected to infinite volume of D,O. 
Ag= 100 A./4.54, percent yield of photo-neutrons relative to the total yield of delayed neutrons. 
Ay=2X0.01 Az/100, absolute yield of photo-neutrons per fission using the value: fast neutrons/fission = 2. 


—-=number of y-ray photons per photo-neutron in an infinite volume of D,O, assuming each collision reduces energy 


; _ below threshold. 
C=correction factor allowing for photons not reduced below threshold in first impact. 


A. =absolute number of photons per fission. 


r 1 A 
Fy ( 1/Fy 


tT A Ey B Ae Aa 10+ 


76hr. 0,00038 (3) 
6.3hr. 00119 (3) 
24hr. 0093 2.62 

39 min. 0082 _2.60 

6 
2 
4 








>) 


Ref. 





1.085 0.00074 0.0163 0.0326 476 
1.085  .00232 .0515 101 476 
2 1.052 .0168 .370 742 670 
2 1.056 .0149 328 656 643 
11 min. 0124 3 1.085 .0242 534 1.065 476 
34min. .0276 2 1.054 .0504 1.11 2.23 650 
59 sec. .090 2 1.01 .147 3.24 6.50 2450 
3.6 sec. .225 3 ‘ ‘ 1.11 469 10.31 20.7 389 
Total—16% 3.2X10-% 
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Fic. 9. Thirty- 
second exposure 
volume _ integrals. 
No. 1 is the total 
activity with D,O. 
No. 2 is the total 
activity without 
DO. No. 3 is the 
activity caused by 
photoneutrons 
alone. No. 4, show- 
ing a half-life of 
2.5 seconds is ob- 
tained from No.3 
by subtracting 
the activities due 
to the 7 longer 

periods. 
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Figure 7, curve No. 1, shows the total activity 
with D,O following a 3-min. exposure. Curve 
No. 2 shows the photo-neutron activity alone, 
after subtracting from No. 1 the contribution of 
the delayed neutrons, determined with the 
sphere empty. A formula representing the 5 
previously determined photo-neutron periods was 
fitted to curve No. 2 at t=53 min. and then 
subtracted off, giving curve No. 3. A sixth 
period is indicated, 7,;=2.4 min. Extrapolating 
all yields back to t=0 and infinite exposure, the 
yields of the six photo-neutron periods could 
then be expressed in terms of the delayed neu- 
tron yields. These values are given in Table II. 
In Fig. 8, curve No. 1 shows the total activity 
with D,O and No. 2 the activity without. D,O, 
for an exposure of 30 sec. and a 0.15-g sample of 
enriched material. The delayed neutron curve 
has been corrected for epicadmium activity, 
which amounted to about 2 percent. The dif- 
ference between curves 1 and 2 amounts to 
only 10 percent of the total activity at t= 30 sec., 
so the experimental accuracy in determining the 
photo-neutron yields for the short periods is con- 
siderably less than for the longer periods. After 
subtracting off the longer periods, a component 
of half-life 41 sec. is indicated by curve No. 4. 
Figure 9 gives the result of a 30-sec. exposure 








i i i 1 jl 
7 i) 7 + T 1 








with a 1.5-g sample of unenriched uranium. 
Here the decay curve is extended back to within 
0.5 sec. after the end of the exposure. A period 
7,4=2.5 sec. is found. 


B. Estimated Accuracy 


The accuracy to be expected in obtaining the 
area under one of the curves of the type illus- 
trated in Fig. 3 is not better than one percent, 
owing to statistical probable errors alone. Hence, 
the yields given. in column A, Table II, for the 
two shortest periods where the photo-neutron 
activity is of the order of 10 percent of the total 
are uncertain by about 15 percent. Except for 
the 53-hr. period, the yields of the longer periods 
should be considerably more accurate. 

The half-life and yield of any short period 
depends on the previously determined constants 
of all longer periods, giving a cumulative error. 
The possibility cannot be ruled out that some of 
our shorter periods are fictitious, that perhaps 
some of them cannot be identified with individual 
fission isotopes. Indeed, the suggestion is made 
below that the shortest periods given in Table II 
cannot represent a single isotope but must be 
due to many half-lives too close together to be 
resolved. However, the periods and yields given 















do represent the photo-neutron activity as a 
function of time fairly well. 

If the heavy water were as pure as stated by 
the manufacturer, the absorption of neutrons in 
it would have been negligible. Since contamina- 
tion might have been introduced during the 
course of the experiment, the neutron absorption 
of the heavy water was tested in the following 
manner: An Sb-Be photo-neutron source was 
placed at the center of the sphere and the neutron 
distribution curves in the oil were obtained, with 
and without heavy water. The y-rays from Sb 
have energies less than 2.17 Mev and cannot 
create any photo-neutrons in the heavy water, 
so the areas under the (counting rate) X R? vs. R 
curves should be equal in the two cases. Actually 
the areas were found to be the same within about 
2 percent. 





C. Corrections to Experimental Yields 


In Table II, column A represents photo- 
neutron yields (saturated activities) for our 
sphere of heavy water relative to the 22-sec. 
delayed neutron period yield as 1.00. In order 


large sphere of D,O and complete y-ray absorp- 

tion, it is necessary to know the energy of the 
, ‘y-rays. Conclusions about the initial energy of 
\Y the photo-neutrons could not be made from the 
data of this experiment. The large amount of 
heavy water moderator slowed the neutrons 
down to such an extent that the shape of the 
distribution curve in the oil, where the measure- 
ments were taken, was no longer a measure of 
the initial neutron energy. Hence, we have tried 
to correlate our periods with those found by 
W. D. B. Spatz, D. J. Hughes, and A. Cahn in 
unpublished work on the same subject. Their 
values for the y-ray energies are given in column 
E,, Table II. For purposes of calculation we have 
arbitrarily assumed an energy of 3 Mev for the 
two longest periods, which were not reported by 
Spatz, Hughes, and Cahn. 

In D.O and for y-ray energies between 2 and 
5 Mev the photoelectric effect and pair produc- 
tion are negligible in comparison with the 
Compton effect. If we assume that every Comp- 
ton scattering of a y-ray quantum by an electron 
reduces the quantum energy below the y,n 
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{ to correct these yields to the case of an infinitely . 





threshold in D,O, the fraction of the quanta 
absorbed in a sphere of radius R cm, from j 
central point source, is F;= 1—exp — C(us+m,)R] 
where y, is the linear “absorption” coefficient per 
cm for Compton scattering and yy is the coef. 


ficient for the y, reaction in D.O. For this ex. 


periment the effective thickness of heavy water 
is 23.3 cm. us=o,X3.35 X 1078 per cm? for D.O, 
where o,, the Compton scattering cross section 
per electron varies from about 1.27 10-5 om: 
at 2.5 Mev to 1.06X10-** at 3.5 Mev. u,=67 
X10" Xo, per cm’, where o,, the y,m cross gee. 
tion per deuteron varies from zero at 2.17 Mey 
to about 14X10-* cm? at 3.5 Mev. In Table J] 
are listed values of 1/F; the correction factor by 
which the yields A must be multiplied to give 
the yields for an infinite sphere of D,0O, for the 
assumed energies E,. 

A small correction, a=1.06, is listed in the 
next column of Table II. This compensates for 
the y-ray absorption in the 0.050-inch wal] 
thickness of the pneumatic tube, the fixed 
0.125-inch thick Al tube through the center of 
the sphere, and in the ;s-inch thick wall of the 
plastic rabbit. 

Column £8 in Table II gives an additional cor- 
rection necessary because not every Compton 
scattering process reduces the energy of'a y-ray 
quantum below the y,m threshold, so that the 
flux of y-rays of energy above the threshold is 
everywhere higher than that calculated using the 
full Compton scattering cross section. 8/F; is 
then the ratio of the number of y,” neutrons 
produced by multiple scattered photons in an 
infinite sphere of D,O to that produced by 
multiple scattered photons in a sphere of radius 
23.3 cm. The theory of the multiple scattering 
correction was developed for us by Mr, E. 
Greuling and Mr. H. Soodak. 

Column A, in Table II is A Xa8/F;, the cor- 
rected yield of photo-neutrons in an infinite D,O 
medium, relative to the yield of* the 22-sec. 
delayed neutron period. Ag= 100 XA./4.54 is the 
percent yield relative to the total delayed neu- 
trons (see Table I, column A). Column As=2 
X0.01XA./100 is the absolute yield of photo- 
neutrons per fission, calculated from the data in 
column B, Table I, using the arbitrary value 2 
for the number of instantaneous neutrons 
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emitted per fission of a U5 nucleus and 0.01 for 
the ratio‘ of delayed to instantaneous neutrons. 


D. Gamma-Ray Yields 


We now compute the number of y-rays per 
fission necessary to produce the photo-neutrons 
of each period, with yields A; in Table II. The 
fraction of the total y-ray quanta from a point 
source which are absorbed in the y,” reaction in 
an infinite D,O sphere is 


Fy =y/ (Hythe), 


. provided each scattering process is assumed to 
reduce the energy below the threshold. If we 
take into account the fact that some scattered 
quanta still have energy above the threshold, 
the actual photo-neutron production will be 
higher by a factor C, which will be small near 
the threshold and increase with the initial energy 
of the y-ray. In Table II are listed, for each 
period, the values of 1/F,, and of C. The absolute 
yield of y-ray quanta per fission, A,, will be 


A,=Ay;/F,C. 


V. DISCUSSION OF RESULTS 


A. The Effect of Possible Parent-Daughter 
Contributions 


In the analysis of the data as described above 
it has been assumed that the hard y-ray emitters 
are initial fission fragments. The interpretation 
of the data is made uncertain because of the 
possibility that the emitter of the photon giving 
the photo-neutron may be a daughter of the 
initial fragment. In such cases the data would 
yield the true period, but the amplitude deduced 
by extrapolating the decay back to zero time 


‘H. D. Smyth, Atomic Energy for Military wie ag 


(Princeton University Press, Princeton, 1945), Appendix 3; 


or Rev. Mod. Phys. 17, 459 (1945). 
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may be in error. If the daughter has a very much 
longer life than its parent, the error will not be 
large. If the two have comparable lifetimes the 
saturated activities given may be a considerable 
overestimate. 


B. Number of Photons per Fission 


From the number of photo-neutrons per fission, 
the number of hard photons per fission was cal- 
culated from a consideration of the relative 
probabilities of Compton and y, collisions. The 
results are given in column A, of Table II. The 
sum of the values in this column is 2.5 photons 
per fission of energy above 2.17 Mev. The 
highest single yield is 1.6 for the 41-sec. period. 
This value, however, is subject to great uncer- 
tainty. The assumed energy E,=2.25 Mev is 
very near the threshold, where the y,” cross 
section is small and changing rapidly. If the 
true energy is only slightly higher, the calculated 
y-ray yield A , would be much reduced. However, 
the highest observed yields of individual fission 
products are only about 6 percent, so the y-ray 
yields of the 2.5-sec. 41-sec., and 2.4-min. 
periods are all much too high to be caused by 
single isotopes. These high photon yields can be 
explained by supposing that there are a large 
number of high energy y-emitters with half-lives 
too close together to be resolved in this experi- 
ment. 


C. Correlation of the Gamma*Ray Periods with 
Known Fission Products 


We have been unable to correlate without 
ambiguity our periods with those given for the 
known fission products. It seems that in order to 
accomplish this correlation it would be necessary 
to repeat the experiment with individual fission 
products or with groups of products each con- 
taining a small number of members. 
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A Calorimetric Determination of the Energy Produced by Plutonium (239) 
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The rate of energy production by alpha-disintegration 
of plutonium (239) has been measured calorimetrically. 
The method was to measure the rate of evaporation of 
liquid nitrogen by a 120-gram sphere of metallic Pu and 
compare this rate with that produced by the introduction 


of electrical energy. The power produced by Pym r 
1.923+0.019 abs. watts-gram™. Combining this With the 
energy per alpha-particle determined from 
range-energy measurements, one obtains 2.411 x 1@ Years 
as the half-life of Pu® with respect to alpha emission, 





HE energy produced by the alpha-particles 
from plutonium (239) is sufficient to cause 
a sample of the metal to feel slightly warm to 
the touch. There was available at the Los 
Alamos laboratory a very pure sample of Pu®® 
in the form of a 120-gram sphere. The heat 
generated by this sphere was about 3 calories 
per minute. This paper describes a calorimetric 
experiment designed to measure the heat to an 
accuracy of about 1 percent. 

The power generated per atom of Pu™® can be 
combined with the value of the energy per alpha- 
particle to give the number of alpha-particles 
per unit time and consequently the half-life of 
Pu*® with respect to alpha-disintegration. In 
the sample we used, all sources of energy other 
than the alpha-disintegration were extremely 
small in comparison with the experimental error. 
The calorimetric value of the power, the meas- 
urements of which are described below, gives an 
accurate determination of the Pu*® half-life 
which is independent of the methods ordinarily 
used. 

The sphere of Pu®® was 2.25 cm o.d. and had 
a central spherical hole 0.8 cm in diameter. It 
was made in two hemispheres. The surface, 
including the hemispherical hole and the equa- 
torial planes, was coated with a 0.0025-inch film 
of silver. The sphere with silver coating weighed 
121.53 grams and contained 120.06 grams of 
Pu®®, 

The method adopted was to immerse the 
sphere in a calorimeter containing liquid nitrogen 
and measure the rate of evaporation of nitrogen 
gas due to the heat generated. Calibration experi- 
ments were done in which electrical energy was 


* Present address: Institute for the Study of Metals, 
University of Chicago, Chicago 37, Illinois. 


introduced at a rate approximating that of the 
sphere. 


DESCRIPTION OF THE APPARATUS 


The calorimeter was a silvered Pyrex glag 
Dewar vessel. It is shown in Fig. 1. Above the 
Dewar is a ground glass joint which was sealed 
with rubber cement diluted with toluene. This 
type ‘‘of seal’ has been previously described! a5 
suitable for use at low temperature, and we found 
it satisfactory. In every case the joint was tested 
and found to be tight after being cooled tj 
liquid-air temperature and before beginning 
measurements. At the bottom of the calorimeter 
was an electrical heater, D, made of 275 ft. of 
No. 40 D.S.C. copper wire wound on a cylinder 
1% inches in diameter and 2? inches high. Two 
leads of No. 30 copper wire, C, were soldered to 
the inside of the tungsten wires, B, sealed through 
the glass. Attached to the outside of each tungsten 
wire was a current and potential lead of No. 0 
copper wire. These go to the electrical system. 

The electrical energy was supplied by lead 
cells. The current was measured by reading with 
a 100,000 microvolt Leeds and Northrup White 
Double Potentiometer and H. S. Galvanometer 
the potential drop across a calibrated standan 
one-ohm resistance and the voltage by meansd 
the same potentiometer used with a calibrated 
volt-box constructed from standard 100 and 
10,000-ohm resistances. The accuracy of th 
electrical measurements is about 0.02 percent, 
and their contribution to the over-all error d 
the measurements is negligible. Corrections wet 
made in the usual manner? for the energy dist 
- 1J. W. Stout and W. F. Giauque, J. Am. Chem. Se 
63, 714 (1941). 


‘2W. F. Giauque and R. Wiebe, J. Am. Chem. Soc. & 
101 (1928). 
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pated in the volt-box and in the single leads to 
the heater. 

The sphere was supported in a cylindrical 
cradle, E, }§-inch diameter made of brass foil 
and closed at the bottom with a copper screen. 
The cradle containing the two hemispheres, F, 
each with its equatorial plane upwards, was 
lowered to the bottom of the calorimeter, and 
later removed, with a heavy thread, G. Since 
the measurements are made at constant temper- 
ature, the heat capacities‘ of the sphere and 
cradle which are, of course, not present during 
the calibration runs, do not introduce any error. 

The calorimeter was suspended by metal sup- 
ports inside a large glass Dewar vessel 5} inches 
id. and 54 inches long. A glass filling tube, 0.6 
cm i.d., extended through the closure of the 
Dewar vessel and connected the top of the calo- 
rimeter to stopcocks for evacuation, filling with 
nitrogen, and metering of the nitrogen evapo- 
rated. The tube was bent in a spiral a few inches 
above the calorimeter in order to intercept room 
temperature radiation. The large Dewar vessel 
was filled with liquid nitrogen-oxygen mixture 
of such composition as to be 2° to 5°K warmer 
than the calorimeter: 

In making a measurement, the calorimeter was 
assembled (with or without the sphere), and the 
large Dewar containing liquid nitrogen-oxygen 
mixture was slowly raised in place. More liquid 
nitrogen was added to the large Dewar through a 
tube in the top closure so as to bring the liquid 
level well above the bend in the glass tube. The 
large Dewar was then evacuated until the temper- 
ature of the liquid in it fell below the boiling 
point of nitrogen at local atmospheric pressure. 
Meanwhile, the calorimeter was evacuated and 
flushed with dry nitrogen from a cylinder. All 
nitrogen condensed in the calorimeter was from 
this cylinder, which contained 99.8 percent pure 
N3. It was dried by passing through two drying 
tubes filled with magnesium perchlorate and then 
went through a gas meter. The nitrogen passed 
through the meter at atmospheric pressure and 
then was condensed in the glass filling tube, 
from which the liquid fell into the calorimeter. 
About two hours were required to condense 110 
cc of liquid nitrogen in the Dewar during the 
sphere measurements. When sufficient liquid was 
condensed in the calorimeter, dry air, nitrogen 
or oxygen was bubbled through a tube which 
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extended to the bottom of the large Dewar 
vessel, until the liquid in this vessel was at 
its boiling point and uniform in temperature 
throughout. The rate of evaporation of nitrogen 
from the calorimeter due to the heat of the 
sphere or to electrical energy was measured with 
a second gas meter. A period of several hours 
after condensing the nitrogen was allowed in 
order for the thermal system to come to equi- 
librium. Since all gas entering or leaving the 
calorimeter was measured, the amount of liquid 
nitrogen in the calorimeter was known at all times. 

The temperature of the bath was measured 
with three copper-constantan thermocouples at- 
tached to the outside of the calorimeter at points 
marked H in Fig. 1. They were read on the 
White potentiometer with an accuracy of about 
0.1°K. The thermocouples were calibrated in a 
bath of liquid oxygen. Fhe temperature of the 
liquid in the calorimeter was calculated from its 
pressure and the vapor pressure equation of 
Giauque and Clayton.* 
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*W. F. Giauque and J. O. Clayton, J. Am. Chem. Soc. 


55, 4875 (1933). 
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The meters used were Precision wet test 
meters, which had been filled with butyl phthal- 
ate to eliminate the need for correcting for the 
vapor pressure of the sealing fluid. The meter 
hand made one revolution for every 0.1 cubic 
foot of gas. The time for a complete revolution 
was measured with a calibrated stopwatch. If 
care is taken to time only complete revolutions, 
the meter is reproducible to about +0.5 percent. 
In the sphere measurements about 36 minutes 
were required for one revolution of the meter 
hand (0.1 cubic foot). 

The temperature of the gas entering the meter 
was measured with a mercury thermometer, and 
the pressure was calculated from the barometric 
pressure combined with the pressure of a butyl 
phthalate manometer which measured the pres- 
sure difference between the gas entering the 
meter and the atmosphere. All pressure readings 
have been converted to international centimeters 
of mercury, using 979.1 cm sec.~ as the local 
acceleration of gravity and 980.665 cm sec.~ as 
the standard acceleration. 


HEAT LEAK CALIBRATIONS 
Measurements were made of the rate of evolu- 
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tion of nitrogen when no heat other than that 
due to heat leak was being introduced inty 
the calorimeter. These measurements 
throughout the period of calibration and 
runs and show no significant change with time 
At first the gas evolved was measured with the 
meter, but since the rate was only a few percent 
of that during the calibration or sphere 
the accuracy was not great. Later the rate of gas 
evolution during heat leak measurements Was 
determined by allowing the gas to collect in an 
evacuated rubber prophylactic sack. The volum 
of gas was then determined by measuring with 
a buret the equivalent volume of water. 

The heat leak was found to be a function g 
the height of liquid in the calorimeter. Most ¢ 
the leak is due to conduction by residual gas 
across the evacuated space between the walls gf 
the calorimeter Dewar. Measurements were mage 
at different temperature differences between bath 
and calorimeter and at different heights of liquid 
in the calorimeter so as to cover the range in the 
calibration and sphere experiments. It was foun 
that the heat leak obeyed Newton's law, 


dn/dt = B(T path 7 T calorimeter); 


where dn/dt is moles-min.— of nitrogen evapo 
rated; B is coefficient of heat leak, moles-min7 
-deg.!°K; T is temperature in °K. The valug 
of B are plotted as a function of height of liquid 
in Fig. 2. This curve was used in calculating th 
calibration and sphere measurements. Points 
indicated by open circles were taken with the 
rubber sack and are given more weight tha 
those indicated by filled circles which wer 
measured with the gas meter. 


TABLE I. Electrical calibration measurements. 








Corrected 
rate of 
evolution 
of nitrogen 
moles min.~! 
X10 


Heat leak 


Obs. flow e 
correction, 


aa Tbath 


min.~! —Teal.. 
x10 °K 


2.186 4.41 
2.487 4.51 
2.722 4.51 
3.053 4.55 
2.552 2.26 
2.574 5.48 
2.501 5.43 
2.770 2.39 
2.653 2:51 


Average Average 
temp. 


pressure 
of nitrogen, of nitrogen, 

int. cm at meter, 

mercury °K 


576.5 296.3 
576.7 295.2 
$77.3 294.7 
577.6 293.8 
579.8 295.9 
575.9 295.4 
575.6 295.5 
582.0 291.5 
582.7 293.5 


Electrical 


power, 
int. 
watts 


0.19891 
0.22672 
0.24869 
0.27878 
0.23540 
0.23898 
0.23188 
0.25239 
0.24248 


Duration 
of obser- 
vations, 


Z 
9 





2.119 
2.422 
2.659 
2.992 
2.514 
2.498 
2.428 
2.727 
2.615 
Mean 
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TABLE II. Measurements of heat evolved by Pu sphere. 

















Corrected Heat to 
Average Average Obs. flow Heat leak rate of vaporize 
Duration pressure temp. of nitrogen, r correction, evolution nitrogen 
of obser- of nitrogen, of nitrogen moles bath Av. ht. moles of nitrogen int. Heat evolved 
vations, int. cm at meter, min.~! —Teal., of liquid, min.~! moles min.~} joules by sphere, 
No. min. mercury “— x 108 °K cm x<108 x10" mole=! int. watts 














144 584.7 294.9 2.501 2.01 


2 250 585.3 294.8 2.515 5.16 
3 216 584.5 294.4 2.511 2.93 











11.47 0.031 2.470 5629 0.2317 
9.08 0.070 2.445 5629 0.2294 
11.18 0.044 2.467 5629 0.2314 
Mean 0.2308+-0.0010 
(av. dev.) 











F CALIBRATION MEASUREMENTS 


The data for the electrical calibration meas- 
urements are given in Table I. These measure- 
ments were taken both before and after those on 
the sphere. The agreement is good except for 
runs No. 6 and 7 which are unaccountably high. 
Since we have not been able to find any error in 
these two measurements, they have been in- 
cluded in the average. The top of the heating 
coil was 7.2 cm above the bottom of the Dewar 
so it was always covered with liquid. 

The average value of the heat to vaporize 
nitrogen is 5629+48 (av. dev. from mean) int. 
joules-mole“ or 1346+11 cal.-mole~. (One cal. 
=4,1833 int. joules.) It is interesting to compare 
this with the more accurate heat of vaporization 
of nitrogen measured by Giauque and Clayton.’ 
They obtained 1332.9+1.0 cal.-mole at the 
normal boiling point, 77.32°K. Taking AC, =6.6 
cal.-mole— one obtains for the heat of vaporiza- 
tion of nitrogen at 75.17°K, the temperature 
corresponding to our measurements, the value 
1347.1+1.0 cal.-mole—'. The agreement of our 
value with this is evidence that there are no large 
undetected errors in our measurements. 


THE HEAT PRODUCED BY THE Pu** SPHERE 


The data on the runs with the Pu*®® sphere 
are given in Table II. The height of the topmost 


‘portion of the sphere was about 3.8 cm; so it 


was well covered with liquid at alltimes. The 
average value of the heat produced is 0.2308 int. 
watt or 3.310 cal.-min.—'. It is believed that this 
value is correct to within +1 percent. 


THE HALF-LIFE OF Pu*®* WITH RESPECT TO 
ALPHA-DISINTEGRATION 


The weight of Pu®® in the sphere was 120.06 g. 
The atomic weight of Pu®® is taken as 239.08. 





Other physical constants are from Birge’s 1941 
tables.‘ The energy given out by Pu®® is, there- 
fore, (0.2308 X 1.00020) /120.06 = 1.923 K 10 abs. 
joule-sec.'-gram™, or 7.633X10—* erg-sec.— 
-atom™, Possible sources of energy due to other 
nuclear reactions and to impurities have been 
considered. These are*small in comparison with 
the experimental error of our measurements, and 
may be neglected. 

The range of Pu®* alpha-particles has been 
measured by Chamberlain, Gofman, Segré, and 
Wahl as 3.675 cm relative to polonium as 3.842 
cm. Taking the energy of the polonium alpha- 
particles (5.2984 Mev) and the range-energy 
relationship from the paper of Holloway and 
Livingston® one obtains 5.144 Mev as the energy 
of the plutonium alpha-particle. The recoil of 
the residual nucleus contributes 0.086 Mev so 
the total energy per alpha-disintegration of Pu™® 
is 5.23 Mev or 8.38 X10~—* erg-atom—. 

The half life of Pu®* for alpha-emission is then 


T, = (0.69315 X 8.38 X 10-*) /(7.633 K 10-8 
X 3.15569 X 10") = 2.411 X 10* years. 


We wish to thank Professor E. Segré for many 
helpful discussions about this problem. 

This paper is based on work performed under 
Contract No. W-7405-Eng-36 with the Man- 
hattan Project at the Los Alamos Scientific 
Laboratory of the University of California. The 
information contained in it will appear in Divi- 
sion V of the Manhattan Project Technical 
Series as part of the contribution of the Los 
Alamos Laboratory 


*R. T. Birge, Rev. Mod. Phys. 13, 233 (1941). 

* Chamberlain, Gofman, Segré, and Wahl, Phys. Rev. 
71, 529 (1947). 

* Holloway and Livingston, Phys. Rev. 54, 18 (1938). 
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From the beta-activity induced in beryllium by (1) fast deuteron bombardment, and (2) 
slow neutron capture, the half-life of Be” is determined. The value is between 10* and 10? 
years. The beta-end-point is roughly 650 kev, and no gammas were observed. These results 
confirm those of McMillan and Ruben. The B"(n,p)Be” cross section for thermal neutrons 
is computed by comparison with the known value for B"(n,a)Li’, assuming disintegration 
to the ground state possible for the former reaction. 





INTRODUCTION 


— and Ruben! have recently pub- 
lished measurements made in 1940 on the 
half-life and beta-end-point of a long-lived beryl- 
lium activity. A mass-spectroscopic separation 
by Pierce and Brown? has established that the 
active isotope has mass 10. In 1944 we measured 
the activity produced by Be® (d,p) and (n,7) 
reactions; and from these yields estimated the 
Be’® half-life. We have also measured the beta 
end-point, and studied the gamma-activity. 
These measurements are not as complete as we 
would like, owing to our transfers to other 
sections of the Manhattan Project. 


CHEMICAL SEPARATIONS 


Our strongest sample of Be!® came from the 
thick beryllium target of the Washington Uni- 
versity cyclotron, which had been activated by 
about 600,000-microampere hours bombardment 
by 10.45-Mev deuterons.’ The target was etched 
to the calculated depth of deuteron penetration, 
giving us 1.7 grams of beryllium. The beryllium 
was purified by a chemical procedure similar to 
that used by McMillan and Ruben. The initial 
separation consisted of a chloroform extraction 
of beryllium basic acetate, and purification by 


t This document is based on work performed under 
Contract No. W-7401-eng-37 for the Manhattan Project 
and the information covered in this document will appear 
in Division IV of the Manhattan Project Technical Series, 
as part of the contribution of the Argonne National 
Laboratory. 

* Now at Cornell University, Ithaca, New York. 

** Now at Washington University, St. Loyis, Missouri. 
( 946) M. McMillan and S. Ruben, Phys. Rev. 70, 123 

1 ' 
(1946) K. Pierce and F. F. Brown, III, Phys. Rev. 70, 779 

3 Based on a recent measurement of the deuteron energy 
by Kurie. 


precipitation of the copper group of sulphide, 
Six different samples were then prepared fo 
counting the beta-activity, using different addi. 
tional purifications in preparing each sample, as 
follows: (1) five beryllium hydroxide precipita. 
tions, and precipitation of the zinc group of 
sulphides, to remove contaminants; (2) one 
beryllium hydroxide precipitation; (3) three 
beryllium hydroxide precipitations; (4) two 
beryllium hydroxide precipitations, and an jp. 
complete berryllium ammonium phosphate pre. 
cipitation ; (5) a beryllium ammonium phosphate 
precipitation of the material not precipitated in 
the previous procedure; and (6) a chloroform 
extraction of beryllium basic acetate. The aver 
age specific activity of these 6 samples was 
measured on a GM tube as 17.5+1 betas/sec. 
mg of beryllium. This small difference amiong 
the six samples established that the activity was 
beryllium. | 


HALF-LIFE OF Be” 


Let Y’ be the d-p yield for a thick beryllium 
target activated by 10.45-Mev deuterons. The 
total deuteron flux was 1.310”. From the 
specific activity given above, the activity pro 
duced in the whole target is calculated as 29,000 
betas/sec. This figure is certainly low, as we 
know that much of the beryllium was sputtered 
off the target during the bombardment.‘ The 
upper limit on the Be!® half-life is then 


1.3 X 10** X 0.693 Y’ 
<< 


= =1.010"°Y’ years. 
29,000 X 3.16 x 107 





‘ This spalling off of active material has recently beet 


confirmed by measurements of E. , 
other _—- University target (private commutt 
cation). 
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McMillan and Ruben measured the Be’® 
activity produced by 12,000-microampere hours 
bombardment of beryllium by 16-Mev deuterons. 
They give the half-life as 6.4 10*°Y years, where 
Y is the d-p yield for 16-Mev deuterons. For 
comparison with their results, we may estimate 

/Y' from the ratio of the deuteron ranges as 
about 2. Then our measurement gives a half-life 
eight times the McMillan and Ruben determi- 
nation; but, as already pointed out, our results 
can only be considered as an upper limit. 

Assuming that the d-p cross section for high 
energy deuterons is somewhat less than half the 
geometrical cross section, and that the effective 
nuclear radius is about 5 X 10-" cm, we calculate 
a d-p cross section of 0.3 barn (+0.3 b., —0.15 
b.). This agrees in order of magnitude with the 
cross section as estimated from the thick target 
d-n yield for 1-Mev deuterons.’ By use of this 
cross section, Y’ is calculated as one nucleus of 
Be! formed per 1000 deuterons, and Y as 1/500. 
The Be’® half-life from our results is estimated 
as less than 10’ years; and from McMillan and 
Ruben’s measurements as about 1.3 X 10° years. 

Since the beryllium cross section for capture 
of slow neutrons is known,® the Be"® half-life can 
be determined from the specific activity induced 
in beryllium by a known neutron activation. 
Beryllium metal was activated in the Clinton 
pile for several months. After chemical separa- 
tions similar to those given above for purification 
of the Washington University cyclotron target, 
the beta-activity was only 10 counts/min. We 
are uncertain whether this low activity is due to 
contamination, or to Be!®. The upper limit for 
the specific activity due to Be’® is taken as 0.02 
beta/sec. mg of beryllium; and the Be’® half-life 
is calculated as greater than 10° years. This value 
is in fortuitously good agreement with our calcu- 
lation from McMillan and Ruben’s measure- 
ments; and is consistent with the upper limit on 
the half-life set by our measurement of the d-p 
yield. 


$6 AND y ACTIVITY 
Deflection of the particles by a small electro- 


- magnet placed between the active beryllium from 


*H. A. Bethe and M. S. Livingston, Rev. Mod. Phys. 9, 
330 (1937). 

*M., Goldhaber, Phys. Rev. 70, 85 (1946); and work at 
the Metallurgical Laboratory, to be published in the 
Plutonium Project Reports. 
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the cyclotron target and the GM tube showed 
that the particles were negatrons. Absorption 
measurements in aluminum from 0.2 to 300 
mg/cm? were made for both a thin and a thick 
source. The thin source gave an inspection end- 
point of roughly 200 mg/cm’; and an end-point 
from Feather analysis of 225 mg/cm’, or 650 kev. 
Owing to the low beta-intensity and high gamma- 
background, the Feather analysis was carried 
only to 0.6 of the range, and the extrapolated 
range is therefore inaccurate. The inspection 
end-point for the thick source was 200+20 
mg/cm?. These measurements are in fair agree- 
ment with the end-point of 180+20 mg/cm’, or 
550 kev, given by McMillan and Ruben. Both 
measurements are consistent with the Be'®—B'° 
mass difference of 550+110 kev.’ 

Initially a gamma-ray of intensity roughly 4 
quanta/beta was observed in the active material, 
and chemically determined to be.a beryllium 
activity. The half-life from decay measurements, 
and energy from lead absorption for this gamma- 
activity were both in good agreement with the 
kriown values for Be’. No gamma-activity due 
to Be'® was observed, in agreement with the 
results of McMillan and Ruben. 


PRODUCTION OF Be” 


Samples with higher specific activity of Be'® 
might be produced by a Li’(a,p)Be!® reaction 
using cyclotron alphas; by a C(n,a) Be’ reaction 
with fast neutrons (since this reaction is endo- 
thermic by about 3.5 Mev); or by a B'°(m,p) Be'® 
reaction in a chain-reacting pile. Bretscher* was 
unable to find Be’® activity produced by this last 
reaction, using deuterons on lithium as a neutron 
source. He sets the upper limit for the fast neu- 
tron cross section for this reaction as 510-7; 
barn, where 7, is the Be’® half-life in years. 
Using our estimated half-life of 210° years, 
this sets the upper limit on the cross section as 
1 barn. Cross sections which are thousands of 
times smaller than this can be measured using 
the high neutron fluxes now available. For suffi- 
ciently high neutron energies we would expect 
the B'°(n,p)Be'® cross section to be of the same 

7B. J. Isaacs, Cornell Thesis, ‘‘A Revision of the Isotopic 
Mass le, (1942)"; H. A. Bethe and M. S. Livingston, 
Rev. Mod. Phys. 9, 245 (1937); and E. Pollard, Phys. Rev. 


57, 241 (1941). 
* E. Bretscher, Nature 146, 94 (1940). 
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TABLE I. Comparative data on B!°(",p)Be?® and B’°(n,a)Li’ 




















reactions. 
B(n,p)Be® 
McMillan and 
Measurements Ruben* This paper B(n,a)Li? 
beta-end-point Be!® 550 kev 650 kev 
reaction ene 230 kev 130 kev 2500 kev 
Coulomb barrier 1850 kev 1850kev 3700 kev 
penetration factor 5x 10-3 3x10- 0.45 
cross section (thermal 
neutrons) 8 barns 0.5 barn 700 barns 
* Reference 1. 


order of magnitude as the N'(n,p)C™ cross 
section, which is about 0.01 barn. 

The beta-decay theory’® predicts a large spin 
for Be!®; while we would expect a spin of zero, 
since this is an even-even nucleus. Measurements 
of the ratio of the B'°(n,p)Be!® and B'°(n,a)Li’ 
cross sections, especially for slow neutrons, might 
be of interest in establishing selection rules, 
leading to some knowledge of the spin of the 
ground state of Be!®. Comparison of the (n,p) 
cross section for thermal neutrons with that for 
fast neutrons might well be made by comparing 
the beryllium activity at the surface of a boron 
lump activated in a pile, with the beryllium 
activity in the inside of the lump. 


oH. H. Barschall and M. E. Battat, Phys. Rev. 70, 245 
1946). 
10 F, J. Konopinski, Rev. Mod. Phys. 15, 222 (1943). 
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For thermal neutrons we may estimate the 
B'°(n,p)Be!® cross section by comparing the 
penetration probability for the protons to escape 
through the Coulomb barrier, with the penetra. 
tion factor for the much more energetic alpha- 
particles emitted in the competing B'(n,q)Lj7 
disintegration, of known cross section." We as. 
sume here that the compound nucleus is per- 
mitted to disintegrate to the ground state of Bel. 
while the (”,a) reaction goes to a level of Li? 
about 0.5 Mev above ground.” Penetrations are 
calculated using the Gamow theory.” 

Table I shows that in the absence of selection 
rules the B'°(n,p) Be’® reaction could be observed 
by the beryllium activity produced by activation 
of boron by thermal neutrons in a pile; and that 
the cross section for this reaction is sensitive to 
the beta-end-point of Be’®. 
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The Transmission of Slow Neutrons through Microcrystalline Materials* 
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E. Fermi,** W. J. Sturm, AND R. G. Sacus 
Argonne National Laboratory, Chicago, Illinois 
(Received January 27, 1947) 


The transmission of monochromatic slow neutrons through microcrystalline Be and BeO 
has been determined. The source of neutrons was the Argonne heavy water pile. These neutrons 
were monocromatized by means of a mechanical velocity selector for low energies and a neutron 
crystal spectrometer for higher energies. The results are in excellent agreement with the 
theory of elastic scattering from crystals. It is found by comparison of the results on BeO 
with the theory that the scattering amplitudes of Be and O have the same sign. This method 
may be used to determine the relative scattering phases of other pairs of nuclei which can be 
combined to form a crystalline material. The sample must consist of crystals smaller than a 
micron in linear dimensions. Other possible sources of disagreement between theory and 


experiment are discussed in Section 5. 





1. INTRODUCTION 


HE measurement of the effective scattering 

cross section of crystalline materials for 
slow neutrons is of considerable interest in con- 
nection with problems of neutron diffusion. If the 
material consists of nuclei having a small ab- 
sorption cross section, the scattering cross 
section can be determined by measuring the 
transmission coefficient for slow neutrons. It is 
known that the scattering cross section deter- 
mined in this way is a very sensitive function of 
neutron energy which shows violent fluctuations 
for small changes in the energy. These fluctua- 
tions are associated with the appearance of 
Bragg reflections from the appropriate planes in 
those of the microcrystals which are properly 
oriented. The shape of the transmission curve 
has been determined theoretically by Halpern, 
Hammermesh, and Johnson! and by Weinstock? 
for materials containing one atomic species. The 
generalization to multi-atomic matetials is 
simple and will be given in Section 4. 

With the intense monochromatic neutron 
sources which are now available, it is possible 
to obtain a rather precise measurement of trans- 
mission curves as a function of energy. Measure- 
ments of this kind have been carried out for 

* This document is based on work performed under the 
auspices of the Manhattan Project at the Argonne Na- 
tional Laboratory, This paper was submitted for declas- 
sification on November 25° 1946. 

** Also a member of the Institute for Nuclear Studies, 
University of Chicago. 

+0. Halpern, M. Hammermesh, and M. H. Johnson, 


Phys. Rev. 59, 981 (1941). 
*R. Weinstock, Phys. Rev. 65, 1 (1944). 
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microcrystalline Be and BeO. It is the purpose 
of this paper to present the results of these 
measurements and to compare them with theory. 
In particular, it will be shown in Section 4 how 
the results can be used to determine the relative 
phases of the Be and O scattering. 

In Section 2, the experimental methods for 
obtaining monochromatic neutron beams are 
discussed and in Section 3 the experimental 
results obtained with these beams are’ given. 
Section 4 contains a discussion of the theory with 
particular reference to its application to ma- 
terials containing more than one atomic species. 
The limits of applicability of the theory are con- 
sidered in Section 5. 


2. EXPERIMENTAL METHOD 


The monoenergetic neutron beams required 
for the measurements were obtained by means 
of two velocity selectors in conjunction with the 
heavy water pile at the Argonne Laboratory. 
The first was the mechanical velocity selector 
first used by Fermi, Marshall, and Marshall as 
modified by Brill and Lichtenberger,* which by 
use of a motor driven cadmium shutter and 
proper timing devices can be applied to the 
measurement of cross sections for the range of 
neutron energies between 0.004 and 0.20 electron 
volt. The second, which extended the measure- 
ments in one case to about 1.0 electron volt, 

* This velocity selector will be described in forthcoming 
pose by E. Fermi, L. W. Marshall, and Z. Marshall and 

y T. Brill and H. D. Lichten . This description of a 


use of the instrument before publication of their papers 
is presented with the permission of these authors. 
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Fic. 1. The microcrystalline beryllium cross section. The 
solid line represents the cross-section function as calculated 
in Section 4. Triangles at the base of the curve indicate 
the resolution of the instrument in the various ranges of 
measurement. The cross sections, ¢, are given in units of 
10-* cm?/atom (barns). 


made use of the monoenergetic beams of neutrons 
diffracted from LiF (100) in the neutron crystal 
spectrometer which has been previously  de- 
scribed.‘ 

The bursts of neutrons periodically released 
by the rotating shutter of the mechanical ve- 
locity selector passed through a neutron detector 
after traversing a measured path between the 
shutter and detector. The output of the detector, 
recorded during short measured intervals after 
the release of the initial burst and during the 
time the burst was passing through the counter, 
made it possible to record only neutrons in a 
very small velocity range. By introducing a 
sample of a microcrystalline substance into the 
periodic beam released by the shutter, it was 
possible to measure the transmission of the 
substance as a function of energy by measuring 
the fraction of neutrons scattered by the sample 
in each small energy interval. The value of the 

_ effective. scattering cross section was calculated 
for this interval from the measured transmission 
data. 


4W. J. Sturm and S. H. Turkel, Phys. Rev. 70, 103 
(1946); W. H. Zinn, Phys. Rev. 70, 102 (1946). 


The neutron source was a beam from the 
graphite thermal column of the chain reactor 
collimated by a series of three absorbing slits tp 
a 2X3” area. The spectrum of neutrons emitted 
from this source had an approximate Maxwellian 
distribution of velocities with its peak at about 
0.04 electron volt. The intensities in the range 
between 0.004 to 0.20 electron volt were suf- 
ficient for these measurements; outside this 
range intensities were generally too low to be 
used. Passing through the shutter the inter. 
rupted beam fell upon the absorber and reached 
the detector, 150 cm from the rotating cylinder. 
The detector, mounted with its axis normal to 
the beam, was a BF; proportional counter, 38 
cm in diameter and 11 cm long, having a central 
wire 0.002 inch in diameter. The filling gas was 
enriched in the isotope B'° to increase its counting 
efficiency. In order to minimize the effect of 
background neutron radiation, the whole counter 
except for a port to admit the beam was syr- 
rounded by a 2.5-cm shield layer of boron 
carbide. 

Measurement of transmission was made by 
comparing the beam intensity at the various 
energies with the scattering sample interposed 
in the beam to the intensity of the unfiltered 
beam. A second series of similar measurements 
was taken at each energy to evaluate the inten- 
sity of neutrons leaking through the shield plus 
those fast neutrons which passed through the 
closed shutter. This measurement was made by 
stopping the shutter and setting it at the posi- 
tion which it assumed at the time the signal from 
the detector was sent to the recorders. In this 
position only fast neutrons could pass through 
the shutter to reach the detector. In all cases 
this background intensity was less than 5 percent 
of the total flux in each energy interval, reaching 
this fraction only in the extremes of the energy 
range. The transmission was then calculated 
from the relation 


T=R(a)/R(0), 


where R(a) is the counting rate in the filtered 
beam corrected for background, and R(0) is the 
counting rate of the open beam similarly cor 
rected. 
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3. RESULTS 
A. Beryllium 





The sample for the measurements was pre- 
from a block of very pure microcrystalline 
beryllium metal reduced to chips by turning on 
a lathe. These chips, subsequently further pul- 
yerized by trituration in a mortar and pestle 
made of beryllium metal, were reduced to their 
final size in a beryllium ball mill. A 3.25” 
diameter cylindrical disk of the material having 
a surface density of 2.45 grams/square centi- 
meter was prepared by pressing in a die. The 
cross section per atom shown plotted in Fig. 1 
was calculated from the relation 


o=—In7T/N 


in which N is the number of atoms of Be per 
cm? and T is the transmission measured in the 
manner discussed above. Solid lines represent 
theoretically predicted values. Uncertainties in 
the cross section were evaluated by assigning a 
statistical uncertainty equal to the square root of 
the number of counts to each of the four measured 
components of the transmission value. Because 
the intensities were lowest at the extremes of the 
energy range, cross-section measurements at 
these energies show the greatest statistical un- 
certainty. 


















B. Beryllium Oxide 


A sintered block of BeO, prepared from very 
fine powder crystals was used to measure the 
scattering cross section for this compound. Cross- 
section values were calculated from transmission 
measurements and, as can be seen from Fig. 2, 
agree very well with the theoretical treatment 
given below. 

In the case of the beryllium oxide curve, 
mechanical velocity selector measurements were 
repeated on the neutron crystal spectrometer in 
the range between 0.04 and 0.20 electron volt 
and extended to 1.0 electron volts by this 
method. 






4. THEORY 





Although the theory of the coherent scattering 
of neutrons by microcrystalline materials has 
been given,'? it seems worth while to indicate the 
physical factors which go into the results. For 
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this purpose, we consider a beam of neutrons 
moving in a given direction incident upon a 
single microcrystal. The wave-lengths of the 
neutrons in the beam are assumed to be dis- 
tributed uniformly between A and A+AA. The 
orientation of the microcrystal is such that one 
of the wave-lengths in the interval AX undergoes 
Bragg reflection from some particular set of 
lattice planes. Then, since the scattered wave is 
undergoing constructive interference, the scat- 
tered intensity is proportional to the square of 
the number of scattering centers. Thus, if No is 
the number of nuclei in the microcrystal and F 
the scattering cross section of a nucleus, the 
scattered intensity is proportional to N¢F. 
Since the crystal is finite in size, the scattered 
beam is spread by diffraction over a solid angle 
of the order of (A/1)?/cos@ where | is the average 
linear dimension of ‘a microcrystal and @ is the 
angle between the scattered beam and the 
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Fic, 2. Scattering cross section of BeO. Alternate theo- 
retical curves for same and opposite phase of neutron 
scattering are shown as solid and dotted lines, respec- 
tively. Results indicate that the Be phase is the same as 
the oxygen phase. 
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normal to the crystal plane responsible for Bragg 
reflection. Since cos@ = b\/2 where b is the product 
of the order of the reflection and the reciprocal 
of the spacing between lattice planes, the inten- 
sity of the scattered beam integrated over its 
spatial width is proportional to 


\NeF/Pb. (1) 


Only neutrons within a small range of wave- 
lengths in the interval AX contribute to this 
scattering. The fraction of the neutrons with the 
proper wave-length can be obtained from the 
resolving power of the crystal which is d\/A = 1/01 
where b/ is equal to the product of the order of 
the reflection and the number of crystal planes 
contributing to the reflection. The fraction of the 
incident neutrons which are scattered is therefore 
d\/AX\=/bIAX, so the ratio of the intensity 


scattered to the incident intensity is propor- 


tional to 
NF /Pb AX. (2) 


This scattered intensity is now to be averaged 
over all orientations of the microcrystal. Since 
the crystal gives rise to a Bragg reflection only 
if its orientation is such that the Bragg condition 
is satisfied for some wave-length within AX, a 
contribution to the scattering is obtained only 
in the solid angle 2xA(cos@)=2xbA\/2. The 
average over-all orientations, therefore, gives 
rise to the factor A(cos@)/2=bAd/4. The final 
result for the scattering cross section per nucleus 
corresponding to scattering of a particular order 
from a particular set of lattice planes is, apart 
from numerical factors, 


on = NF?/b, (3) 


where N is the number of nuclei per unit volume. 
This cross section is to be summed over all sets 
of crystal planes and all orders of reflection which 
are consistent with the Bragg condition; i.e., for 
which 6<2/)X. 

From this result it can be seen that for very 
low neutron energies of wave-length large com- 
pared to twice the spacing between any adjacent 
lattice planes, no Bragg reflection can occur so 
there is no coherent scattering. When the energy 
is increased to the point at which the Bragg con- 
dition is just satisfied for the most widely spaced 
pair of lattice planes, there is a discontinuous 
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jump in the scattering given by Eq. (3), Then 
the cross section decreases proportionally to the 
square of the neutron wave-length until the 
Bragg condition is satisfied for the next mogt 
widely spaced pair of lattice planes. At that 
point, another term of the type Eq. (3) is intro. 
duced so the cross section again increases discon. 
tinuously. This behavior repeats itself as the 
energy is increased so a jagged curve of the form 
shown in Fig. 1 is obtained. 

The exact form of the elastic scattering crogg 
section with the correct numerical factors jp. 
cluding a factor which takes into account the 
zero point and thermal oscillations of the crystal 
is? 


Oe = > 9<2/,F NX exp(— wh) /8xb, (4) 


where w is a constant which depends on the 
Debye temperature and temperature of the 
lattice and is to be obtained from reference 2, 
Eq. (30). F is to be interpreted as the form factor 
of the unit cell in the crystal and N the number 
of unit cells per unit volume. Then «, is the 
cross section per unit cell. 

Since the Be lattice is hexagonal close packed, 
the form factor is given by 


Fre = 2ene[ 1 +cosr(2/+4m+3n)/3 ]/us, (5) 


where ope is the cross section of the free Be 
atom, wpe is the reduced mass, in units of the 
neutron mass, of the neutron and a free Be atom, 
and /, m, n are the products of the order of the 
reflection and the Miller indices of the plane 
leading to the particular Bragg reflection under 
consideration. The appearance of the factor 
1/ux.? is owing to the very large effective massof 
a Be atom which scatters neutrons with a 
energy small compared to the binding of the Be 
in the lattice.6 The theoretical curve obtained 
from Eq. (4) for Be at room _ temperature 
(7 =293°K) is given in Fig. 2. The cross section 
ope has been taken to be 6.110-*4 cm*. This 
value is obtained from the mean experimental 
cross section for a neutron energy of the ordet 
of several volts (an energy at which the Be atom 
may be treated as free). The Debye temperature 
of Be has been taken to be® © = 1000°K. 


5 E. Fermi, Ricerca Scient. 7, No. 2, 13 (1936). See also 
R. G. Sachs and E. Teller, Phys. Rev. 60, 18 (1941). 

6 Mott and Jones, Properties of Metals and Alloys (Or 
ford, 1936). 
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The generalization of these considerations to 
polyatomic crystals leads again to Eq. (4). In 
this case, however, the form factor depends on 
the positions of the different atomic species in 
the unit cell and on the relative signs of the 
amplitudes of the waves scattered from the dif- 
ferent nuclei. If our considerations are limited to 
biatomic crystals, the form factor is given by 


F=4n|A1 Dose exp[2ri(la+mb +nc) 
+As Das, exp[2ri(lat+mB+ny)]|*, (6) 


where abc are the positions in terms of the 
primitive translations of the lattice of the nuclei 
which scatter with amplitude Ai, and ay are 
the corresponding positions of the nuclei which 
scatter with amplitude Az. For the BeO crystal 


this reduces to 
Fpeo =44|ApetAo exp3rin/4|*F pe 


where Fz, is to be obtained from Eq. (3). In 
particular, for m= +1, the amplitude dependent 
factor is 


Aze+Ao?—V2ApAo. 


Since Age and Ao are of the same order of mag- 
nitude, this quantity depends rather sensitively 
on the relative signs of the amplitudes Ax», and 
Ao. It is larger if the signs are opposite. Thus, 
for the reflections from the twelve sets of planes 
(I, m) =(+1, 0), (0, +1), (1, 1), (—1, 1) with 
n= +1, there will be a considerable difference in 
the strength of the Bragg peak between the two 
cases of like sign and opposite sign. This makes 
it possible to determine experimentally by a 
transmission measurement whether the signs 
Ag, and Ao are equal or opposite. 

The theoretical curves for the cross section 
as a function of energy are given in Fig. 2 for 
both possible choices of sign. In obtaining these 
curves, the magnitudes of the scattering ampli- 
tudes were obtained from 


4x|Ao|?=c0/uo*, (7) 


where go is the scattering cross section of a free 
O atom (found as before from the scattering of 
one volt neutrons to be 4.110-** cm?) and yo 
is the reduced mass of the neutron plus a free O 
atom. The Debye temperature of BeO was 
estimated from the velocity of sound and found 
to be @=1200°K. Comparison of the experi- 


4n|Ane|?=one/uBe’, 
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mental points and the theoretical curves shows 
clearly that the signs of the amplitudes of the ‘Be 
and O scattering are the same. 


5. CONCLUSION 


The good agreement between theory and 
experiment indicates that the application of Eq. 
(4) to other crystalline materials should give an 
adequate representation of the scattering cross 
section as a function of energy. There are, how- 
ever, a number of conditions that may lead to 
deviations from this result. 

Probably the most important of these would 
be the occurrence of such large microcrystals in 
the sample that the peaks in the scattering are 
wiped out by extinction in the crystals. An 
estimate of the maximum size of the micro- 
crystals to be used can be obtained by deter- 
mining what fraction of those neutrons of a given 
énergy which hit a particular crystal at the 
Bragg angle are scattered. This fraction can be 
obtained by dividing fhe expression (1) by the 
projected cross-sectional area of the crystal 
which is given by ? cos@=Pb\/2. If we introduce 
the correct numerical factor which can be ob- 
tained by comparing Eqs. (1)—(4), the fraction 
of the neutrons scattered by a particular set of 
crystal planes turns out to be 


NP F/rb’. (8) 


Since it is desirable to have not more than 1 
percent of the beam scattered by one micro- 
crystal, the upper limit on the crystal dimen- 
sions may be taken to be 


1< (x)'b/10N(F)}. (9) 


If b=3X10' cm, F=10-** cm? and N=4X 10%, 
Eq. (9) becomes /< 5 X 10-5 cm. Thus, for crystals 


‘of linear dimensions smaller than a micron, one 


would not expect an appreciable reduction in the 
peaks because of extinction. These dimensions 
refer to the prefect microcrystals of which the ma- 
terialiscomposed. Considerably larger crystallites 
would be acceptable if the crystallites consist of 
many such microcrystals oriented at small angles 
with respect to each other. In order to be sure 
that this condition was satisfied for the Be 
crystals on which were made the measurements 
discussed in Section 3, it was necessary to break 
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the material up mechanically ; for BeO, there was 
no such difficulty. 

Other sources of deviations from the theory 
may be the following: 


(1) The microcrystals not randomly oriented. This may 
occur when the sample is prepared by extrusion. 

(2) The temperature of the sample not small compared 
to its Debye temperature. Then an appreciable amount of 
inelastic scattering would be expected. This would render 
the peaks less sharp and would introduce appreciable 
scattering below the first Bragg limit. 

(3) The crystal contains two or more isotopic constitu- 
ents which have appreciably different scattering properties. 
The incoherent scattering background would then be 
appreciable. However, it is necessary that there be a very 
strong dependence of the scattering amplitude on the 
isotopic identity of the nucleus for this to be a large effect. 

(4) The nuclei of the sample have a spin different from 
zero and the scattering is strongly spin dependent. This 
would again lead to an incoherent background which 
would be large only if the spin dependence of the scattering 
were very strong. 

(5) If the sample used is made up of a very fine powder, 
the surface area will be so large that there is a chance 
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that adsorbed surface films may contribute appreciably tp 
in 





the cross section. This may lead to an erroneous j 
tation of the results, particularly if one is trying to distin. 
guish between two possibilities Corresponding to 

different combinations of the signs of the sca “= 
amplitudes, If a large amount of water, for a 
adsorbed on the crystals, there will be a smooth ; 
cross section rising rapidly at low energies which could 
easily transform the lower curve in Fig. 2 into a curve 
that could hardly be distinguished from the upper curye, 














If care is taken to avoid the difficulties men- 
tioned above, it is believed that the transmission 
method offers a good procedure for determining 
the relative phases of the scattering of pairs of 
nuclei which can be compounded to form a micro. 
crystalline material. This would then make jt 
possible to determine the absolute sign of the 
scattering phases of a series of nuclei if the phage 
were determined for one particular nucleus, 

The numerical work required for determining 
the theoretical curves in Figs. 1 and 2 was carried 
out by M. G. Goldberger. 



























PHYSICAL REVIEW 


A method is described for determining the half-lives of metastable ions by varying the 
voltage which sweeps the ions from the ionization chamber of a mass spectrometer. The half- 
life of the state corresponding to the transition C,H1o‘>C;H:*+[CHs] has been found to be 
2.0X10-* second and the value found for the transition CyHiot—>C;He++[CH,] was .1.7 
X10-* second. In answer to the objection that the diffuse peaks should not be prominent 
in the 180° instrument if they are to be attributed to spontaneous dissociation during transit, 
it is shown that this result is to be expected for metastable transitions. The relative number 
of ions in the metastable states initially is estimated. 


INTRODUCTION 


HE diffuse peaks appearing in the mass 
spectra of hydrocarbons have recently been 
attributed to the spontaneous dissociation of 
metastable ions during transit.4? Several workers 
in the field have objected to this interpretation 


* Present Address: National Bureau of Standards, 
Washington, D. C. 
a ob. Hipple and E. U. Condon, Phys. Rev. 68, 54 
J. A. Hipple, R. E. Fox, and E. U. Condon, Phys. Rev. 
69, 347 (1946). 
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on the basis that the appearance of diffuse peaks 
in the 180° instrument is unexplained—the peaks 
having roughly the same appearance as in the 
sectored-field instrument such as that employing 
a deflection of 90°. In the latter instrument there 
is a considerable distance between the ion source 
and the deflecting magnetic field where. the 
dissociation could occur without an appreciable 
deflection of the ions, and this situation does not 
obtain in the former case. This paper will 
attempt to show that a diffuse peak observed 
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METASTABLE 


with the 180° instrument will be quite similar in 
ition and contour with the corresponding peak 
with the sectored-field instrument. In order to 
do this, it will be necessary to determine approxi- 
mately the half-life of the metastable ion. The 
contour will then be constructed for sim- 
plified conditions to ascertain whether the actual 
contour can be explained reasonably by dissoci- 
ation during transit. 


DETERMINATION OF HALF-LIFE 


Table I is a reproduction of a portion of Table 
II from an earlier paper.? The more rapid increase 
with ion-draw-out voltage (@=Vi—Vo) of the 
peaks at the apparent mass numbers m* = 39.2, 
31.9, and 30.4 relative to the rest of the spectrum 
of n-butane was a potent argument for their 
origin being the spontaneous dissociations : 


C;H;* (43) —-C;H;* (41) 
+(H2] (2) m*=39.2, (1) 


CyHiot (58)—>C;H;* (43) 
+[CH3] (15) m*=31.9, (2) 


CyHio* (58)—>CsHe* (42) 
+[CH,] (16) m*=30.4. (3) 


As before, brackets are placed around the neutral 
fragment since the mass spectrometer cannot 
determine the nature of the neutral fragment or 
fragments. The metastable ions will decay ex- 
ponentially with time, and the higher values of 8 
will cause them to reach the region where they 
may be detected as diffuse peaks earlier in their 
lifetime leading to a higher peak. A study of 
Table I in conjunction with the geometry of the 
instrument should make possible an estimate of 
the decay constant, X. 

Referring to Fig. 1, let it first be assumed that 
all ions that dissociate later than ¢, will con- 
tribute to the height of the diffuse peak. This 

TaBLE I. Relative peak heights at various voltages 


(Vi— Vo)=8 (from Table II of Hipple, Fox, and Condon, 
Phys. Rev. 69, 347, 1946). 











a’ a’ , , ul 
6 s SS we Bo We 
10 1.00 1.00 1.00 1.00 1.00 
20 173 157 222 221 3.82 
30 1.78 167 262 2.72 5.26 
40 175 169 311 297 5:52 
50 175 169 333 3.29 6.07 








IONS IN 





MASS SPECTRA 


oe 
































E 
1, Va trs23--------- —-- 
% ty», Cc slo oil . 
Ee j ® R Se 
= 8 -{---- 


“ar te, © 















iittala lc sill eee ia al 


Fic. 1. Schematic diagram of ion source to relate the 
time of dissociation in terms of the geometry and the 
accelerating fields. 


should hold fairly accurately for the sectored- 
field instrument if the half-life is quite short 
compared with the transit time through the 
analyzer. In this discussion only the field EZ, will 
be changed so #2(8) is defined as the value of ¢, 
when V;— Vo=8, and n°(8) is the number of ions 
in a particular metastable state passing through 
the slit in electrode 3 when V,;— V,=8. Con- 
sidering the peak at m* = 31.9, the collected ion 
current of mass 43 having 43/58 of the original 
energy will then be given by 


n°(8) = mo exp —([d:(8) ], (4) 


n°(8) may be interpreted either as the number 
of ions of mass 58 left in a metastable state after 
passing through the slit in electrode 3, or the 
number of ions of mass 43 having 43/58ths of the 
original energy received at the ion collector 
(assuming all dissociate before reaching the col- 
lector and all these are collected). mo is the 
number of ions of mass 58 originally in the beam 
in the metastable state responsible for the diffuse 
peak at m*=31.9. As the normal stable ions 
increase in intensity with 8, the metastable ions 
will also have this normal increase a(8) and 
hence 





n°(8)a(8) =n(8). (S) 
Equation (4) now becomes 


n(8) = a(8)mo exp —[Af2(8)]. (6) 
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Fic. 2. Decay curve for the metastable state responsible 
for the diffuse k at m*=31.9 in n-butane. This is a 
plot of the data of Table I and Table II. For A mass 40 
was used in normalizing and for B mass 58. 


Equation (6) was derived with the assumption 
that all metastable ions passing through the slit 
in electrode 3 would contribute to the height of 
the diffuse peak regardless of the point at which 
they dissociate. Now suppose that, because of 
the geometry of the instrument, those ions 
which have travelled a distance greater than L 
in the analyzer will not contribute to the peak. 
Since the ion accelerating field E2 is held con- 
stant, and most of the energy is acquired be- 
tween electrodes (2) and (3), the velocity of the 
metastable ions remains essentially constant. 
Therefore, it can be assumed that all the ions 
contributing to the ion current being measured 
will dissociate between #2(8) and #.(8)+-7 where 
t is a constant. From Eq. (6) this number is 
given by 


An(B) = —a(8)mo exp—([A[t2(8)+7]] 
+a(8)mo exp—[Al2(8)]. (7) 


Normalizing to 8B=Bo 
An(8) a(8)mo exp—[A42(8) J[1—exp—[ar]] 
An(Bo) a(o)moexp—[s(80) [1 —exp—[Ar]] 


An(8) (8) : 
= N(t2(Bo) —te 8 
AntBe) (Be) exp[A(t2(80) —42(8)) ] (8) 


From Eq. (8) it is evident that the half-life 
may be determined without regard to those ions 
which dissociate at a time later than /2(8)+7, 
i.e., considering only ions dissociating between 
te(8) and #.(8)+ 7. This, in turn, means that the 
half-life may be determined by varying 6 even 
though the half-life is not short compared with 








the transit time. The values of 
n' (8) =An(8)/An(Bo) 


and a’ (8) =a(8)/a(8o) for By=1 volt are given in 

Table I. It remains to determine ¢,(8), 
Referring to Fig. 1, it can be shown that », 

V2, ti, te and ts are given by the following relations. 


Y= (2a151)!, (9) 
V2 = (2a/2+2a151)!, (10) 
i= (251/a1) i, (11) 
te = (ve—01) /ag th, (12) 
and 
t3=13/ve+te, (13) 


where a, = E,e/300m and a2= E.¢/300m, and B, 
and Ez are expressed in volts per centimeter, ¢js 
the charge of the ion in e.s.u., and m is the mag 
of the ion in grams. In the apparatus used fo 
obtaining the data of Table I, 14:=0.39 om 
1,=1.0 cm, and s; is estimated to be 0.25 om, 
The ion accelerating voltage was 600 volts and, 
therefore, 


E:=600/1.0 = 600 volts/cm. 


The values of ¢2(8) for various values of 8 ar 
shown in Table II together with the com 
sponding values of n’(8)/a’(8). In column 3 th 
values of a’(8) for mass 40 are used and in 
column 4 the values for mass 58. The values d 
n'(B)/a'(8) are plotted in Fig. 2. It would seem 
more reasonable to use the a’(8) for mass 58if 
the variation in height of a normal peak with 
changes in £ is to be attributed to discriminating 
effects occurring before the ion enters th 
analyzer. The experimental results indicate that 
the change in the height of mass 58 due to th 
variation in the transit time of those ions d 
mass 58 in a metastable state is small. From 


TaBLE II. Normalized peak height (m’(8)/a'(6)) d 
m* = 31.9 for various values of 8 (ion-draw-out vol 











n') ") 
a (8) 26) 
31.9 319 
B t2(8) X 108 sec. ry 38 
1.0 3.86 1.00 1.00 
2.0 2.86 1.41 1.28 
3.0 2.40 1.63 1.53 
4.0 2.21 1.76. 1.70 
5.0 1.94 1.94 1.88 
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METASTABLE IONS IN MASS SPECTRA 


Fig. 2, it is seen that the half-life is approxi- 
mately 1.6X10-* second when correcting with 
mass 58 and 2X10~* second with mass 40. 


REMEASUREMENT OF THE DECAY CURVE 


The results obtained with the data previously 
published suggested the desirability of repeating 
the experiment for a more extended range of £. 
This has been done with another mass-spec- 
trometer tube of the same type. The main ion 
accelerating voltage was raised to 1000 volts and 
8 varied from 0.18 to 10.0 volts. The electron 
accelerating voltage was 75 volts. The data is 
summarized in Table III and points on Fig. 3 
are plotted from this data in the same manner 
as curve B of Fig. 2 with the exception that it 
is assumed that the measured voltage scale must 
be corrected by adding 0.2 volt to the measured 
values. This is the same type of correction em- 
ployed in the measurement of appearance poten- 
tials with the mass spectrometer and could be 
attributed to contact potentials, polarization 
effects, or incomplete shielding of the ionization 
region from the strong fields between neighboring 
electrodes. In the present example this added 
voltage is only 0.02 percent of the main acceler- 
ating voltage and it is reasonable to expect the 
field E, to be shielded incompletely because of 
the slit in electrode 2. This correction is only 
important for small values of the ion-draw-out 
voitage. The experimental points are quite close 
to a stright line on the semi-log plot over a 
range from 1.5 to 6X10~* second. The half-life 
from this curve is 2.0 10~* second. Probably the 
greatest source of error is caused by the difficulty 
in determining s; accurately. 

TABLE III. Data similar to Table II but obtained with 


a different mass-spectrometer tube using 1000-volt ions 
and an extended range of 8. 








n’(8)/a'(B) 
31.9 30.4 
58 58 
185.6 . 522.0 ¥ 1.00 
393.5 ° 1068. ‘ 1.41 
530.0 d 1410, 1.84 
551.0 J 1436, 2.27 
547.0 J 1410, 2.63 
556.0 d 1442. 4.09 
.20 d 560.0 — 1386. 4.39 . 
10.20 1.41 444.0 91.0 1006. 5.1 5.33 4.58 
The values below were obtained several days later and although a’(8) 
and »’(8) are different, their ratio should fall closely on the curve if 
the method of correcting for the “normal” increase with ion-draw-out 
voltage is effective. 
1,70 . J 97.0 1154. 3.0 1.7 2.92 2.34 
2.45 2.53 482.0 98.8 1200. 3.9 2.3 3.75 3.13 


#2(8) X 108 
(sec- 


Peak height (arbitrary units) 
onds) 58 57 29 31.9 0. 
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Fic. 3. Semi-log plot as in Fig. 2 but with new data over 
a more extended range. 


PEAK CONTOUR IN 180° INSTRUMENT 


Having determined the half-life, the next step 
is to relate the time at which the dissociation 
occurs with the apparent mass registered by the 
instrument. The simplifying assumptions will 
first be made that all the ions enter the analyzer 
normal to a line drawn between the entrance and 
exit slits, and the width of the entrance slit will 
be neglected. The metastable ion of mass mo 
then enters the analyzer and first travels in a 
circular path of radius ro; after an angular 
deflection @ or a time TJ after entering the 
analyzer, it dissociates and the ionic component 
of mass m then travels a circular path of radius 
ro* = (m/mo)ro (Fig. 4). When dissociation occurs 
at the moment of entry into the analyzer then 
r* =ro*. The time T is given by : 





T =— cos™! 
V2 


ro = — my) | 


2(2m** — mo*) (mo! — mo*?) 
(14) 
In this equation mo*, mo, v2, and ro are known 
from the experimental conditions and therefore 
T is known as a function of the apparent mass. 
Since the half-life has been determined, the 
intensity as a function of T and through (14) as 
a function of the apparent mass m* is known. 
In determining the peak contour of the diffuse 
peak at 31.9 in n-butane, mo=58 and my* 
= (43)?/58. Thus, the number of ions between 
m* and m*-+Am* will be 
n(m*) —n(m*+-Am*) 
= amy exp — [Ae |[exp — [AT (m*) ] 
—exp[AT(m*+Am*)]}]}. (15) 
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Fic. 4. Dissociation during transit through 180° ana- 
lyzer, The dissociation may occur after appreciable travel 
in the analyzer without a large shift of the apparent mass 
from mo*. 


In order to plot the peak contour the term 
any exp—[Ate] may be taken as unity. If the 
half-life is 2X10-* second then \=0.346X10** 
sec. Taking intervals of Am* the corresponding 
values of the time T are obtained from (14). If 
we assume that the tube is designed for a radius 
of curvature in the analyzer of 12.5 cm (ro* = 12.5 
cm) then r79=16.9 cm. For ions of mass 58 and 
kinetic energy of 600 ev, v2=4.46X 10° cm/sec. 
The values of »(m*) —n(m*+Am*) for successive 
values of m* are given in Table IV and from this 
data the histogram of Fig. 5 was plotted. This 
shows the general shape of the peak contour for 
the diffuse peak at mass 32 in n-butane which 
would be expected when this molecule is studied 
with the 180° instrument and the beam does not 
hit any part of the analyzer during transit. It is 
clear that the maximum appears at essentially 
the same position for the 180° instrument as for 
the sectored-field instrument. For a smaller 
value of the half-life, the ion current would have 
been even more concentrated about m* = 32. 
Now in the 180° instrument it is conventional 
to place fairly wide slits at various points in the 
analyzer to prevent the formation of ‘‘ghosts’”’ by 
reflections from the walls. These slits will cut 
out all those ions in Fig. 5 except those having 
an apparent mass close to 31.9. For instance, 
suppose the excursion from the median ray (ro*) 
is limited to a distance g in Fig. 4 by a slit at an 
angular deflection of 90° for this ray. Then for 
metastable ions dissociating after a small angular 
deflection @ in the analyzer, g will be given by 


&=(ro—10*)(9). (16) 
If g=0.6 cm, then @=8°. In terms of m* this 


means that ions of m*>31.95 (Fig. 5) wil} be 
eliminated. 

It is apparently true that ions should be 
recorded only very close to m* = 31.9 in both the 
180° instrument and the sectored-field instry. 
ment. The intensity should be somewhat 
in the sectored-field instrument. The diff 
however, is not as great as one might expect, Jy 
the example just given, approximately 17 peroens 
of the ions dissociating after entering the 1g 
analyzer would be collected. With the sectored. 
field instrument the ions are rapidly defocusged 
when the dissociation occurs in the magnetic 
field. If a 90° instrument with a radius of curya, 
ture of 12.5 cm is considered (such as that used 
for obtaining the data used in this paper), the 
distance of travel between the ion source and the 
magnetic field region is about 12.5 cm when the 
fringing field is neglected. However, since the air 
gap is ~2.5 cm, the fringing field is large fora 
considerable distance from the magnet and jt 
would be unreasonable to assume that all ions 
resulting from dissociation in the first 12.5 cm 
of path in the analyzer will be properly focussed 
at m* = 31.9. If this distance is assumed to be as 
great as 8 cm, then 7=1.8X10-* sec. and 4 
percent of the ions will be collected, which is only 
about 3 times the value for the 180° instrument, 
The diffuseness of the peaks in each case could 
be attributed to a mutual kinetic energy im 
parted to the fragments at dissociation, or the 
defocussing of the paraxial rays as a result of 
dissociations occurring in the magnetic field. 

However, some of the ions will strike the tube 


TaBLE IV. Number of ions as a function of m*. 











n(m*) 
—n(m* +An*) 


6 
(radians) (microseconds) n(m*) 





0 0 1.00 
0.075 0.284 .906 0.094 

.139 527 833 
182 .690 788 
217 822 752 
246 932 724 
272 1.031 .700 
297 1.125 .678 
340 1.288 .640 
377 1.429 610 
412 1.561 583 
443 1.679 559 
473 1.792 538 
.500 1.895 519 
527 1.997 502 
551 2.088 486 
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APPARENT MASS (M*) 
Fic. 5. Histogram illustrating that a large number of 
the ions dissociating after entering the analyzer of the 180° 


instrument will have an apparent mass close to mo* 
(=31.89). Ions having m* very much different from mpo* 


will be lost in the analyzer. 


walls in passing through the 90° analyzer so the 
percentage getting through would be less than 
46 percent. The value obtained would depend 
somewhat on the position of the main magnet 
relative to the mass-spectrometer tube. 

In discussing the peak contour the width of 
the entrance slit was neglected and the distri- 
bution in intensity was obtained along the focal 
plane of a normal ion beam (not metastable). 
Actually the data was taken with an exit slit at 
a fixed position and the spectrum was scanned 
by varying the magnetic field. This makes no 
essential difference in regard to the general 
conclusions. Clearly, of course, the exit slit must 
be sufficiently wide to allow the passage of ions 
in the range between m* = 31.89 and m* =31.95 
(see Fig. 5) in order that the peak height will be 
a measure of the number of metastable ions 
that have passed through the analyzer without 
being lost. This requirement is fulfilled in the 
analytical instruments of the 180° type in 
general use today. The entrance slit to the 
analyzer is generally much smaller than the exit 
slit and for this reason will not materially affect 
our conclusions. 


RELATIVE NUMBER OF IONS IN THE METASTABLE 
STATE INITIALLY 


Since the half-life has been determined, it 
should now be possible to estimate the relative 
number of ions of mass 58 originally in the meta- 
stable state responsible for the diffuse peak at 
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Fic. 6. Decay curve for metastable state responsible for the 
diffuse peak at mv* = 30.4 in n-butane. 


m* = 31.9. For 8=3.2 volts, the peak height of 
m* = 31.9 is 4.9 divisions from Table III. If it is 
assumed on the basis of the previous discussion 
that only 40 percent of the ions that dissociate 
in the analyzer contribute to this peak height, 
then the peak height must be multiplied by 2.5 
to obtain the value of °(8), i.e., 2°(8) = 12.2. All 
the quantities in (4) are now known except mo 
and this may be calculated: 


mo = 12.2 exp[ (0.346 X 10°) (2.25 K 10-*) ] = 26.7. 


This is only about 5 percent of the number of 
ions of mass 58 which reach the ion collector 
without dissociating. 

The decay with time of the ions in the meta- 
stable state responsible for the diffuse peak at 
m* =30.4 is shown in Fig. 6. The transition 
involved here is 

CyHipt->CsHo+ + [CH]. 


The data were obtained from the same charts 
used for the points in Fig. 3. This peak was 
smaller than that at m*=31.9, and the current 
for 8=0.38 volts was too small to include in the 
measurements. From this curve the half-life is 
1.7 X 10~ second and the corresponding value for 
mo is 18.8 divisions estimated in the same way 
as described for the other metastable transition. 
This is between 3 and 4 percent of the number 
of ions of mass 58 which reach the ion collector 
without dissociating. 
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With a resonance method, the velocity of the temperature waves (second sound) in liquid 
helium II has been determined as a function of temperature from about 1.4°K to the A-point. 
The velocity is found to be zero at the \-point, reaches a maximum of 20.46 m/sec. at 1.7°K 
and decreases to 19.80 m/sec. at 1.42°K. The accuracy of the measurements is at least +0.5 
percent. The results are in good agreement with the. theoretical predictions of both Tisza 
and Landau but tend to favor the former’s theory over the latter’s. The fact that second sound 
generated in the superfluid is able to register on a microphone diaphragm immersed in the 


fluid is discussed. 





INTRODUCTION 


N an attempt to explain the many anomalous 

properties of liquid helium at temperatures 
below 2.19°K (liquid helium II), F. London has 
proposed a model for this fluid which is based 
on an ideal gas obeying Bose-Einstein statistics. 
In these statistics a peculiar “‘condensation”’ 
effect takes place at a certain critical temperature 
T. whereby a certain fraction of the atoms go 
into a highly degenerate ground state, the 
remainder remaining in excited states. In par- 
ticular the number of-degenerate atoms Np is 
given by 


men{i-(Z)'], rere 


N being the total atoms present per unit volume 


and 
h n i 
T.= ( ) ; (2) 
2armk \ 2.6124 


By putting the appropriate values for helium 
into Eq. (2) London finds T,=3.14°K which is 
surprisingly close to the \-point (2.19°K) con- 
sidering that we are dealing with an ideal gas 
and not a liquid. Below 7, the entropy of the 
system is caused solely by excited atoms, the 
degenerate ones have zero entropy. 

Using this admittedly rough model, London 
has been able to show that a number of known 
properties of liquid helium II are theoretically 
derivable, at least qualitatively. L. Tisza has 
improved the London model in a number of 
particulars. He supposes liquid helium II to be a 
“‘mixture’”’ of two kinds of atoms, superfluid of 





density p, and normal of density pn». These cor. 
respond to the degenerate and excited atoms jp 
London’s Bose-Einstein gas and again the 
entropy of p, is zero. A further assumption made 
by Tisza, which appears somewhat difficult to 
justify, is that the interaction between p, and », 
is negligible and indeed p, atoms do not exchange 
momentum with either p, or other atoms. Tisza 
further supposes that the density function for 
degenerate atoms is 


Ps=p ae ’ <lLe 
T. 


P= Pst Pn, T.=2.19°K, 


and o&5 to 6 as against $ in the ideal Bose gas, 
The justification for this latter modification lies 
partly in the better fit obtained with the exper- 
mental results for the specific heat discontinuity 
at the A-point and partly because of the fact that, 
as London had observed, it is probable that the 
introduction of interaction between the atoms 
would modify the density function in this 
direction. 

From this modified theory and a consideration 
of the hydrodynamical equations of the two 
phases, Tisza in 1940' predicted that temperature 
differences in liquid helium II should be prope 
gated as a wave motion with a characteristic, 
temperature dependent velocity. 

The experimental proof of this prediction was 
made in 1944 by V. Peshkov,? and these tem 
perature waves were named “second sound’ 


1L. Tisza, J. de phys. et rad. [8] 1, 164, 350 (1940). 


2V. Peshkov, J. phys. et rad. U.S.S.R. 8, 381 (1944). 
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Peshkov generated these waves by means of a 
wire resistance coil in which a variable tem- 
perature was maintained by means of an alter- 
nating current, and detected them from the 
indications of a small phosphor bronze resistance 
thermometer. By measuring the wave-length of 
running or standing waves (produced by reflec- 
tion from the liquid surface) and knowing the 
heat waves to be double the frequency of the 
exciting a.c. he could determine the velocity 
which turned out to be of the order of 20 m/sec., 
i.e., around 10 times less than the velocity of 
ordinary sonic waves. 

In Tisza’s theory the velocity of second sound 


is given by 
pn 9 f/1\T? 
30) 
p,"dT S 


where S is the entropy of helium II, which has 
been accurately determined as a function of tem- 
perature from measurements on the specific heat. 
Accordingly measurements of C as a function of 
temperature would yield information on the 
temperature variation of the important ratio 
p:/pa. The purpose’ of the following experiment 
is therefore clear. In addition to this we must 
mention that the existence of second sound was 
predicted independently by L. Landau,’ some- 
what later than Tisza, by use of an entirely dif- 
ferent kind of theory. We shall discuss the 
bearing of the experimental results on the two 
theories at the end of this paper. 


EXPERIMENTAL DETAILS 


Following discussions with L. Onsager, of this 
laboratory, we decided to attempt to measure 
second sound velocities with the following kind 
of arrangement. The second sound would be 
generated by an a.c. heater of the same general 
type as that used by Peshkov, and standing 
waves would be produced in a column of liquid 
helium by reflections from the free surface. If 
this free surface were an anti-node, the temper- 
ature of this surface would fluctuate at double 
the frequency of the a.c. source, giving rise to a 
fluctuating evaporation of the liquid at this 
surface. In consequence small periodic fluctua- 





*L. Landau, J. Phys. U.S.S.R. 5, 71 (1941). 
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tions in the pressure of the vapor above the 
liquid surface would occur (first sound). These 
fluctuations of pressure would register on a 
microphone immersed in the helium vapor, and 
hence second sound resonances in the liquid 
would be detectable. 

Figure 1 is a schematic diagram of our ap- 
paratus. He and Ni are, respectively, the liquid 
helium and liquid nitrogen Dewars. The latter 
is of Pyrex with a highly evacuated interspace. 
The helium Dewar, on the other hand, is of soft 
glass and the interspace is filled with air at a 
pressure of about 2 millimeters mercury. This 
has the following purpose. Before liquid helium 
is admitted, the apparatus in this Dewar is 
cooled to the temperature of liquid nitrogen, the 
air in its interspace acting as the heat transfer 
medium. A large fraction of the heat content of 
this apparatus is thereby removed and the load 
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Fic. 1. Schematic{drawing of the apparatus. 









LANE, FAIRBANK, AND FAIRBANK 








™ RESONANCES IN LIQUID 
1-57 °K 

















Fic. 2. Second sound resonances as registered on recording potentiometer at 
1.57°K. The microphone is in the vapor above the liquid surface. 


on the liquid helium accordingly much reduced. 
This is important in view of the very small latent 
heat of liquid helium. When the liquid helium 
is now admitted to the Dewar, the air in the 
interspace freezes and a very high vacuum is 
produced, thereby providing gdod thermal insula- 
tion for the liquid.‘ Both Dewars were fully 
silvered except for a pair of narrow longitudinal 
slits in each permitting visual observation of the 
the interior of the helium Dewar. The latter was 
about 18 inches long and 1 inch internal diameter. 
L is a thin walled Lucite tube approximately 15 
cm long and 2 cm in diameter. It is provided 
with four small holes (G) to permit filling from 
the Dewar. H is a constantan heater, wound 
with 3-mil wire, a detail of which is shown 
separately. This heater is fed by number 40 
copper wire (W) twisted to eliminate magnetic 
coupling into the receiver circuit. Several Lucite 
spacers (F) are used to center the cavity in the 
Dewar. Fitted tightly into the Lucite cavity is 
the microphone (M). The specifications for this 
device were severe. It had to be small, light in 
weight, and capable of operating at temperatures 
near absolute zero. The operating temperature 


4 This technique originated at the Royal Society Mond 
Laboratory, Cambridge, and was communicated to us by 
Professor Onsager who learned of it on a recent visit to 
that laboratory. We found it to be immensely useful in all 
sorts of work with liquid helium and also liquid hydrogen. 


ruled out the usual Rochelle salt-crystal type 
and we finally made use of the receiver of a 
Sonotone hearing aid used in “‘reverse’’ which is 
a magnetic type® instrument. The whole as- 
sembly is supported from the top of the Dewar 
by a small, thin-walled, stainless steel tube (5), 
This tube contains an insulated inner conductor 
and the pair form the shielded microphone leads, 

The heater is fed continuously with 1000-cycle 
a.c., of good wave form, from a stable tuning 
fork oscillator. The power fed into the heater, 
which had a resistance of about 60 ohms at 
liquid helium temperature, was roughly 1/25 
watt. The microphone output was amplified by 
an amplifier roughly tuned to 2000 cycles. This 
latter was provided with both a.c. and rectified 
output. The rectified output was fed intoa Brown 
high speed recording potentiometer with a chart 
speed of 2 inches per minute. The Lucite cavity 
was illuminated through the slits in the Dewars 
by means of an external fluorescent lamp and the 
liquid helium level in the cavity was observed 
with a cathetometer. The vapor pressure in the 
cryostat was read on an absolute manometer 
filled with diffusion pump oil and was controlled 
manually by means of a needle valve. The 1937 


5 We are indebted to Dr. L. Grant Hector of Sonotone 
Corporation for his gift of several of these devices. 
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Leiden Scale was used in computing the cor- 
responding absolute temperatures. 

The procedure in taking data was as follows. 
As the liquid helium slowly evaporated and the 
level in the cavity fell, cathetometer readings 
were taken as a function of time. This clock was 
synchronized with the recorder chart, which 
recorded a series of resonance peaks. By plotting 
the liquid level as a function of time (this was 
always a straight line when the temperature was 
steady) the height of the liquid column in the 
cavity at each resonance was known. The dif- 
ference in height between two successive reso- 
nance peaks being, of course, one-half wave- 
length. Usually about 6 to 8 peaks were taken 
at each temperature. Figure 2 shows a typical 
series of such resonance peaks at 1.57°K. One 
pair of plates of an oscilloscope were connected 
to the heater oscillator, the other pair to the a.c. 
output of the microphone amplifier. The resulting 
Lissajous figure confirmed that the sound fre- 
quency was double the electrical heater fre- 


quencyse 
RESULTS 


Using the previously described method we 
determined the second sound velocity at a 
number of temperatures from about 1.4°K to 
the A-point. Since we believe our accuracy is 
better than +0.5 percent we present the data in 
both tabular (Table I) and graphical (Fig. 3) 
form. It will be noticed that the velocity extra- 
polates to zero, very accurately, at the \-point 
as the theory requires. 

It is apparent from our data that the liquid 
column, at resonance, is probably oscillating 
with a temperature anti-node at both the free 
surface and the heater surface. This fact enables 


TABLE I. Second sound velocity versus temperature. 








= meters/sec. 

19.80 
19,86 
20.23 
20.38 
20.46 
20.16 
20.08 
17.12 
14.13 

9.40 





1.453 
1.570 
1.607 
1.780 


1.994 
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Fic. 3. Second sound velocity as a function of tem- 
perature. The heavy line with circles is the experimental 
curve. The lighter broken line is the theoretical curve com- 
puted from Tisza’s theory of liquid helium II. 


us to compute the wave-length with great 
accuracy. Thus, at a given temperature, having 
roughly determined the half-wave-length from 
the distance between successive resonance peaks 
it is possible to assign a definite number of half- 
wave-lengths (m) to the whole liquid column at 
each resonance. If then successive resonant 
lengths of liquid column be designated by Lo, 
Ly, «+ +L, we found the following very accurately 
true: 


and therefore 


ae 
(s+1)(n—s/2) 


In most of our velocity determinations S was of 
the order of 6 or 7 and m was usually about 30 
or so and hence a very accurate determination 
of \ was possible. 

In one run we held the temperature constant 
for about 23 successive resonances. The reso- 
nance peaks varied widely in height and width 
and appeared to be modulated by a much longer 
wave-length. We were able to determine this 
modulating wave-length roughly and the result- 
ing velocity was of the order of the velocity of 
normal sound in helium vapor at the prevailing 
pressure and temperature in the vapor space of 
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Fic. 4. Second sound resonances as registered on recording potentiometer with the microphone 
completely immersed in the liquid. 


the Lucite cavity. All our runs showed this 
modulating feature to a greater or lesser extent. 

Much more rarely we observed a double 
resonance peak. Such a one is shown on the 
third peak from the left in Fig. 2. We only ob- 
served two such twin peaks in all our work. In 
both cases only one of the twins fitted our 
second sound data. Assuming the spurious peak 
was resonance of ordinary sound in the liquid 
we could show that it was in accord with the 
known value of this velocity. 

A rather curious effect was invariably observed 
when the liquid surface in the cavity was about 
even with the top of the heater wires. We always 
obtained a pair of second sound resonance peaks 
close together at this stage. 

One run was made which was a departure 
from our standard technique. In this case the 
microphone was completely immersed in the 
liquid helium II and there was no observable 
vapor space in the Lucite cavity as observed by 
the cathetometer telescope.* The temperature 
was set at about 1.75°K and allowed to rise 
slowly and steadily up to the A-point. Figure 4 
shows the recorder chart record in this case. 
Examination of these data, with the curve of 
Fig. 3, leaves little doubt but that we are ob- 
serving true second sound resonances in this 

*In view of the enormous heat oatety of liquid 
helium II no “vapor lock” would be expected. Vapor locks 
in liquid helium I and liquid hydrogen always occurred in 
similar circumstances. 


case. Here, of course, both the frequency and the 
length of the resonant cavity are constant anda 
peak will be obtained whenever the temperature 
is such that the second sound velocity has.a valye 
capable of rendering this fixed system resonant, 
If L is the length of the cavity and fo the steady 
heater current frequency, then we shall observe 
a resonance peak when the second sound velocity 
V(T) is given by 
4 
=, 


n 


where m is a whole number. Accordingly, starting 
at some known temperature (and velocity), near 
the maximum in the curve, and proceeding to 
the A-point we should get a large number (theo 
retically an infinite number) of resonance peaks 
corresponding to velocities given by m, m+, 
n+2, --+- in the above formula. As the tem- 
perature approaches the \-point the peaks will 
be closer together due to the steepness of the 
velocity-temperature curve in this region until 
the speed of the recorder is no longer sufficient 
to record all the peaks. This is approximately 
what we observe. 
DISCUSSION 

As mentioned previously, a theoretical relation 

for the second sound velocity versus temperatutt 


curve has been deduced by Tisza depending oa 
(1) the experimentally determined relationship 
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between the entropy and the temperature, (2) 
p./ pn a8 a function of temperature. If one takes, 
following London,® as an average in the range 
from 1°K to the d-point, 


S=2.1X10575-* erg/°K g, 


7 5.6 
.= oA - ’ 
(—) 


then we obtain 


c=2.6X 104 —I1 —_ (—) |] cm/sec. 


This function is plotted as a broken line in Fig. 3 
and the agreement is seen to be very good. As 
mentioned in the introduction, the Landau 
theory predicts that the temperature variation 
of velocity shall follow a somewhat different 
course than that predicted by the Tisza theory. 
Thus at absolute zero the velocity according to 
Tisza should be zero, according to Landau it 
should be finite. The two theories, therefore, 
diverge in the neighborhood of 1°K, Landau 
showing a rising velocity, Tisza a falling one as 
the temperature decreases. Unfortunately we are 
unable to achieve a temperature below 1.4°K 
with our set-up but in view of the accuracy of 
our measurements and the fact of a distinct 
maximum in our curve, Tisza’s theory seems to 
be somewhat favored over Landau’s. 


* We are uaty indebted to Professor Tisza for drawing 


our attention to his articles in Journal de Physique, copies 
of which had not reached us, and for making available to 
us an advance copy of F. London’s paper before the 
Cambridge conference in 1946. The above formulas are 
drawn from this latter source. Mr. Tisza also communicated 
to us Peshkov’s results, given at the same conference, 
which are in good agreement with ours. 
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There remains the reason for the observance of 
second sound resonances in the case where the 
microphone is totally immersed in liquid helium 
II. There is no reason to expect, from the theory, 
that the temperature waves could give rise to 
pressure fluctuations in the liquid and therefore 
register an effect upon the microphone dia- 
phragm. Professor Onsager has suggested to us 
that the mechanism may lie in the nature of the 
cavity which we are using. Referring back to 
Fig. 1 the four small filling holes G will be noted. 
The temperature fluctuations of the second 
sound inside the cavity will, according to Tisza’s 
theory, produce a flow of p, atoms into the 
cavity through the holes and a balancing flow 
of p, atoms out through the holes. The flow of 
p, atoms will be frictionless, but the p, atoms, 
possessing viscosity, will probably experience a 
pressure drop in their passage through the holes. 
Since the flow of the p, and p, will be fluctuating 
we could by this mechanism have fluctuations of 
pressure in the cavity, induced by the second 
sound, which would be capable of registering on 
our microphone diaphragm. 

It will be observed from Fig. 2 that many of 
the resonance peaks are extremely sharp, indi- 
cating a high “Q”’. This suggests that both the 
transmission coefficient across the liquid surface, 
and the dissipation in the liquid, are small. The 
absence of end effects in the resonant column 
also confirms the low transmission coefficient. 

Finally we should like to express our thanks 
to Professor Lars Onsager for his many valuable 
discussions throughout the course of this work. 
We are also indebted to Mr. Robert T. Webber 
for his aid in taking data. This work was finan- 
cially supported by the U. S. Navy under con- 
tract N6ori-44, Task Order III. 
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Theory of the Propagation of Shock Waves* 


Stuart R. BrinKLeEy, Jr.** AND JOHN G. KIRKwoop 
Department of Chemistry, Cornell University, Ithaca, New York : 


(Received January 22, 1947) 


A new theory of propagation of one-dimensional shock waves is described. The partial P 
differential equations of hydrodynamics and the Hugoniot relation between pressure and 
particle velocity are used to provide three relations between the four partial derivatives of 
pressure and particle velocity, with respect to time and distance from the source, at the shock 
front. An approximate fourth relation is set up by imposing a similarity restraint on the shape 
of the energy-time curve of the shock wave and by utilizing the second law of thermodynamics 
to determine, at an arbitrary distance, the distribution of the initial energy input between ‘ 











dissipated energy residual in the fluid already traversed by the shock wave and energy available Ls 
- for further propagation. The four relations are used to formulate a pair of ordinary differential 
equations for peak pressure and shock wave energy as functions of distance from the source. 
The theory takes proper account of the finite entropy increment of the fluid produced by the 
passage of the shock and permits the use of the exact Hugoniot curve of the fluid in the nu- 
merical integration of the basic equations. 
INTRODUCTION based upon the numerical integration of the 
N discussing earlier work by Poisson,! Stokes? Partial differential equations of hydrodynamics, 
appears to have first suggested the possibility However, the labor involved in this approach is 
of the propagation at greater than acoustic %° teat as to limit it to special applications, - 
velocity of discontinuous pressure waves. Earn- The need for a more flexible and rapid theoretical in 
shaw? and Riemann‘ have discussed the laws of ™eéthod, based if necessary on well-defined fl 
propagation of waves of finite amplitude and the approximations, is therefore ey Such a method “ 
building up of the discontinuity. The conditions as been previously developed for shock ware as 
for a wave of permanent type have been investi- i water from spherical charges of explosive er 
gated by Rankine’ and Hugoniot* who have Underwater explosion waves are simpler to treat pe 
provided statements of the conditions for the than blast waves in air, since the relatively small on 
continuity of mass, momentum, and energy ¢®tropy increment produced at the shock front 
across the moving discontinuity. Rayleigh’ has permits the use of the approximation of adiabatic 
shown how to solve the hydrodynamic equations ow. Failure of this he geen = a has 
for plane shock waves when the pressure and Util now prevented the formulation of an 
density are connected by the adiabatic law. adequate theory of blast peng ee g has Ei 
An exact solution for the cases of spherical A theory of eee’ hi ” aaa teal of 
and cylindrical shock waves cannot be given, due Tecently been presented by Osborne and Taylor. in 
to the spherical and cylindrical divergence terms heir theory is based upon the acoustic approx 
of the equation of motion. A straight forward ™ation and is therefore strictly valid only for no 
attack on the mathematical problem may be small excess pressures at large distances from for 
the source. an 
* This paper is a report on a portion of work done at In the present communication, we describe a Eu 
ana ee Rat ees yooh meg et 4 theory of the propagation of one-dimensional to| 
** Present address: Central Experiment Station, U.S. —that is, plane, cylindrical, and spherical— on 
Hu 


Bureau of Mines, Pittsburgh, Pennsylvania. 
1M. Poisson, J. de I’école polyt. 7, 319 (1808). 
2G. G. Stokes, Phil. Mag. 33, 349 (1848). 
*S. Earnshaw, Phil. Trans. A150, 133 (1850). 
4B. Riemann, Nachr. Ges. Wiss. Géttingen 8, 43 (1860). 
5 W. J. M. Rankine, Phil. Trans. A160, 277 (1870). 
a 388} ugoniot, J. de l’école polyt. 57, 3 (1887); 58, 1 
7Lord Rayleigh, Proc. Roy. Soc. A84, 247 (1910). This 
paper gives a full discussion of earlier work. 


shock waves which is valid both in air and water. 

The theory takes account of the finite entropy 

increment in the fluid resulting from the passage 

of the shock wave. It also permits the use of the 
8 J. G. Kirkwood and H. A. Bethe (1941). 


*F. M. Osborne and A. H. Taylor, Phys. Rev. 70, 322 
(1946). 
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exact Hugoniot curve for the fluid and in this 
respect is superior to previous treatments of the 
shock wave in air which have been based upon 
ideal gas adiabatics with constant heat capacity,” 
an approximation which fails badly near the 


explosive source. 


THE PROPAGATION EQUATIONS 


In discussing the propagation of one-dimen- 
sional shock waves, it is convenient to write the 
equations of hydrodynamics in the form, 


where u is the particle velocity, p the pressure 
in excess of the pressure po of the undisturbed 
fluid, p the density, po the density of the undis- 
turbed fluid, ¢ the time, and r the Euler coordi- 
nate at time ¢ of an element of fluid with La- 
grange coordinate R. The Euler sound velocity 
c is equal to [(0p/dp) s |. The coefficient a takes 
on the values, 


a=0 for a plane wave, 
a=1 for a cylindrical wave, 
a=2 for a spherical wave. 


Equations (1) are supplemented by the equation 
of state of the fluid and the entropy transport 
equation 0S/dt=0, the latter of which we shall 
not explicitly use. Equations (1) are of a hybrid 
form in that we use the Lagrange coordinates R 
and ¢ as independent variables but retain the 
Euler equation of continuity. Equations (1) are 
to be solved subject to initial conditions specified 
on a curve in the R, ¢-plane and to the Rankine®- 
Hugoniot® conditions at the shock front, 


D = pou U, 
p(U—u) =poU, 
AH =(p/2)(1/po+1/p), 


”G, I. Taylor, Proc. Roy. Soc. A186, 273 (1946). 
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where AH is the specific enthalpy increment 
experienced by the fluid in traversing the shock 
front and U is the velocity of the shock front. 
Equations (2) constitute supernumerary bound- 
ary conditions which are compatible with the 
differential equations and specified initial condi- 
tions only if the shock front follows an implicitly 
prescribed curve Ro(t) in the R, t-plane. 

We denote a derivative in which the shock 
front is stationary by 


d y:. Big 


> <M, (3) 
dR @R Udt 

If the operator d/dR is applied to the first of 
Eqs. (2) and if Eqs. (1) are specialized for the 
shock front, r=R, three relations are obtained 
between the four partial derivatives dp/dt, 
Op/AR, du/dt, du/AR. 


p ou 


1 dp au 
Po OR pc? at 


R 


Ou 


1 ap 
——=0, 


4 
po OR «) 


Ou 


at 


Ou 


All of the coefficients in Eqs. (4) can be expressed 
as functions of pressure alone by means of the 
Hugoniot conditions, Eqs. (2), and the equation 
of state of the fluid. If Eqs. (4) can be supple- 
mented by a fourth relation between the partial 
derivatives, it would be possible to solve for 
each of the four derivatives as a function of p 
and R, and with the aid of Eqs. (3) to formulate 
an ordinary differential equation, 


‘uth F(p, R) 
aR OR UH.’ 


dp ap 


(5) 


for the peak pressure p of the shock wave as a 
function of the distance R, and in-addition to 
obtain the initial slope, —1/0,"of_the Euler 





608 Ss. 


pressure-time curve of the wave, 
1 10p ap dp 


—-=-——- — — (6) 


6 pdt ppoaR 


The physical basis for the supplementary 
relation to be established lies in the fact that the 
non-acoustical decay of waves of finite amplitude 
is closely associated with the entropy increment 
experienced by the fluid in passing through the 
shock front and the accompanying dissipation 
of energy. As a shock wave passes through a 
fluid, it leaves in its path a residual internal 
energy increment in each element of fluid deter- 
mined by the entropy increment produced in it 
by the passage of the shock front. As a result, 
the energy propagated ahead by the shock wave 
decreases with the distance it has traveled from 
the source. 

The adiabatic work wo per unit area of initial 
generating surface done on the fluid exterior to 
the surface is given by 


Wolo? = f pot o*E[ (ro) jdro 


ao 
« 


a reu'(p’+po)dt, (7) 


to(R) 


where u’ and p’ denote particle velocity and 
excess pressure behind the shock front (the 
unprimed quantities being reserved henceforth 
for quantities at the shock front), to(R) is the 
time of arrival of the shock front at R, E(p) is 
the specific energy increment of the fluid at 
pressure po and for an entropy increment corre- 
sponding to shock front pressure p, and dp is the 
Lagrange coordinate of the generating surface 
of the shock wave. Now, 


f por*u'dt 
teo(R) 
B / po 
= poA V+pof (=- 1 rors (8) 
ao Pp 


where p is the final density of the fluid and AV 
is the volume swept out by the generating surface 
per unit area of initial generating surface. Com- 
bining Eqs. (7) and (8) and introducing A(p) 
= E+ A(1/p), the specific enthalpy increment 
of an element of fluid traversed by a shock wave 
of peak pressure p after return to pressure Po 
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along its new adiabatic, one obtains 


R 
Wolo* = PyAV+ f por o*h[ (ro) |dro 
> rou’ p'dt. 
to(R) ) 
The time integral may be assumed to vanish at 
R=. If one subtracts from Eq. (9) the expres. 
sion obtained from that equation at R= ©, there 
results the relation, 
ru’ p'dt, 
to(R) 


D(R) = (10) 


where 


D(R)= f perohLP(re) dre (At) 
R 


The energy of the shock wave at the point R js 
by definition the work done on the fluid exterior 
to R. Thus, the shock wave energy at R per unit 
area of initial generating surface is given by the 
quantity ao~*D(R). 

Our dissipation assumption breaks down if the 
first shock wave can be overtaken by second 
shocks built up in its rear. This will not be the 
case if the pressure-time curve is initially mono- 
tone, decreasing with asymptotic value pp. If 
the excess pressure p’ has a negative phase, a 
second shock will develop in the negative part 
of the pressure-time curve but cannot overtake 
the initial positive shock. In this case, our theory 
will apply to the positive phase if the time 
integrals of Eqs. (7) to (10) are extended not to 
infinity but to the time at which the excess 
pressure in the positive phase vanishes. The 
general theory of shock waves is not sufficiently 
developed to permit one to say that there is 
proof for these statements, but they can never- 
theless be accepted with some assurance as 
plausible. 

The energy-time integral can be expressed in 


reduced form, 
D(R) = R*pupy, 


1 (—2=") 1d0p 104 au 
B ot t=to(R) pot u ot Kay 

bed t—to(R) 
=f f(R, 1dr, r= ’ 
0 


yp 
S(R, 1) =rep'u'/ Repu. 





dt. (9) 


nish at 
expres. 
>, there 


(10) 


(11) 


it R is 
xterior 
er unit 
by the 
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mono- 
po. If 
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> part 
artake 
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time 
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ever- 
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1s 


R 
(12) 














The function f(R, r) is the energy-time integrand, 
normalized by its peak value R«pu at the shock 
front, expressed as a function of R and a reduced 
time r which normalizes its initial slope to —1 
if ~ does not vanish. We also assume f to be a 
monotone decreasing function of 7. Elimination 
of » between the first two of Eqs. (12) yields the 


desired fourth relation supplementing Eqs. (4) 


between the partial derivatives at the shock 
front. This set is exact, involving integrals of 
Eqs. (1) for the knowledge of the reduced energy- 
time function f(R, 7). However, if f(ao, r) is 
initially a monotone decreasing function of r, 
f(R, 7) will remain so, and in fact will at large R 
become asymptotically a quadratic function of 
r corresponding to the linear form of the pressure- 
time curve that has been shown to be stable at 
large distances. This means that » is a very 
slowly varying function of R, for which suff- 
ciently accurate estimates for many purposes 
can be made without explicit integration of the 
hydrodynamic equations. The assignment of a 
value independent of R to »v is equivalent to 
imposing a similarity restraint on the energy- 
time curve. This type of approximation is equiva- 
lent in principle to that underlying the Rayleigh- 
Ritz method, although we do not include a 
variational procedure to carry the result to any 
desired degree of approximation. 

The initial pressure-time and energy-time 
curves of an explosion wave are rapidly decreas- 
ing. An expansion of the logarithm of the func- 
tion in a Taylor series in the time, the well-known 
peak approximation, is appropriate for an initial 
estimate of v. This corresponds to an exponential 


f(r), 
f(r) =e, 


and results in the value, y= 1. For the asymptotic 
quadratic energy-time curve, corresponding to a 
linear pressure-time curve of the positive phase 
of the wave, 


f(r) =(1—1/2)?, 
f(r) =0, 


which leads to the value, y=3. As a convenient 
empirical interpolation formula between the two 
extreme values, we have employed the relation, 


v=1—} exp—[p’po] (13) 


rQ2, 
t>2, 
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in a series of calculations" of the peak pressure 
distance curves of shock waves from explosive 
sources. 

Elimination of u between the first two of Eqs. 
(12) and combination with Eq. (4) yields a set of 
four equations which may be solved for the four 
partial derivatives, and an ordinary differential 
equation for the peak pressure p as a function of 
the distance R may be formulated with the aid 
of Eq. (5). A second ordinary differential equa- 
tion relating D to R may be obtained by differ- 
entiation of Eq. (11). The resulting expressions 
may be written in the form, 








dD RL 
a. 7 p), 

(14) 
as kaka ye ) 
n° CDS OR 

where 
L(p) = poh(p), 
1 
M(p)= , 
@) poU? 2(1+2)—G 
(py wee TAU = eel OVG 
2(1+g)—G 
G=1—(poU/pc)? = 1 = 
= Po pce)", g= U dp 


The functions L(p), M(p), and N(p) can be 
evaluated as functions of the pressure by means 
of the equation of state of the fluid and the 
Hugoniot relations, Eqs. (2). It may be remarked 
that Eqs. (14) are independent of any assumption 
regarding the equation of state of the fluid and 
that they take proper account of the finite 
entropy increment of the fluid produced by the 
passage of the shock. 

The functions L(p), M(p), and N(p) are most 
conveniently expressed as functions of the pres- 
sure in tabular form, and they may be evaluated 
by numerical methods from a tabular presenta- 
tion of the exact Hugoniot curves for the fluid. 
Equations (14) can then be integrated numeri- 
cally by the use of standard methods.” 


u J. G. Kirkwood and S, R. eee 6 r. (1945). 


‘umerical Mathe- 


12 See, for example, J. B. Scarboroug 
(fohns Hopkins Press, Baltimore, 1930), 


matical Analysis 
pp. 218 ff. 
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THE ASYMPTOTIC PROPAGATION EQUATIONS 


It is of interest to examine the asymptotic 
form of the solutions of Eqs. (14) for small 
excess pressure p. If the exterior medium is 
air, the ideal adiabatic equation of state, 
P= pol. (p/po)” —1], may be employed in the limit 
of small excess pressure, where y is the ratio of 
the heat capacities of p) and at the temperature 
of the undisturbed fluid. Employing the acoustic 
approximation in the evaluation of the limiting 
shock wave velocity, it is easy to show that 


y+1 p* 
12? po?’ 


‘ y+1 p 
Lim(p—0) M(p) = —, 
8y? Po? 
Lim(p—0) N(p) =1, 


Lim(p—0) v= §, 


Lim(p—0)L(p) = 


and the asymptotic equations are 


y+1 Rp? 


dp y+1 Rep 


dR 2R 12? Dpy? 


Equations (16) have the solutions, 


Rp =P,(logR/R:)“, 
D=[(y+1)/6y*po? ]P?Rp, 


for the spherical wave, a=2, 


V/Rp=Pi[2(/R-V Ri) J? 
D=[(y+1)/67°pe JP? Rp, 


for the cylindrical wave, a=1, and 


p = P,(R-—R;)“}, 
D=[(y+1)/6y°p0? ]P xp, 


for the plane wave, a=0, where P; and R, are 
constants. Equations (15) to (17) are valid when 
the exterior medium is water if, in these expres- 
sions, the pressure fo is replaced by the char- 
acteristic pressure B of the Tait! equation of 


(17b) 


(17c) 


% R. E. Gibson, J. Am. Chem. Soc. 56, 4 (1934); 57, 284 


{ga} See also A. Wohl, Zeits. f. physik. Chemie 99, 234 
(1921), and H. Carl, ibid. 101, 238 (1922). 
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state, which may be written for use alo a 
adiabatic in the form b=B[(p/po)*—1], and y 
is interpreted not as the heat capacity ratio but 
as the exponent of p/po in the Tait equation. 
The asymptotic form for the peak pressure of 
the shock wave with spherical symmetry js jp 
agreement with the results of the theories of 
Kirkwood and Bethe* and with the limiting 
theory of Osborne and Taylor.® The pressure, 
in this -case, decays non-acoustically as the 
slowly varying factor (logR/R,)-. 


THE IMPULSE 


The impulse J delivered by the shock wave at 
a point of fixed Euler coordinate r is 


I= p'dt (18) 


to(R) 


along a path of constant r. If the excess pressure 
p’ has a negative phase, the positive impulse is 
obtained if the time integral is extended not to 
infinity but to the time at which the excess 
pressure vanishes. The pressure-time integral can 
be expressed in reduced form in a manner 
analogous to the reduction of Eq. (10). 


I=v*pé, 


1 (- logp’ ) 
0 ot r, on 


t—to(R) 


w= (@'/pde, = 
0 


The initial slope, —1/0@, of the Euler pressure- 
time curve is expressed in terms of the Lagrange 
partial derivatives at the shock front by Eq. (6), 
and these partial derivatives were obtained as 
functions of peak pressure by the solution of 
Eqs. (4) and (12). The following expression is 
obtained for 6: 


1 Ufe p p 1 dp 
—-=—};—+4+] — -—— — —}. (20 
6 GIR [Fato+(1 "Jef. = (2 


For an exponential pressure-time curve, con- 
sistent with the peak approximation and the 
exponential energy-time curve, the reduced 
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pressure-time integral »* is equal to 1. For the 
asymptotic linear pressure- -time curve, consistent 
with the asy mptotic quadratic energy-time 
curve, we have »*=4 for the positive phase of 
the wave. As an empirical interpolation formula 
between the two values, we have employed the 


relation,” 
y*=1—} exp—[p/po]}'. (21) 


We have found that Eq. (21) leads to satisfactory 
agreement between calculated and experimental 
values of the positive impulse when used in 
conjunction with Eq. (13) for the reduced 
energy-time integral. A more detailed analysis of 
the relation between v and »* does not seem 
justified in view of the empirical nature of 


Eq. (13). 


INITIAL CONDITIONS FOR SHOCK WAVES FROM 
EXPLOSIVE SOURCES 


The two constants of integration may be 
determined from the thermodynamic properties 
of the explosive and those of its products, either 
in the Chapman-Fouget detonation state or in 
the instantaneous detonation state corresponding 
to adiabatic isometric conversion of the entire 
explosive charge to its decomposition products. 
Since the instantaneous detonation state may be 
expected to give the better average representa- 
tion of the behavior of the detonation products 
in the generation of a shock wave in an exterior 
medium, we shall employ it in the discussion of 
the initial conditions for the shock wave. 

If at the initial instant of time, a charge of 
explosive has been converted to products at a 
uniform pressure p,, a shock wave will advance 
into the exterior medium and a rarefaction wave 
will recede into the detonation products. The 
initial excess pressure 1, and particle velocity 








u%,, continuous at the boundary, are related by 





o*(pi)+u1=0, u=p/poU(p), via 
. " pte dp 
w=) ) 


where o* is the Riemann‘ function in the explo- 
sion products and U, the shock velocity in the 
exterior medium, is determined as a function of 
the pressure by the Hugoniot conditions, Eqs. (2). 
The asterisked quantities refer to the explosion 
products. The first of Eqs. (22) expresses the fact 
that the Riemann r-function‘ initially vanishes 
in the receding rarefaction wave. In order to 
solve Eqs. (22), it is necessary to use the Flugo- 
niot tables for the exterior medium, the energy 
of explosion, and the heat capacities and equi- 
librium constants for the substances making up 
the explosion products; and to evaluate o* it is 
necessary to use an equation of state for those 
products. 

In the development of the propagation equa- 
tions, the rate of energy delivery has been 
approximated by an exponential function of 
time. For shock waves in air, it may be assumed 
that the integral of this exponential function is 
equal to the total energy of explosion, since 
experimental evidence suggests that there is little 
residual energy available for second shocks. For 
shock waves in water, experimental evidence 
suggests that approximately one-half of the 
energy of explosion is delivered to the first 
shock. The initial value of the quantity D is 
readily estimated from these considerations, and 
the disadvantages of the approximate nature of 
this procedure are minimized by the circum- 
stance that except in the immediate vicinity of 
the explosive charge, the shock wave parameters 
are not very sensitive to the initial energy. 
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As in the Wigner-Seitz and Slater cellular approximation 
it is assumed that the valence electron potential is spheri- 
cally symmetrical within a lattice cell. Then the Schréd- 
inger equation becomes separable in spherical coordinates 
so that the eigenfunction can be expanded in terms which 
are products of spherical harmonics and solutions to the 
radial equation. The method differs from the previous 
ones in that: (1) eigenfunctions belonging to special wave 
vectors are first constructed by suitable linear combinations 
of spherical harmonics so as to satisfy symmetry require- 
ments and thus, essentially, including more terms in the 
expansion without increase in labor and (2) surface bound- 
ary conditions are satisfied exactly, in effect, at points 
which are more representative of the cell surface. The 
method was tested by applying it to the Shockley empty 


lattice (body centered cubic) for which the eigenvalue 
are known exactly. The lowest four eigenvalues : 
to the reduced wave vector (0,0,0) and the lowest three 
belonging to (0,0,*/a) showed errors of one percent @& 
less in energy using from two to four terms in the eigen. 
function expansion. When applied to sodium reguly 
showed that electrons in the first few Brillouin zones ap 
indeed very nearly free, having free electron energies 
within a few percent even at points near the center ang 
corners of reduced wave vector space. Furthermore og 
the boundary of reduced wave vector space in the (0,0,1) 
direction there is no energy gap between the first and 
second Brillouin zone eigenvalues since they “stick 
together” at that and equivalent points. 





INTRODUCTION 


N the W-S!? cellular method because of the 

periodicity of the electron-potential it suffices 
to focus attention upon a single cell .of the crystal 
lattice. Schrédinger one-electron solutions for 
the valence electrons are sought satisfying the 
BC (cell boundary conditions) imposed by sym- 
metry and periodicity. The interactions of va- 
lence electrons with nuclei and core electrons of 
the cell are represented in this approximation by 
an ordinary, appropriately chosen, potential V. 
Because of high symmetry and the essentially 
neutral electrical nature of surrounding cells, 
their contribution to this potential is assumed 
to be zero. 

In the interest of mathematical simplicity and 
justified’? on grounds of high symmetry the 
potential, V, is assumed to have spherical sym- 
metry within a cell. The Schrédinger equation 
then becomes separable in polar coordinates and 
~ * The method and results were presented in condensed 
form at the New York Meeting of the American Physical 


Society, January 14-15, 1944. 
** Now at the Clinton Laboratories, Oak Ridge, Ten- 


nessee. 

1A. Sommerfeld and H. Bethe, ‘‘Electronentheorie der 
Metalle,” Handbuch der Physik (Verlagsbuchhandlung, 
Julius Springer, Berlin, 1933), Vol. 24/2; F. Seitz, Modern 
Theory of Solids (McGraw-Hill Book Company, Inc., New 
York, 1940); N. F. Mott and H. Jones, Theory of Properties 
of Metals and Alloys (Oxford Press, New York, 1936). 

* E. Wigner and F. Seitz, Phys. Rev. 43, 804 (1933) and 
46, 509 (1934). 


the one-electron eigenfunction, ¥, can be ex 
pressed as 


V=ZA[EmCim¥i"($, 6) IRAE, p) —(t) 


where the Y," are SH (surface spherical har- 
monics), and the R,(E, p) are the solutions to 
the radial equation with an energy parameter, 
E. The constant coefficients Az, Cim and the 
energy parameter E can all be chosen in principle 
so as to fit the appropriate BC at the surface of 
the cell. (The split between A; and Cin is of 
course largely arbitrary but convenient for some 
of the following.) 

The labor without limit implied in (1) was 
reduced by Slater*® to practical limits by the 
approximation of assuming all Cin’s to be zero 
for / greater than a certain fixed value. Then, of 
course, the BC cannot be fitted over the entire 
surface of the LC. Slater therefore selects certain 
points on the surface of the LC at which he 
satisfies the BC by appropriate choice of the 
Cim’s which are not assumed zero. The fitting is 
done in this way for any arbitrary wave vector* 
k. However, Shockley has shown that Slater's 
procedure leads to large errors for an empty 

3 J. C. Slater, Phys. Rev. 45, 794 (1934). 

‘It is assumed that these eigenfunctions, ¥, have been 
brought into the Bloch form viz. y=ux(@)expik-@ where 
k is a reduced wave vector and ux(@) is periodic with the 


period of the crystal lattice. 
5 W. Shockley, Phys. Rev. 52, 866 (1937). 
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lattice for which the value of the energy & as a 
function of k is known. He concluded that a 
much larger number of Cim's will have to be 
used than was done in Slater’s calculations. 

It is the aim of this paper to include Cim's up 
to much larger /’s without increasing the labor. 
This can be done for eigenfunctions belonging to 
certain points in k space (reduced wave vector 
space) of especially high symmetry, i.e., points 
having the complete symmetry of the cubic 
lattice, or lying on symmetry lines or planes. 
The symmetry of the wave functions at such 
special points in k space was discussed by BSW.® 
Our method consists now of constructing, once 
and for all, combinations 2CimY:" which have 
the correct symmetry required according to 
BSW; these we call Kubic Harmonics (KH). 
Then choice of the A; in (1) serves to fit the BC; 
note that the BC need to be fitted only on a few 
of the surfaces of the LC, symmetry then makes 
them automatically satisfied on a number of 
other surfaces. By using only 3 or 4 different A,’s 
we can generally include values of / up to 6 which 
would require 49 coefficients Cim in Slater’s 
method. 

Actual calculations in the present work were 
confined to the specially symmetrical points 
k=(0,0,0) and k= (0,0,7/a) of k-space belonging 
to the body centered cubic lattice. The method 
was tested by applying it to the limiting concept 
of Shockley’s' empty lattice for which the value 
of the energy E is known. For the eigenfunctions 
tested, the results of the ELT (empty lattice 
test) were encouraging beyond expectation; the 
corresponding eigenfunctions and energies were 
therefore obtained for sodium. 


METHOD OF CONSTRUCTING EIGENFUNCTIONS: 
KUBIC HARMONICS 


For any distinct wave vector’ k, there is a 
certain group of symmetry operations on the 
coordinates which do not change k. This “group 
of k” is the full cubic symmetry group for the 
two values of k here considered, (0,0,0) and 
(0,0,7/a). The eigenfunctions can then be classi- 
fied according to the irreducible representation 


*L. Bouchaert, R. Smoluchowski and E. Wigner, Phys. 
Rev. 50, 58 (1936). 

’ For the meanings of “distinct k” and “group of k” the 
reader is referred to BSW, reference 6. 





EIGENFUNCTIONS AND EIGENVALUES IN SOLIDS 





613 


of the group of k which describes their transfor- 
mation under the symmetry operations of the 
group ; each irreducible representation represents 


” 


a “‘type.’’ For the cubic group, for instance, 
there are 10 different types (see below). For each 
type, we can find various “‘sets’’ of KH which 
transform according to the representation char- 
acteristic of the type. By a set we mean functions 
which transform into each other by the sym- 
metry operations of the group of k. A set will 
contain only functions of the same order /, and 
of course of the same type. Not all types occur 
for a given /, nor all /’s for a given type. 

It is the aim now to show how each correct 
linear combination of spherical harmonics 
YmCim Yi" belonging to a type, for each order / in 
which the type occurs, may be constructed. 
This construction proceeds in three steps: (1) 
constructing a “group table’ showing the be- . 
havior of each type under transformation of each 
element of the group of k, i.e., whether the 
element of the type is unchanged, designated by 
+, or changed in sign only, —, or whether it is 
transformed into a linear combination of other 
members of its type, 0. (2) constructing ‘“‘char- 
acteristic polynomials,’’ CP, in x, y, z for each 
type and order / which have the transformation 
properties of the group table constructed above, 
and (3) deriving the correct linear combinations 
of the SH by dividing the CP by p'=(x*+~,y" 
+2*)"2 and orthonormalizing. This process is 
carried out below for the types belonging to 
wave vectors (0,0,0) and (0,0,7/a) of the body- 
centered cubic lattice. These correct linear com- 
binations for these wave vectors are the KH 
sought. 

The group of these wave vectors is the full 
cubic symmetry group.® It consists of 48 ele- 
ments, 24 rotations forming a normal divisor, 
N, and these 24 rotations followed by inversion 
about the center forming a coset, JN. N and JN 
each contain five classes and, therefore, there 
are 10 irreducible representations and ten types 
of KH a, 8, vy, 5, €, and a’, B’, y’, 8’, é& corre- 
sponding to the five “‘positive”’ irreducible repre- 
sentations I';, I's, I's, I's, I's, and the five ‘‘nega- 
tive’ irreducible representations I';’, T's’, I';’, T's’, 
r,’, with dimensions 1, 1, 2, 3, 3, respectively. 

It is useful to know the type of KH and the 
number of sets of a type that can be constructed 
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TaBLE I. Group Table. Behavior of Kubic Harmonic types under operations of the Cubic Symmetry Group, Plys indi. 
cates invariance under the element; minus means change of sign only; zero indicates change to one of the other funcs 
with which the given function is degenerate. If the type is degenerate the behavior of that function is listed which has the 


2 axis as a unique axis. 














a 

i 2— 3 4 5 6 7 8 9 10 11 12 
Class Operation Rotation a a’ 8 8’ 7 YY’ 6 a r 2 
aan ae. —— ee 
E zs ¥ 8 none > + + + ¢ + + ae + 
C2 —x—-y 2: z through + ; + £.t ** 2 -* 22 
x—-y -—z x through x + + + + + + |= = oe 
—x yous y through r + + + + + + | = oo 
C; —y x s through +2/2 + + - - + + + + 
yrmx s z through —x/2 + + - - +- +- + + ae 
x-s y x through 2/2 + + - = 0 0 0 0 O 9 
x ey x through —2x/2 + + —=+ _ 0 0 0 0 0 9 
2 syo-x y through x/2 + + - = 0 0 0 0 6c 
—s y y through —2/2 + + - = 0 0 0 0 0 4g 
Cs y x—-s x=y through r + + - - + + - = 2 
s-y = x=z through x + + - = @ 0 0 0 0 9 
—x 3s y y=s through r + + - = 0 0 0 0 0 9 
—y -x —s x= —y through r + + - - + $+ © = 
—s —-y —x x= —g through r + + - = 0 0 0 0 6a 
—x —-s -y y=-—z through + + + —- = 0 0 0 0 0 9 
Cs ss 9 x=y=<s through+272/3 + + + + O 0 0 0 0 9 
y ss 2 x=y=zs through —2zx/3 + + + + 0 0 0 0 0 90 
s—x —y x=—y=z through +22/3 + + + + O 0 0 ®o tan 
—y-s « x=—y=z through —2x/3 + + + + 0 0 0 0 oO 9 
—s-x y x= —y=-—zs through 22/3 + + + + 0 0 0 0 8 @ 
—y s-x x= —y=—z through —27/3 + + + + 0 0 0 0 oO 9 
—z «x -y x=y=—z through 27/3 + + + + 0 0 0 0 0 90 
yr-s-x x=y=—s through —27/3 + + + + 0 0 0 0 6 
J —x -y —2 none + - -—- + + = - + + -= 
JC: x yous z through x t+ - - + + - - + 559 
—< y 3 x through x + - + + - + —- = & 
JC; y-x —s s through +2/2 + + - + =- + - 385 
—x s-y x through +2/2 + —- + -=- 0 0 0 0 0 0 
JIC, —-y-x 32 x=y through r + - + - + - + - + & 
—s yx x=z through r + - + - 0 0 0 0 0 0 
JCs —s —-x -y x=y=z through 27/3 + - — + 0 0 0 0 0 0 








for each order /. One of us® has derived this and 
the results are: for: /=0, one a type; /=1, one 6; 
l=2, one each of y and e; /=3, one each of 
B, 6, e’; l=4, one each of a, y, ¢, 6’; 1=5, one 
each of y’, e’ and two 5; ]/=6, one each of a, y, 
B’, 5’ and two e. 


Construction of the Group Table 


The behavior of each type of KH under the 
transformations of the crystal group is shown 
(for one function of a set of each type) in Table I. 
This is obtained by a study of the lower order 
SH using the above results as a guide. For /=0 


8H. Bethe, Ann, d. Physik 3, 133 (1929). 


there is the non-degenerate type a and there is 
one SH which is spherically symmetrical and 
remains invariant under ‘all transformations. 
The behavior of this type is shown in column 4. 
For /=1 there is the triply-degenerate 6 type. 
The three SH are x/p, y/p, z/p. The latter one 
transforms according to column 10, and only 
into either itself or one of the remaining two SH. 
For /=2 there are the doubly degenerate y type 
and the triply degenerate « type. The SH are 
[2*—1/2(x*+-y*) ]/p?; (y? —x*)/p?; 2x/p?, zy/p* and 
xy/p*?. The behavior of the first of these is shown 
in column 8. In those cases in which it does not 
transform into itself it transforms only into a 
linear combination of itself and the second listed 
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above. It is therefore double degenerate and is 
identified with the y type. The behavior of 
xy/p? is shown in column 12; it is degenerate 
with the remaining two, and is identified with 
the « type. For /=3, there is among others the 
non-degenerate 6 type. Among the SH of /=3 
there is found the function xyz/p*. This obviously 
satisfies the requirements for a non-degenerate 
function and its behavior is listed in column 6. 
The remainder of the table corresponding to the 
behavior of the a’, 8’, y’, 4’, e’ types can now be 
determined with ease by observing that their 
behavior differs only in sign from corresponding 
unprimed types whenever an inversion is in- 


volved. 


Construction of Characteristic Polynomials 


For this purpose it is noted that any power of 
p is invariant under the transformation of each 
element of the cubic group. Thus not only the 
SH x/p, y/p, 2/p form functions of type 6 but 
also the first order polynomials x, y, z or the third 
order polynomials xp*, yp’, zp*. All sets which 
differ only by a power of p are therefore con- 
sidered “identical” for the present purposes. 

The problem of constructing CP is the problem 
of finding linear combinations of terms of the 
form x?y%" (l= p+ q+r) which are linearly inde- 
pendent of other chosen polynomials of the same 
order (or identical polynomials of lower order) 
and which behave according to their type (Table 
1). The results of reference 8 quoted above aids 
in pointing out types for which a search is 
instituted. 

As an example, the procedure for obtaining 
the CP of order four is briefly outlined, assuming 
CP of lower order have already been obtained. 
This order contains one set each of types a, y, € 
and 6’. These are nine linearly independent func- 
tions to be constructed from the fifteen linearly 
independent polynomials: three of form x‘, six of 
form xy’, three of form x*yz, and three of form 
x*y*, However, these fifteen are effectively re- 
duced to nine because of the six fourth order 
polynomials obtained by multiplying identical 
lower order polynomials by powers of p? viz: 


a type, P:: p', 


y type, Pe: p?(x*—y*), p*(y?—2*), 
e type, P3: p*xy, p*yz, p*zx. 
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With an eye on the Table I there are constructed 
the following : 


Polynomials obviously of the a type 


Py: x*+y'+2! 
Ps: x*y’+y"2?+2*x? is not new but is equal 
to $(P,—P,) 


Polynomials of y type 


P,: x4—y, y—2 
Pz: (x*—y*)z*, (y?—2*)x* not new viz. P2—P, 


Polynomials of ¢ type 


Ps: xyz", y2x*, sxy" 
Py: *y+xy', y's+yz2", 2x+2x* not new vis. 
P;—P3 


Polynomials of the 6’ type 
Pio: *y—xy', y’z—y2*, 28x —2x* 


These exhaust the possible number of fourth 
order linearly independent functions. The three 
fourth powers are exhausted in the three function 
P, and Px, the six of form x*y in P, and Pip, 
the three of form x*yz in Ps, and the three of 
form x*y* in Ps and P;. 


Construction of Kubic Harmonics from the 
Characteristic Polynomials 


Dividing the CP of order / by p' yields func- 
tions, of the angles alone. These functions, Fy, 
have the correct symmetry properties of type s 
but do not obey Legendre’s differential equation 
of order 7. However, from their nature it follows 
that they are linear combinations of SH of order 
l and lower. In order to eliminate all lower order 
harmonics F;, can be made orthogonal to all SH 
of order less than /. Because of the orthogonality 
of functions belonging to different irreducible 
representations, it is only necessary to make F;, 
orthogonal to all KH of type s and lower /. This 
is done by expanding F;, in terms of KH of the 
same type and orders / and lower, and obtaining 
the KH, Ky, in the usual way. One of us® has 
obtained all the KH up to /=6; these are listed 
in Table I]. 


*H. A. Bethe, unpublished paper (1935). 














_ TaBeE II. List of the Kubic Harmonics classified according to symmetry properties. All functions normalized to de, 


Factors of p~' are omitted throughout, i.e., x, y, 2 are written for x/p, y/p, z/p. Functions in square brackets indicate 
functions with normalization factors omitted. 
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t For types 3 and ¢, which are triply degenerate, the other two functions are found by cyclic interchange of coordinates. 














BOUNDARY CONDITIONS AT THE SURFACE OF 
THE CRYSTAL CELL 


The BC on the surfaces of the crystal cell are 
fixed by the periodicity of the eigenfunction 
which in turn is the property of its wave vector. 
A translation from a point A on one face to the 
point B perpendicularly opposite on a parallel 
second face is a translation of a lattice-point 
vector T. It is well known® that y(e+T) 
=¥(9)expik-T so that, in general 


ve =vaexptk-T. (2) 


These BC can be put into a generally more usable 
form by employing the symmetry properties 
(listed in a group table) of the eigenfunction 
type. For the KH this is now done. 

The unit crystal cell for the body centered 
cubic lattice is shown in Fig. 1. It has two non- 
equivalent face forms: six square faces on planes 
x=+a, y=+a, z=+a; and eight hexagonal 
faces whose normal vectors are (+a/2, +a/2, 
+a/2). The translations between pairs of square 
faces are (0,0,2a), (0,2a,0), and (2a,0,0) and those 
between the hexagonal faces are (+-a,+a,a). For 
the wave vectors (0,0,0) and (0,0,7/a) (to which 
the KH belong) the periodicity conditions are 
particularly simple. From (2) they are: 


fork=0,0,0 Eigenfunctions are periodic 
with the periods of all 
pairs of parallel faces of the 
cell 
Eigenfunctions are periodic 
with the period of pairs of 
square faces 

Eigenfunctions are anti- 
periodic with the period of 
pairs of hexagonal faces 


(3a) 


for k=0,0,7/a 
(3b) 


(By “anti-periodic” is here meant that the func- 
tion changes sign only upon a displacement by 
the “period” stated.) Correspondingly, eigen- 
functions belonging to k=0,0,0 are here said to 
be periodic ; those belonging to k =0,0,x/a, anti- 
periodic. 

As an example, the working BC are now de- 
rived for the anti-periodic 4-type function with 
2 axis preferred. Non-equivalent pairs of faces 
are considered separately. For the square faces 
= +a: Periodicity, (3a), and symmetry oper- 
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(3c) . 


ator x, y, —2 (Table I), give successively 
v(x, y; a) =y7(x, y; —a) = —y(x, y, a) 


so that the eigenfunction vanishes over this face. 
For the square faces x = +a (or y= +a): Perio- 
dicity and symmetry operator —x, y, z gives 


V(at+e, y, 2) =¥(—a+t+e, y, 3) =¥(a—e, y, 2) 


so that the normal derivative dy(a, y, z)/dx, 
vanishes over this surface. For the hexagonal 
surfaces +x+y+z=+(3/2)a: Anti-periodicity, 
(3c), and symmetry operator —y, —x, —z gives 


v(x, y, 2) = —¥(x—a, y—a, s—a) 

= +y¥(a—y, a—x, a—z) 
so that the function is symmetrical about the 
line ¢=0 (Fig. 1). Similarly, the normal deriva- 
tive is anti-symmetric about the same line. 
Table III lists working BC similarly obtained 
for all types of functions for that function of a 
type which has its z axis preferred. 


Fitting of the Boundary Conditions— 
Empty Lattice Test 


In the notation of Kubic Harmonics the candi- 
dates for eigenfunctions (1), have the form 


Ve = LwAuK wRiALZ, p), 


z — a? 


(4) 











Fic. 1. Unit lattice cell for the cubic body centered 
crystal lattice. The lattice constant, 2a, for sodium, is 
shown as 8.138 Bohr radii. The circle on the hexagonal face 
is the trace of an equal volume sphere and it was over this 
—_ a the boundary conditions were approximately 
satished. 
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TaBLE III, Boundary conditions on y¥ and its normal 
derivative on the faces of the unit cell. For nerate 
types the BC on that eigenfunction is listed which has its 
z axis preferred. 0 = function vanishes over surface ; a = anti- 
symmetric; s = symmetric. 











Square faces Hexagonal faces 
Either lodic or anti- S abou’ 
( ; : ome t line 
x=tea Anti- 
’ Type | (or y=+a) s=+¢ Periodic 
oy oy oy oy 
_-— © ee @ *_% 
a 0 0 $s a a s 
B 0 0 a s s a 
Y 0 0 $s a a $s 
r 0 0 a s s a 
€ 0 0 s a a s 
a’ 0 0 $s a a $s 
p’ 0 0 a s s a 
7’ 0 0 $s a a $s 
3’ 0 0 a s $ a 
é 0 0 $s a a s 














where / is summed over those orders containing 
the type s, and ¢ is summed over the multiple 
sets of a type which may occur for the higher 
orders of /. Because the Kubic Harmonics contain 
all the requirements imposed by the symmetry 
of the lattice, (4) already satisfy proper BC at 
many points on the surfaces of the cell. It 
remains to satisfy the remaining BC imposed by 
periodicity shown in Table III. 

The radial functions R; can in general only be 
obtained by numerical integration (see next 
section) and depend on the special potential V 
used in the Schrédinger equation. Obviously, 
the coefficients A, will also depend on the 
special problem and cannot be determined ana- 
lytically. To keep the labor within limits, only 
a finite number of terms in (4) must be used. 
Actually at most, four terms were used in which 
case three” constant coefficients Ay, and the 
eigenvalue E can be determined. This allows, in 
principle, the BC on at most" four distinct (not 
symmetrically equivalent) points to be satisfied, 
aside from those points whose BC are automati- 
cally satisfied by the KH. Eigenvalues obtained 
by various choices of sets of points yielded results 
10-100 percent in error” for an empty lattice. 


” The fourth coefficient is determined by normalization. 

1 Most points on the hexagonal faces have BC on both 
y and its normal derivative so that if such a point were 
included the number of points that can be fitted would be 
reduced. 

121t is noteworthy that choices including the central 
point of the hexagonal faces used by Slater gave results 


It became evident that some method Must be 
developed to produce a better “average” fit, 
This is now described. 

To simplify arithmetic and to put (4) ina 
generally more convenient form, polynomials 
related to the Kubic Harmonics are defined 


Pu=p'Kiu 
so that (4) becomes 


¥=ZwAuPi(x, y, 2)Si(E, p), (5) 
with 
S:=Rip". (Sa) 


The functions S,; are obtained by mechanigl 
integration. The form (5) or its normal derivatiy. 
is used directly to fit BC on chosen points of, 
square face (actually at most only one point m 
squares, a central point, was fitted). 

On the hexagonal surfaces an “average” fit 
was desired. For this purpose, a more complicated 
method had to be used: Eq. (5) was transformed 
to cylindrical coordinates, 7, a, § (a=0 in dire. 
tion +n, Fig. 1, and & normal to the surface) 
Each polynomial, P;,, was expressed as a Fourier 
series in a whose coefficients, Pim, are homo 
geneous polynomials in £ and r of order 1. The 
result is 


¥=[2uAwSi(E, p)ZaP itn cosna | 
+[21A wSiZmPiim sinma] (6 


Wherein either ” is summed over zero, and al 
positive even integers, and m over all positive 
odd integers or the reverse, depending upon the 
type. In either case, the first term in (6) 
symmetrical, the second is anti-symmetrical 
about the axis a=0 and correspondingly pis 
said to be composed of an even part, ¥,, and a 
odd part, y_, respectively. The BC on the hem 
gonal surfaces are completely satisfied when % 
(or y_) and the normal derivative of y_ (or ¥) 
vanish over the hexagonal surfaces corresponding 
to the BC (Table III) that y is antisymmette 
(or symmetric). 

For a symmetrical eigenfunction the exat 
boundary conditions on the hexagonal surface 
may be stated in the form. 


which were particularly bad. This is understandable 
because the central point has a much smaller distance from 
the origin than an average point, and also has a mua 
higher symmetry which makes the BC degenerate. 
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For all m ZuAwSilE, p)Pum(r, —) =0 (7) 
and for alln (0/0§)[2uA w5rP itn] =0, 


respectively. (When y is anti-symmetric m and n 
are interchanged in (7).) For fitting with the 
limited number of Ax, it was assumed, as an 
approximation that the radius, r,, of the trace 
of the equal volume sphere on the hexagonal 
surface was a good “‘average”’ radius and that 
boundary requirements (7) corresponding to the 
higher harmonics in the expansion (6) could be 
neglected. The actual numerical procedure was, 
using the experimental lattice constant, Fig. 1, 
to make a likely choice of E and determine the 
disposable constant coefficients, A, with a com- 
patible number of BC, (7), of lowest m and n. 
The choice of E is then tested by use of a BC on 
a square surface, and then the procedure is 
repeated until the consistent value for £ is found. 

The method was tested using the ELT for 
each of the solutions carried out. In the limit of 
vanishing potential V, the radial functions S; 
are functions akin to the Bessel’s functions 
Jis(K, p)/p'. Using these functions, the eigen- 
values E are obtained from the above procedure 


of fitting the BC, and can then be compared 
with their known true values, K?. (K is the “‘free 
electron” wave vector which depends on k and 
on the number of the Brillouin zone.) The results 
are given in Table IV. 


Discussion of the Method 


Table IV lists the error in the computed empty 
lattice eigenvalues belonging to k=0,0,0 and 
0,0,*/a as a percentage of their true values for 
the lowest lying levels (not including the ground 
state for which the test is not applicable since 
¥=constant is an exact solution). The number 
of terms used in the expansion of each of the 
eigenfunctions is also listed. This was arbitrarily 
limited by including Kubic Harmonics only up to 
1=6, 

The results of the empty lattice test are most 
gratifying. At the upper end of the first Brillouin 
zone (k=0,0,7/a), the error in the energy is one 
quarter of a percent or less, as compared with 
errors up to 40 percent in Slater’s method (see 
Fig. 2, reference 5). Even at the upper end of 
the second zone, the error does not exceed 1 
percent (Slater’s method 35 percent). This im- 
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provement has been achieved while taking, in 
each individual case, only 2 to 4 terms in the 
expansion while Shockley took 8 terms in testing 
Slater’s method. Even considering the greater 
complication of our Kubic Harmonics as compared 
with Legendre polynomials, and of our average 
fitting of the BC as compared with Slater's 
fitting at a point, the total labor in our method 
is probably no greater than in his. In fact, it 
would not involve excessive labor to include a 
considerably greater number of terms in our 
expansion. This would permit an extension of 
the method to higher Brillouin zones. 

Since the eigenvalues for the empty lattice are 
in error by one percent or less it seems likely 
that the accuracy for filled cubic lattices will be 
limited by the accuracy of the potential, V. It 
is believed, therefore, that, at least for the wave 
vectors here considered, this is a one-electron 
approximation for cubic solids which compares 
favorably with methods for treating atoms. 

The method is not confined to the investigation 
of one-electron states belonging to k=0,0,0 and 
0,0,7/a. It appears to be a straightforward 
problem to obtain the energy spectrum for other 
symmetrical wave vectors without recalculation 
of the radial functions. This would involve con- 
structing appropriate harmonics corresponding 
to the new irreducible representations and apply- 
ing the more general periodicity boundary con- 
ditions compatible with the new wave vectors. 
Once the energy spectra of a number of sym- 
metrical wave vectors has been determined the 
eigenvalue hypersurfaces can be traced out by 
interpolation guided by the BSW compatibility 
tables.* Behavior near wave vectors of symmetry 


TABLE IV. Percent error in eigenvalues for the empty 
lattice and eigenvalues for sodium belonging to the lowest 
lying levels for k=0,0,0 and 0,0,x/a. Fitting for a and 6 
types belonging to k=0,0,7/a by the authors; remainder 
by Bowers (see reference 13). 
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can be investigated by perturbation methods 
with eigenfunctions obtained above; this has 
been done by Bowers.” 

On the other hand, the harmonics and working 
relations are the same for all crystals having the 
same space group. Thus, having done the work 
for one substance, the problem for a substance 
having an identical lattice structure becomes one 
only of obtaining new radial functions and solv- 
ing elementary sets of simultaneous algebraic 
equation. It must be admitted that thus far the 
method has been tested only for wave vectors of 


1%W. Bowers, Doctoral Thesis, Cornell University 


(1943). 
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; Fic. 2. This shows the 
eigenvalues at k=(0 
and (0,0,x/a) as 
for sodium by the 
method. For compar} 
the eigenvalues in the 
0,0,1 direction of the first 
two Brillouin zones fo 
sodium as computed by 
Slater and for the free 
electron are shown. Detail 
1 and Detail 2 sketch the 
BZ in the (0,0,1) direction 
which join up withthe 
computed values (see 
BS compatibility 
tables). Note that due to 
the degeneracy of the 
y type at k=(0,0,2/a) 
there is no gap in the 
(0,0,1) direction between 
the first and second BZ, 














high symmetry. Accuracy of the method in 
general for crystals having a cell shape and 
symmetry departing widely from that of a sphere 
remains at present in doubt because of the basic 
approximation of separability of variables in the 
Schrédinger Equation. 
RESULTS OF THE METHOD AS APPLIED 
TO SODIUM 


When the method is to be applied to any 
problem other than the empty lattice, the po 
tential V must be appropriately chosen. We 
have used the same potential’ energy function 


“4 W. Prokofjew, Zeits. f. Physik 58, 255 (1929). 
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V as W-S' and Slater.* The radial functions, R:, 

were obtained by mechanical integrationft and 
checked by numerical integration using I.B.M. 
machines. The agreement was satisfactorily 
within the expected limits of error of the differ- 
ential analyzer. Radial wave functions were ob- 
tained for /=0 to 6, and for energies E= +0.6, 
+0.3, +0.1,0, —0.1, —0.3, and —0.6 Rydbergs. 
The eigenvalues were obtained by interpolation. 

Table IV lists the seven lowest lying eigen- 
values for k=0,0,0 and 0,0,7/a as computed for 
sodium. For comparison Fig. 2 shows in addition 
the eigenvalues of the first and second Brillouin 
zones along the line 0,0,1 in k space, for the case 
of the free electron and that of sodium according 
to Slater’s Fig. 3. Detailed sketches in Fig. 2 
exhibit (for sodium) the variation from free 
electron eigenvalues near the center and bound- 
ary of the Brillouin zones. They also show the 
“sticking together” of the eigenvalue hyper- 
surfaces at these points which was pointed out 
by BSW, and indicate the extremity eigenvalues 
to which the lower hypersurfaces join as deter- 
mined by the BSW compatibility tables. The 
notations A, As, etc., indicate that the given 
eigenvalue belongs to the respective representa- 
tion of the group of k (notation of BSW). 

A very significant result of the present calcu- 
lation is that for sodium in the 0,0,1 direction 
the electrons are essentially free. If the electrons 
were in fact free, the eigenvalue at the corner 
(0,0,7/a) of the first Brillouin zone would be 
—0.608+ (2/a)*? = —0.012 rydberg which is very 
close to the eigenvalue for the 5-type (—0.0135) 
actually found there and is in between the other 
two types. However, there is a slight deviation 
from freeness very close to this corner since the 
lowest lying energy hypersurface joins with the 
y type with energy —0.036 rydberg. (The total 
width of the first Brillouin zone in the 0,0,1 
direction is thus 4 percent less than for free 
electrons.) Bowers" has shown that this devia- 
tion extends only a very short distance from the 
corner. This essential freeness of the electron up 

tt In this connection the authors are indebted to Pro- 
fessor S. H. Caldwell of M.I.T. in making available the 
use of the mechanical Bush Differential Analyzer; to Mr. 
P. O. Crawford who as a graduate assistant (1940) made 
the work possible by operating the machine; to Dr. Martin 
Schwarzschild who ran an independent check on two of the 


functions at the Thomas Watson Computing Bureau at 
Columbia University. 


EIGENFUNCTIONS AND EIGENVALUES IN SOLIDS 621 






to a point close to the corner of the first Brillouin 
zone is in agreement with the calculations of 
Bardeen,'* wherein he obtained (near k=0,0,0) 
the ratio of electronic mass to effective mass of 
1.069 indicating a deviation in eigenvalues of 
7 percent from those of the free electron. 
Actually, as shown in Fig. 2, the eigenvalues 
of the higher Brillouin zones also remain close 
to those for the free electron. However, there 
can be little doubt that the gap between the 6 
and ¥ states amounting to 0.022 rydberg is real 
because the empty lattice test indicates an uncer- 
tainty in calculated energies of only 0.26 percent 
of 0.6 Ry, that is, 0.0016 Ry, about 7 percent of 
the calculated difference between 6 and y states. 
Figure 2 shows poor agreement with Slater’s 
calculations for the 0,0,1 direction. It is sus- 
pected that these markedly different results are, 
in part, because of his choice of the centers of 
hexagons as fitting points. Experience gained in 
the present work indicates that these points were 
far from average points yielding particularly 
poor results in the case of the empty lattice.” 
At k=0,0,7/a the lowest lying eigenvalue is 
the y type (Z, = —0.036 Ry) considerably below 
the lowest a type (E.=+0.094 Ry). The type 
A; (notation is that of BSW), belonging to the 
wave vectors lying on the line 0,0,1 is compatible 
with both the a and y types at the line ex- 
tremities. Thus, as we go along the line 0,0,1 in 
k space, the eigenvalues of the first Brillouin 
zone start at the a type ground state and end on 
this y type at k=0,0,7/a. Since the y type is 
doubly degenerate it is also the lowest eigenvalue 
of the second Brillouin zone; if we then proceed 
again along the line 0,0,1 this eigenvalue emerges 
as a A: type from this y type and ends on the 
low lying y type at k=0,0,0. Therefore, due to 
the “sticking together” effect of BSW, there is 
actually no gap in energy at all in the direction 
0,0,1 between the eigenvalues of the first and 
second Brillouin zones. This condition of course 
does allow a change in density of states near the 
corner, and probably also gaps in other direction 
of k space, but at the corner of the first Brillouin 
zone there is no actual discontinuity. 
No experimental verification of sufficient re- 
solving power exists to support these results. 


4 J. Bardeen, J. Chem. Phys. 6, 367 (1938). 
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X-ray absorption measurements of high resolving 
power for potassium, which is expected to show 
wider departure from free electrons, have been 
carried out by Platt.'* The K edge investigated 
shows quite close agreement with his theoreti- 
cally predicted absorption which was based on 
the assumption that the electrons are free. No 
evidence existed to show an energy gap. This, 
of course, is at best supporting evidence of the 
non-existence of a gap in potassium since gaps 


16 J. B. Platt, Phys. Rev. 69, 337 (1946). 


may exist which could be completely masked by 
the eigenvalue dependence upon wave Vector 
direction. 
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The Band Theory of Graphite 
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The structure of the electronic energy bands and Brillouin zones for graphite is developed 
using the “tight binding’ approximation. Graphite is found to be a semi-conductor with zero 
activation energy, i.e., there are no free electrons at zero temperature, but they are created 
at higher temperatures by excitation to a band contiguous to the highest one which is normally 
filled. The electrical conductivity is treated with assumptions about the mean free path. It is 
found to be about 100 times as great parallel to as across crystal planes. A large and anisotropic 
diamagnetic susceptibility is predicted for the conduction electrons; this is greatest for fields 
across the layers. The volume optical absorption is accounted for. 


1. INTRODUCTION 


HE purpose of this paper is to develop a 

basis for the explanation of some of the 
physical properties of graphite through the band 
theory of solids. We shall be concerned pri- 
marily with a discussion of its electrical con- 
ductivity, but the treatment given makes pos- 
sible the explanation not only of the electrical 
conductivity and its anisotropy but also the 
thermal conductivity, diamagnetic susceptibility, 
and optical absorption. 

The electrical resistivity of single crystals of 
graphite is about 4 to 610-5 ohm-cm.' This 
corresponds to a conductivity of the order of 
that of a poor metal. The temperature coefficient 
of the conductivity is negative, as in the case of 


* Now at McGill Serenity. 
1 Given by E. Ryschewitsch, Zeits. f. Elektrochem. ang. 
physik. Chemie 29, 474 (1923), as 3.9-6X 10-* ohm-cm. 


a metal. Polycrystalline graphite, on the other 
hand, has a much higher resistivity which varies 
very strongly according to the type of graphite 
used, and has a positive temperature coefficient 
of conductivity? to about 1400°C, and negative 
thereafter. Since the crystals of commercial 
graphites tend to be of the order of 10~* cm, and 
it is quite porous (density ~1.6 as against 2.25 
for single crystals), it seems reasonable to 
attribute the high resistivity of polycrystalline 
graphite to the crystal boundaries, on which may 
be lodged impurity atoms. The latter would tend 
to be driven off on heating, thus accounting for 
the observed temperature dependence. We shall 
show, however, that the band theory would 
seem to make possible the explanation of the 
conductivity properties of single crystals. 


2C. A. Hansen, Trans. Am. Electrochem. Soc. 16, 329 
(1909) gives 137.5 10-* at 0°C 82.5 10-* at 1400°C. 
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2, ZONE STRUCTURE OF A SINGLE 
HEXAGONAL LAYER 


Since the spacing of the lattice planes of 
graphite is large (3.37A) compared with the 
hexagonal spacing in the layer (1.42A), a first 
approximation in the treatment of graphite may 
be obtained by neglecting the interactions be- 
tween planes, and supposing that conduction 
takes place only in layers. 

Graphite possesses four valence electrons. 
Three of these form tight bonds with neighboring 
atoms in the plane. Their wave functions are of 


the form 


1 
yg Mel) +v2y.(0:2p)) (¢=1, 2, 3), 


where y(2s) is the (2s) wave function for carbon 
and y.(o;2p) are the (2p) wave functions whose 
axes are in the directions o; joining the graphite 
atom to its three neighbors in the plane. The 
fourth electron is considered to be in the 29, 
state, its nodal plane being the lattice plane and 
its axis of symmetry perpendicular to it. The 
three electrons forming co-planar bonds will not 
play a part in the conductivity ; we shall therefore 
treat graphite as having one conduction electron, 
in the 2, state. 

For the hexagonal layer the unit cell, which is 
designated by WX YZ in Fig. 1, contains two 
atoms, A and B. The fundamental lattice dis- 
placements are a@,;=AA’ and a,=AA”; their 
magnitude is 1.42 X v3 = 2.46A =a. The reciprocal 
lattice vectors have magnitude 2/v3a, and are 
in the direction AB and AZ, respectively. From 
this it follows that the first Brillouin zone is a 
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hexagon (Fig. 2) whose sides are distant 1/v3a 
from its center. It is easily shown that this zone 
contains one electron per atom, for the density 
of electron states in k-space is 2A, where A is 
the area of the crystal. The zone in question 
therefore contains 2A X2/(v3a") electron states. 
But the atomic area (area per atom in the layer) 
is V3a?/4, and the number of atoms is 44 /v3a?, 
which is exactly equal to the number of electrons 
in the zone. 

The zone containing two states per atom is 
obtained by extending the sides of the hexagon 
to form a six pointed star. 

Let us now consider the problem from the 
viewpoint of the “tight binding approximation."’ 
If X(r) is the normalized orbital 2p, wave func- 
tion for an isolated atom, then the wave function 
in the tight binding approximation has the form: 


v= gitr¢e (2.1) 
where 


¢1= >a exp[ 2mik- rg |X (r—ra) 
¢2= >> xp exp[ 2rik-rpg |X (r—rz) ; 


The first sum is taken over A and all the 
lattice points generated from it by primitive 
lattice translations; the second sum is similarly 
over the points generated from B. 

From variational principles it is known that 
to get the best value of EZ for this approximation 
we substitute in the wave equation from (2.1), 
multiply by ¢: and ge, respectively, integrate, 
and then eliminate \ from the resulting two 
equations. 

Let us neglect the overlap of the p, wave 
functions centered on different atoms, i.e., let us 


(2.2) 


and 
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fxe-twxe —rg)dr=0. 


Then, substituting (2.1) in 
Hy = Ey 
and proceeding as outlined, we get 
Ay+dHn=ES, 
Ha+hHo2=dES, 


where 


Hu= f 6:*Hexdr, Hu=Hx*= { $s*Héadr, 


A2= fon godt 


S= f ostordr= f os%6sdr. 


Eliminating \ we obtain the secular equation 


Hy—ES Aix 
An H2—ES 


from which it is found that 


=0, 


1 
Bo iutHe 
+ ((Hi1— H22)*+4|Ais| 2)$}. (2.5) 


Now by virtue of the neglect of the overlap 
integrals, S= N, the number of unit cells in the 
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crystal, and by symmetry Hi: =H2:. Introducing 


1 1 


Ay)’ =A’ =—H,=—Hx, (2.6) 
N N 


we have, finally 
E=Hy,;'+ | Hi" |. (2.7) 


The positive sign will apply to the outside of the 
hexagonal zone, the negative sign to the inside. 
The discontinuity of energy across the zone 
boundary is then 


AE=2| Hi,’ |. (2.8) 
Let us now calculate Hy;’=H2.' and Ay’. 
1 * 
A, =— 7 exp[ — 2mik-(ra—ra’) | 
A, A! 


 f X*(e— 14) HX rue 


Omitting the exchange integrals over atoms 
more distant than nearest neighbors amongst 
the atoms A, and writing 


Suni f X*(n)HX (pdr, 


yo = — [ x*e-e)HX (ar, 


where 0’=a; (say) is a vector joining nearest 
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neighbors among the atoms A, we get 


Hy’ =Eo—2y0' (cos2rkya+2 cosrk.av3 cosrk,a). 
If we write 

H=H,+(H-A)), 
where Hy is the Hamiltonian for an isolated 
carbon atom, and put 

H-—-Hy=V-—U<0O, 
where V is the periodic potential of the lattice 


and U is the potential field for an isolated atom, 
we may write, since HpX = EyX, 


E,=£- f X*()(U-V)X(pdr, (2.9) 


and 
site f X*(r—p’)(U—V)X(x)dr>0. (2.10) 


E, is the energy of an electron on the 2, state in 
carbon. 

Let us next calculate Hj,. For this, we shall 
consider only interactions between nearest 
neighbors in the lattice, the nearest neighbors of 
atoms of type A being always atoms of type B 
and vice versa. Writing 


w= f X%-0)(U- V)X(r)dr>0, (2.11) 


where 9 = AB, we obtain 
Hy’ = — yo(exp[ — 2rik.(a/v3) ] 


+2 cosrk,a-exp[2rik.(a/2v3) ]), 
so that 


| Hi2|? = yo?(1+4 cos*rk,a +4 cosrk,a cosrk,V3a). 


The energies at various points may now be 
written down: 


at O: 

at D: 

at C: 

at F, inside: 
outside : 


E=Eo—3y0—6y0, 
E=Eot+3y0—6y0'; 
E=Eot+3y ; 
E=Eo—yot2y0', 
E=Eotyot2y0’. 


Over a side of the zone there is across the bound- 
ary at any point a discontinuity of energy of 
amount 


2y0(2 cosrkya— 1) (2.12) 
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which is greatest at the center and decreases to 
zero at the corners. 

The degeneracy at C and similar points, and 
the zero-energy gap at these points, are pre- 
sumably consequences of the symmetry of the 
lattice, and are independent of the approxima- 
tions considered. 

It is easy to determine the form of the wave 
functions corresponding to different k. Consider, 
for example, the point F. Inside the zone 
boundary \= —1, and consequently the nodes 
of the wave functions bisect the lines joining A 
and B atoms: the sign of the wave function at 
lattice points is designated by + or — signs in 
Fig. 3. It is readily seen how the expression for 
the energy arises. Since an atom has two nearest 
neighbors of the same sign, and one opposite, and 
since (_V—U) is negative, there will be a con- 
tribution to the energy from exchange between 
nearest neighbors of —2y9+70= — ‘yo. Of second 
neighbors, we have two of the same sign and four 
opposite, giving us 20’. 

In Fig. 4 are represented the nodes of the wave 
function corresponding to a point outside the 
zone boundary at A. Since here \=1, the maxima 
and minima correspond to the nodes of the 
previous case. The expression for the energy can 
be verified in the same way. Similar arguments 
may be carried out for other points in k-space. 

Consider next the energy contours. Near the 
bottom of the band they are circular: 


E=Ey—3y0—6y0' 


+4 (yot6yo') (ke? +hk,*)a*?. (2.13) 


The hexagon whose sides are distant 1/2a from 
the center of the Brillouin zone is a surface of 
constant energy. (This conclusion is, of course, 
true only in the approximation which we are 
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here considering.) Curves of constant energy 
are indicated in Fig. 5. Near the corner C 


E=Ey+3yo'+V3ry0|k—k.|a—3x%y0' |k—k,|*a*, 
(2.14) 


and the surfaces of constant energy are again 
circular. 

It should be noted that, when the reduced zone 
scheme is used, the energy contours have the 
same shape for the second zone as for the first. 


3. NUMBER OF FREE ELECTRONS AND CONDUC- 
TIVITY OF A SINGLE HEXAGONAL LAYER 


In this section we shall neglect yo’ relative to 
vo. We then have 


|E—E,| ~vV3ryoa|k—k.| (3.1) 


near the corners of the Brillouin zone. At the 
absolute zero of temperature, the hexagonal zone 
will be completely filled, and the next zone will 
be completely empty. At higher temperatures, 
there will be some thermal ‘overflow’ into the 
outer zone. The extent of this overflow will be 
(in energy) of amount comparable to kT. The 
energy contour E=E,.+kT is approximately 
circular with radius k7/(xvV3y0oa) provided 
kT> vo, which is certainly the case at ordinary 
temperatures. 

Now if N(E)dE is the number of electronic 
energy states in the energy interval dE, 


N(E) =2A _— (3.2) 


|grad,# 

where A is the area of the lattice, and the 
integral is taken over the curve on which the 
energy is E. This gives, since |grad,Z| 


constant, 
NZ) ——————, 


3ryo7a? 


(3.3) 


or, since the area per atom is V3a"/4, 


1 
N(E)/Noa= |E—E£,| = ~ , 
avayo? av3y0" 


(3.4) 








where N, is the number of atoms in the lattice. 
The form of N(Z) near E. is illustrated in Fig. 6. 

Conduction will take place through the elec- 
trons excited into the upper band, and through 


R. WALLACE 





“ 








Fic. 6. 


the positive holes (equal in number) created ip 
the lower band. 

To calculate the number of these, we must 
know the Fermi distribution. For moderate tem- 
peratures, N(E) is even in e=|E—E,| over the 
whole range in which the Fermi distribution jg 
different from its value at absolute zero. There. 
fore, in the Fermi distribution function 





f(/)= : 
exp[(E—$)/kT]+1 
we shall have ¢=£,; in other words, we may 


write 


f(E) =f) =1/(e"*? +1). (3.5) 


It follows that the number of free electrons plus 
positive holes per atom is 


of “Ede =" (—) (3.6 
o N 6v3 Yo 


Now C. A. Coulson’ has estimated that 7 is 
about 20 kcal./mole, or about 0.9 ev. At room 
temperature kT =0.025 ev. Therefore the “effec- 
tive number of free electrons” mer, per atom, is 


Nett = 2.3 x 10-4. 


To determine the conductivity of our graphite 
layer, let us calculate directly the current in the 
presence of a given external field. Explicitly, this 
is given by 





j=2 f evf(k)drs-(2/c), (3.0 


where f(k) is the Fermi distribution in the 
presence of the field, and $c=3.37A is the dis 


3 Private communication to the author. 
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tance between graphite layers. f(k) is given by 


fk) =fo( kB), (3.8) 


where fo is the undisturbed Fermi function (3.5), 
and r is the inverse of the probability of scat- 
tering per unit time. For moderate fields we may 
expand this 

f(k) =fo(k) +2(k) (3.9) 


where 
(3.10) 


er 
g(k)= ——" 
Writing 
end ent dE 
vV=- gra Ta =—- gra 9 
.~ ies ae , 
dE 
dr, =do,——__, 
|gradE | 


do, being element of length on the curve 
E=constant, we have 


adE 
“| f (F- oni, he 
|gradE| 


where the inner integral is over the surface of 
constant energy E. The component of current 
jr in the direction of the field is 


4e*r pdf i (F-gradE)? 
a 5 a 
hc J dELF |gradE | 
Averaging over all directions of the field in the 
plane of the layer we get 


2e?r df 
Gee eer S| ff raat |e. lar. (3.11) 
hc 


ee : 


By virtue of (3.1) this is 


8xe?r 


j= 
htc 


le _ 
=> og ' 
he 


To get a numerical value of the conductivity 
it is necessary to know something about 1, or 
alternatively about the mean free path /=vr. 
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We will not attempt to calculate this quantity 
in the present paper. Let us note, however, that 
experimental evidence requires that the mean 
free path decrease with temperature faster than 
T-., In this connection, one might call attention 
to the sharply increased mobility found in polar 
salts at low temperatures for thermal electrons. 
The mean free path has the form 


T\ 
1=c(—) (e'7 —1),! 
6 


where @ is the Debye temperature and C is a 
constant of the order of one atomic distance. The 
conductivity electrons in graphite are also 
essentially thermal. If the mean free path is 
assumed to be of the same order of magnitude as 
given by (3.13), with C=2X10-* cm and @ esti- 
mated from specific heat evidence as about 
2000°C we will get 


l=3X10-* cm. 


(3.13) 


(3.14) 
This gives rise to a resistivity 
p ~5X10-' ohm-cm. 


It seems therefore, that (3.14) gives at least the 
right order of magnitude for /. This point, and 
the temperature dependence of /, seem worth 
further investigation. 

From the formula 


o = Nee? t/ mers, (3.15) 


we can calculate the ‘“‘mean effective mass,”’ me¢1, 


of our conductivity electrons. This is found to be 


h*kT 
(3.16) 





Mett = 


36 log2 -a*yo? 


which yields, on numerical evaluation, at room 
temperature 


(3.17) 


4. THE BRILLOUIN ZONES OF GRAPHITE 


Mer = 1/18 electron mass. 


The graphite lattice is built up of hexagonal 
layers whose relation to each other is indicated 
in Fig. 7. Dotted lines indicate one layer, and 

*N. F. Mott and R. A. Gurney, Electronic Processes in 
Ionic Cr a (Oxford University Press, Oxford, 1940), 


p. 107. Note also the increase of mobility with decreasing 
temperature in the semiconductor Cu;O shown in Fig. 66, 


p. 168 
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the full lines the one immediately above it. A 
set of basic displacement vectors of the lattice is 


ai= [4v3a, — $a, 0], 


as=[0, a, 0], 
a3= [0, 0, c], 


c being twice the distance between layers of the 
lattice. In Fig. 8 the relation between the layers is 
shown in perspective. The points A, B, C, D are 
the four atoms of a unit cell. 

The volume of the unit cell is 


(4.1) 


Veet = &1° (@2 X @s) “— 


so that the atomic volume is 


vV3a"c 
; (4.2) 


V atomic = 


The intercepts of a possible unit cell on adjacent 
planes are shown in Fig. 8. The cell would extend 
beyond these two planes on each side halfway 
to the next planes, which are not shown. 

Reciprocal lattice vectors corresponding to 
(4.1) are 


E 2 
bi= cae ag 
Paar i V3a 


{ 1 
b: =| — 
 v3a 


' p=! 


So, & 








’ 
" 
* 


It follows tifat the lowest Brillouin zone is 
bounded by six vertical planes forming a right 
hexagonal cylgnder whose sides are distant 1/v3a 
from its axis. it appears at first that the height 
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of the lowest zone is 1/c. This is, however, not 
the case, for the plane at height 1/2c corresponds 
to Bragg reflection of waves whose wave-length 
is 2c. Thus reflections from successive planes 
(distance c/2) differ in phase by z, and therefore 
destroy each other. The structure factor there- 
fore vanishes on the planes k,=+(1/2c), and 
there is no energy discontinuity over them. The 
horizontal boundaries of the first zone are there. 
fore at k,=+(1/c). The lowest Brillouin zone 
so described has the volume 4/v3a*c. But the 
density of states is twice the volume of the 
crystal, that is, 


2X (number N of cells) 
X (volume of cell) = v3at%N. 


Thus there are 4 N electron states in the lowest 
Brillouin zone, or one per atom. Since as in the 
case of a single layer, there is just one “con- 
ductivity” electron per atom, the lowest Brillouin 
zone will be just filled at the absolute zero of 
temperature. 

In Fig. 9 we indicate the upper half of the zone 
in question. 

Let us now consider the problem from the : 
viewpoint of the ‘tight binding approximation.” 
As in the case of a single layer, let X(r) be the 
2p, wave function of an isolated C atom. Then 
we assume for the collective wave function of 
the crystal 
v= > AWM 

tmele++4 
= _ Adie exp[2mik-r*|X(r—r*)), (4.3) 


ime ]+++4 


where i=1, 2, 3,4 correspond to atoms of the 
type A, B, C, D, respectively, and a enumerates 
the different atoms of a given type in the crystal. 





of the 


V3a°CN, 


lowest 
in the 
“con. 
rillouin 
zero of 


1€ zone 


ym the » 


ation.” 
be the 
. Then 
‘ion of 


(4.3) 


of the 
erates 











Proceeding exactly as in the case of the 
hexagonal layer, we are led to a secular equation 


for the energy 
|| —E64,|| =0. (4.4) 
Now the diagonal elements are 
Hy! = Heo! = Hs’ = Hea! = Ep’ — 2y0' (cos2 rk ya 
+2 cosrv3kya cosrkya)=Ho. (4.5) 


Ey’ is not quite the same as the Zo in the 
previous section, due to a difference in the 
potential ; we shall, however, subsequently drop 
the prime for convenience. 

In calculating the non-diagonal terms, we 
shall at first consider first and second neighbors 
in the planes, and also first and second neighbors 


between planes. 
The elements may now be calculated : 


Hy. This is equal to 
—yo(exp[ — 2rik.(a/v3) ] 


+2 cosrkya-exp[wik.(a/v3) }) 
where 


vo= f X*¢—e4s)(U- V)X(r)dr>0. 


Hj. This is equal to 


2; cosrk,c 
where 


n= [ X*—eac)(V— U)X (r)dr. 


71>0 because (V — U) is negative, and 
X*(r— ac), 


X(r) have opposite signs between the planes, 
where the contribution to the integral is greatest. 
Hy. This is equal to 


2y1' cosrk,c(exp[2rik,(a/v3) ] 


+2 cosrk,a-exp[ — rik,(a/v3) }), 
where 


a f X*(r—eap)(V—U)X(#)dr, 


and is almost certainly positive, though it has 
not been evaluated. 

The other élements may be written in terms 
of the above. It may easily be verified that 


Hes= His, Ae,=H23*, Au=A:*. 
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If we put the diagonal elements equal to Ho, 
2 cosrk,c=T and 


exp[ — 2mik.(a/v3)]+2 cosrkya 
-exp[wik.(a/v3)]=S, (4.5) 

the secular equation becomes 
E-H,y -vS wv wTS 
—yoS* E-Hy wTS wTS 
mil miTS E-H, —7o5* 
mTS* ‘ ywTS -yS E-Ho 


=0. (4.6) 








The exact solution of this equation presents 
some difficulty. We shall neglect y;’ at first, and 
later treat it as a perturbation. In the lowest 
approximation, then, we get on solving, 


E=Hot hr +[hrrT? +70?! S| 2}. 


With regard to the sign in front of the $7:I term, 
the choice here is immaterial, a change of sign 
merely corresponding to reflection in the plane 
k,=1/2c of the two halves of the upper half-zone, 
and similarly in the lower half-zone. By choosing 
the negative sign, we adhere to the convention 
that the bottom of the zone (point of lowest 
energy) should be at the center. 

The signs + in front of the square root cor- 
respond to the outside and inside of the lowest 
zone, respectively, i.e. inside 


E=Ho—4}yil —([hyvrT?+707|S|*}4, (4.7a) 
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Fic. 11. Energy contours 
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and outside 
E=Ho— dy +([iv2T*+70?|S|* J! (4.7b) 


The energy inside is always <H)=Epo; it has 
the latter value for |S|?=0 and [<0 (i.e., 
k,21/2c). Now 


| S|?=(1—2 cosrk,a)* (4.8) 


on the sides of the zone, and is zero at the vertical 
edges. Thus the maximum energy inside (Zo) is 
attained on the upper (and lower) quarter of the 
vertical edges of the zone. 

By a similar argument, the energy outside has 
its minimum value, (which is also Eo) on the 
middle half of the vertical edges of the zone. The 
energy gap is therefore zero only at twelve 
points; two on each vertical edge. 

In Fig. 10 we indicate the form of the sections 


-of surfaces of constant energy by a plane such 


as that through OSVQ in Fig. 9. 

In Fig. 11 we represent the cross sections of 
the energy contours by a plane such as the ex- 
tension of PQRS beyond QS. W is at a tip of the 
six-pointed star which forms the upper boundary 
of the second zone, this zone being obtained by 
extending the vertical sides and the top and 
bottom of the first zone. 

The energy contours in planes k,=constant 
will have a form similar to those shown in Fig. 5. 
The hexagonal lines of constant energy will 
appear in each plane, but will correspond to 
different energies in different planes. The zero 
energy gap will of course only appear in the 
planes k,= +(1/2c). 

The strong anisotropy in some of the proper- 


ties of graphite (electrical and thermal op. 
ductivity, magnetic susceptibility) is caused by 
the very anisotropic form of the energy surfaces 
near the corner of the zone. These anisotropies 
will become greater as the temperature js 
lowered, and the surfaces of energy (Ey+kT) 
approach segments of lines. 

Let us now see how this picture is altered if 
we take account of +;’, the exchange integral 
between atoms such as A and D in Fig. 8. If we 
expand (4.6) to terms of the first order in +;' and 
solve, using the relation (4.7), we find that the 
energy is increased by an amount 


vo| S| 2(S+.S*) —2Eun(S?+S**) 


e=yori T , (4 
2Euyil +4y0? | S| 2 





where Eyz is the “unperturbed”’ energy given by 
(4.7). Since in the vertical edges S=0, e=0 there 
both inside and out, and consequently there is 
no overlap on these edges. Also, ¢ vanishes 
identically in the planes k,=-+1/2c. Thus, at 
the only points at which the energy is con- 
tinuous on crossing the boundary of the zone, 
e=0, and it is also identically zero in two diree- 
tions at right angles through this point. This 
alone makes it seem likely that the addition of ¢ 
will not give rise to any overlapping of the zones. 
More detailed consideration shows that such is 
in fact the case. 

It should be noted finally that the two energy 
bands corresponding to the zones considered 
above arise from the same enérgy level (cor 
responding to the 2p-state for C). The ‘‘splitting” 
of the energy level into two bands when the C 
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atoms are brought together to form the crystal 
is caused by the fact that the graphite crystallizes 
in a form having more than one atom per unit 
cell. Thus, the general properties of the zones 
(and in particular their touching) does not 
depend on the details of the potential field, but 
rather on the geometrical form of the lattice. 


5. NUMBER OF FREE ELECTRONS IN GRAPHITE 


The first problem is to find N(Z), the distribu- 
tion of energy levels. The well-known formula for 


this is 
N(E) =2V [> 
|grad,E|’ 
where V is the volume of the crystal, and the 
integral is taken over the surface on which the 
energy has the constant value EZ. The formula 
may be readily transformed to 


N(E) = (5.1) 





ak, 
. vf (0E/ Oks)oustace 


the integral now being taken over the projection 
of the surface E =constant on the (kz—k,) plane. 
Throughout this and the subsequent sections 
we shall consider only the exchange integrals yo 
and ¥; (nearest neighbors in and between planes). 
Let us write (4.7) 
e=E—Ey= —71 cosrk,c 
+ (yi? cos*rk,c +70 | S| 2)4, (5.2) 


Calculating d¢/dk, and eliminating k, we get, for 
2y1>€>0, 


~) 1 €c 
5 bese, vo?| S|?+e? 
X (4712? — (yo? | S| ?—e*)*)!. 


Now we shall be interested in «=kT<vyo. In 
this case we may write 


|S|?=3x2(k.?+k,%)a*=u. 


(5.3) 


The integral (5.1) may be transformed to an 
integral over u: 


4V (yotu+e2)du 


(e2-+2e71) /-y0? 
N(E) = f 
3ra*ce J, (y1%e? — 3 (-yo*u —e*)?)! 
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which yields, on evaluation 


N(E) = 





V 
| ant—e'4 lela 


3ra*cyo 
€ 
+2¢ sin |, (5.4) 
271 


the result has been extended to the case «<0 
by observing that it is, for e«<+yo, even in «. 

The result (5.4) holds only for ¢<2,; Coulson® 
has estimated that y:1 ~ yyv0=0.09 ev. This cor- 
responds to a temperature of nearly 1000°C, so 
that, insofar as calculations of the properties of — 
graphite near room temperature are concerned, 
the expression (5.4) will be adequate. 

It should be noticed that N(Z) does not go to 
zero when e—0, but approaches a finite value 
proportional to 7: (see Fig. 12). Nevertheless, 
as can easily be seen, at the absolute zero of tem- 
perature the number of free electrons is zero. 

The number of free electrons is in fact ap- 
proximately 

271 
Neu =2 N(E)f(E)dE, 


where 





E = ’ 
ae ver 


as in the case of the plane layer. Now with small 
error we may expand N(E£): 


8V1 
N(E) = 24a “<+0(= —)+--.] 
3ra*cy 0? n'* 


and extend the integral from 0 to ©. Doing 


NG) 








Fic. 12, 































PRIA fee) a 





: . i a SI, = = - 3 
















P. 





632 








this we obtain 


_ 32V kTy1 a kT 
nar ah 
EYoo(2Y oJ] 0 


(m+1)"" oo 
$;=log2, s2=2*/12, etc., the s,’s being tabulated 
in Dale’s Mathematical Tables. 

Now the number of free electrons per atom is 
Nett = Nets/Na, Na being the number of atoms. 
Since the atomic volume is V3a*c/8, 


N.=8V/v3a"%. 





where 





== (-1) 


| (5.7) 


#2 an ))* 


=2.25X10" 


at room temperature. This is not at all in agree- 
ment with the value of 2.310-‘* obtained by 
treating a single hexagonal layer; in fact, the 
dependence of mr; on T is quite different in the 
two cases. On the other hand, as we shall see in 
the next section, the conductivities in the plane 
agree to the first order. Thus the discrepancy is 
merely caused by our definition of mer, and is 
compensated by a corresponding change in merr, 
the effective mass. 


6. CONDUCTIVITY IN THE DIRECTION OF THE 
LAYERS AND PERPENDICULAR TO THEM 
By a procedure closely parallelling that used 
in Section 3 for the single hexagonal layer, we 
obtain for the conductivity in the direction of 
the vector u, if we assume that 7 does not depend 
spi on k, 


dfo 
“= Se 


| f (u, a)*|grad,E| 45, bz, (6.1) 
E=const 
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where n is the normal to the surface E= 
and dS, is the element of area in this surface 
Averaging over-all directions in the plane ang 
we dS; = 8rxdx/n,, where «= x2+«,2, 


dfo (dE/dx)? 
as dE/dk, 


h? dE 


ou= 


Sredx [E. (6.2) 


On the other hand, the conductivity perpen. 
dicular to the layers is 


2e*r ¢ dfo dE 
f —| f —Sride XE (6,3) 
h? dELY (x) Ok, 


On the surface ¢=constant, 


51.=- 





de ey y02(32r'a?y)! 
— i ’ (6.4) 
Ox yoy +e? 


where 3x7(k.?+h,?)a?=y. 

Now in the Eq. (5.2), if we remain within the 
limit |¢|<+y:, |S|* may be replaced by y with 
an error of order (y:/7o)?~1 percent. In terms 
of y, the integral in the square bracket of (6.2) 
may be written 











16¢ey04 ‘ ese ydy 
cy, (yet +e) [4ey:?— (ery) 
On evaluation this yields 


8e € € 
—{ #+2 sin-'— —-—__ 
c 271 [yi2—e?]} 
e+yi—[11?—-€ ]}! 
x tog 
etrmtlr’+e])! 
4y2—e }t}—[y2—e }i-—y 
~ie [4y2—e }#-—[1 J |} (6.4a) 
"Tay? —}'+[7:2-—8&}— 1 
Let us note first that if we take only the leading 
term in this expression, we get 


8xe?r dfo 
ou= = fae 














hc dE 
(6.5) 
16me?r 
kT log2 
h’c 


which is precisely what we obtained in (3.12) for 
the case of a single plane layer. It follows from 











(6.4a) 


leading 


(6.5) 


12) for 
vs from 














this and (5.7) that to the first order the effective 
mass for electrons in motion parallel to the 
layers is independent of temperature. Taking 

account of the higher terms does not change the 
result drastically for moderate temperatures; at 
room temperature (7 = 293°C) we get, using the 
numerical values already indicated (including 
the value for / assumed in Section 3), 


Oy) = Py = 5.3 X10-* ohm-cm. 


Turning now to o,, we get for the integral in 
the square bracket in (6.3) 


e 7 
=o <(: ‘ciao a -[1-] 
3 yo?a*Lyi \2 271 4y:? 
ey e+yi1—[yvi2—e*]}! 
+[1-5] (1 
vi? : e+yit+[ri?—e ]}! 
A4y,2—e? }t— 2_ 2 }i-_y 
ie ae ee )] iain 
[4y.°-—e}'+[y?-e]}!'-1 











valid up to e=~71. We shall neglect the con- 
tribution to (6.3) from beyond this point. To 
determine first the order of magnitude of the 
result, consider the leading term only under the 
bracket; this is log y:/e. Then we get 


16 e?r yic v1 ¥1 Ofo 
1. —-— — f: e? log— —-de 
3 h® oa? le| de 
16x? e?r ¢ yi(kT)? 71 
=s—— —- — ———— log—._ (6.6) 
9 ha y¥0? kT 





This expression will give a good approximation 
only at temperatures somewhat below room 
temperature. We will use it, however, to indicate 
orders of magnitude. It gives for the, anisotropy 
factor 


1. c? ykT 71 


og 
Cu 9 log2 a’ vo" kT 





(6.7) 
=0.01 at room temperature. 


Since this does not involve r it is completely 
calculable without assumptions about mean free 
path. 

A comparison of the formulas (5.5) and (6.5) 
gives 


Mest! = 2h*y1/34*a7y0" 
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for the “‘mean effective mass’’ of electrons moving 
parallel to the graphite planes. This is valid so 
long as kT<vyi, and is equal to ~}4. Similarly, 
from (5.5) and (6.6) we obtain for the mean 
effective mass for motion across the graphite 
planes 

24 log2- h? 


re*kT logy:/kT 


Mest —_ 





With decreasing temperature this increases as 
1/T|logT |, becoming infinite as 7-0. At room 
temperature its value is ~25 to 30 electron 
masses. 

More accurate calculation of ¢, by numerical 
integration at room temperature gives 


o,-'=p, =4.8 X10-* ohm-cm, 
and therefore 


o, Pu 


—=—=1.1X10-. 


ou Pr 


It is of interest to note that the anisotropy 
depends on temperature, the ratio p,/p, be- 
coming larger, i.e., the anisotropy becoming 
greater, with decreasing temperature. An attempt 
to confirm this prediction experimentally would 
be interesting; the measurement of ¢, would, 
however, necessarily be difficult. 

In Fig. 13 is plotted a graph of the ratio 
o,/o,, for temperatures up to room temperature. 


7. REMARKS ON OTHER PROPERTIES 





Several qualitative remarks can be made 
about the magnetic susceptibility and the 
volume absorption properties of graphite. Graph- 
ite is strongly diamagnetic, and the diamagnetic 
susceptibility of single crystals can be expected 
to show marked anisotropy. For a field in the 
z-direction, the susceptibility depends on the 
value of the quantity 


Oe Oe ( 0*¢ ) 
Ok? ky? \dk Ok, 
over energy surfaces in the neighborhood of the 


boundary of the Fermi distribution, that is to 
say, near the zone boundary.® This is a large 








5 See Seitz, Modern Theory of Solids, pp. 594-595. 
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quantity, of order of magnitude 9r‘a*yo*/7;?, or 

about 6 X 10‘ times its value, h*/m?, for free elec- 

trons. On the other hand, if the field is in the 

k,-direction, say, the susceptibility depends on 
0% 0% 


0 2 
ak,? ak,? —) 











Unlike the preceding case, this depends on T. 
It, is, in fact, ~3rtacekT y?/y1, or ~kT 1/70? 
tines the value for the field across the graphite 
planes. The difference may be ascribed to the 
difference in effective masses in the quantized 
orpits perpendicular to the field. The anisotropy 
is;therefore of the same order of magnitude as 
that of the conductivity. 

Our other remark concerns the “optical” ab- 
sorption spectrum. It is well known® that elec- 
trons may under the influence of radiation make 
only “vertical” transitions in the reduced-zone 
scheme, that is, transitions in which the wave 
vector changes by one of the vectors of the 
inverse lattice. Let us then consider extreme 
cases. A transition from a point such as V 
(Fig. 9) involves zero-energy change. The great- 
est energy jump is from the center of the zone 


* Reference 5, Section 71, pp. 326-8. 


3e0 «CT °K 


(O in Fig. 9) to a point at the middle of a vertical 
edge defining a “corner’’ of the star-shaped 
second zone. This is of amount 2[y:?+ 9? ? = 6y, 
Transitions involving all intermediate energies 
are possible. 

An energy jump of 6yo corresponds to a fre- 
quency of 1.31 X10'* or a wave-length of about 
2300A, fairly deep in the ultraviolet. Thus, there 
is absorption from the longest wave-lengths 
through the visible spectrum and a substantial 
part of the ultraviolet. At the extreme limits 
absorption will be weak, however, due to the 
small number of states involved. The greatest 
number of states will correspond to energies 
such as those at the mid-points of the sides of 
the zone, which will give rise to transitions with 
energy jumps of the order of 279. Thus we should 
expect the strongest absorption to be in the 
neighborhood of 6900A, i.e., in the red end of 
the spectrum. 
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Two Applications of the Variational Method to Quantum Mechanics 


WALTER KOHN 
Department of Applied Mathematics, University of Toronto, Toronto, Canada* 
(Received September 23, 1946) 












By using the unperturbed wave function with a variable scale factor as trial-function for a 
perturbed Schroedinger equation, a lower bound for the second-order perturbation energy of 
the ground state in terms of the unperturbed and first-order perturbation energy is obtained. 
A minimum property of the sum of the m lowest proper values is utilized in a new method for 
obtaining higher proper values and functions. 
















1. INTRODUCTION 





I‘ most practical problems of quantum me- 
chanics the Schroedinger equation, 


Hy=Ey, (1) 












The first-order perturbation energy, ¢:, of the 
ground state is obtained by substituting its 
unperturbed, normalized wave function, Wo, into 
the expression 


B= f vHyde: (5) 


tical cannot be solved explicitly, and one must take 
taped recourse to approximation methods. 
~ On The present note, which is based on the Thus, if we denote the unperturbed energy by 
= method of Ritz, contains (1) an estimate of the Ep», we have 

second-order perturbation energy in terms of the 
s fre. unperturbed and first-order perturbation energy; f Vo" Hyodx =Eot+kea,, (6) 
‘bout and (2) a new approach to the evaluation of the teen 
there proper values and functions of higher energy 
= states. a= f Yo*v(x)odx. (7) 
in . 
"ad . ee eee ee Now let us take as trial-function the (normal- 
atest : ay ized) function \**/*Yo(Ax), where 3n is the dimen- 
rgies We write the Hamiltonian of a perturbed gion of the configuration space and d is to be 
»s of system as adjusted to minimize EZ. Making use of (3) and 
with H=T(x)+ V(x) +ko(x), (2) (4) we obtain 
ould where 

the x represents the coordinates of the E= f Wo" (Ax) [T (x) + V(x) +hv(z) JWo(Ax)d(Ax) 

d of system; 

Tz)=—2 = & is the — energy a - repre- ° ss J 

ts) sentin Masses of the rticies ; 
V(x) oa is the cena energy anne of the = f wn(e)| r(-) + v(-) +kv (=) roca 
unperturbed system ; and (8a) 

t of k(x) is the perturbing potential energy or 
and operator. 
sity It follows from the definition of T(x) that E=»! f Yo (x) T(x) volx)dex 
7 T(x/d)=T (x), (3) 





and if we confine ourselves to electrostatic 
potentials, then 


.V(x/r) =A V(x), 





v(x/A) =do(x). (4) 





*At present in the Department of Physics, Harvard 
University, Cambridge, Massachusetts. 
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+n f Yo (x) V(x) yo(x)dx 


+yk f Yo*(x)0(x)Yo(x)dx. (8b) 


1 The integral to be taken over the configuration space. 
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The first and second integrals in (8b) represent, 
respectively, the average kinetic and potential 
energy of the unperturbed system. According to 
the Virial theorem in quantum-mechanics,’ these 
are equal to —E» and 2K, respectively. Thus 
(8b) may be rewritten as 


E= —Ep\*+(2Eo+he)a. (9) 
Minimizing with respect to \ we obtain 
A= 1+he:/2Eo, E=Eot+ke+h’'e2’, (10) 


where 
€2' =¢;"7/4Ep. (11) 


The term ke; is just the first-order perturbation 
energy. The next term, ke’, is negative and 
therefore improves our estimate of E. 

Furthermore, ¢2’ is (in absolute value) a lower 
bound of the second-order perturbation energy, 
€2. For, if we subtract the true perturbed energy, 


E,=Eotkeit+h'e2.+kh*e3+:-- (12) 
from E, as given by (10), we find 
E—E,=k(€s' —€2) —k*es—---, (13) 


and remembering now that E is an upper bound 
of E, we obtain, on letting & tend to 0, €2' 2 és. 
Since €2’ is negative this proves our statement. 
In summary we may say that if one chooses as 
trial-function for a perturbed system the unper- 
turbed wave function of the normal state, but 
with an adjustable scale-factor, the variational 
method leads to an energy which agrees with the 
perturbed energy to the first order and in addi- 
tion contains a correction in terms of the first 
order and unperturbed energy (see (11)), which 
is a lower bound of the second-order perturbation 


energy. 
Example: The Helium 1'S Term 


Hylleraas* derived the expression 


1 
E,= ~224+1.282*(—) 


1 
-031492:(—) +--+ (14) 
Z? 


2A. Sommerfeld, Wave Mechanics (Methuen and Com- 


pany, London, 1930) p. 248. é 
*E. A. Hylleraas, Zeits. f. Physik 60, 624 (1930). 


for the normal state of a helium-like atom of 
nuclear charge Ze. 

Here 1/Z corresponds to our k go 
Eo= —22?, ,=1.25Z*, and ¢.= —0.3149Z2 Our 
correction term, 


€2' = (1.25Z?)?/( —8Z?) = —0.1953Z2, (15) 
accounts for 62 percent of €:. 


3. HIGHER ENERGY STATES 


In applying the Ritz method to a higher 
state only functions orthogonal to all lowe 
proper functions may be admitted for competi 
tion. Unfortunately, the task of assuring ortho 
gonality is, in general, prohibitively complicated 

Hylleraas and Undheim‘ have used the fag 
that the mth root of a secular determinant is ay 
upper bound of the mth proper value, to approxi. 
mate the 24S term of helium. However, their 
method is in practice quite cumbersome and 
only slowly convergent. 

As an alternative we shall derive an expression 
for the mth energy state whose properties ar 
analogous to (5): Its minimum value is the ath 
proper value and the function which leads tp 
this minimum gives us the mth proper function, 

Our method is based on the following 


Lemma‘ 


Let ¥1, v2, ***Wn be the m lowest normalized 
proper functions of (1) so that 


Bim [ vtHvds, 4=1,2,---n. (16) 
Further, let 
A= ¢*Hedx, i= 1, 2, or (1) 


where ¢1, ¢2, ***¢, are any normalized ortho 
onal functions. Then 


E,+£2+---E,=min.(Hu+Heet:° Han), (18) 


where degenerate proper values are to be counted 
multiply. The minimum is attained whenever 


‘E, A. Hylleraas and R. Undheim, Zeits. f. Physik @ 
759 (1930). 

5 Courant-Hilbert, Meth. d. Math. Phys. (Verlagsbude 
oy - Julius Springer, Berlin, 1931), second editio, 
p- 399. No proof is given there but it can be easily su 

y expanding the g—s in terms of the y—s and 
deriving expressions for H;; in terms of the E;. 
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VARIATIONAL METHOD 
vy o% *** x arise out of Yi, ¥2, --*¥. by an 


Geanied transformation. 
An application to the estimation of higher 


proper values suggests itself immediately. 

Let ¢1 be a normalized approximation to the 
lowest proper function of (1) and let ¢: be an 
arbitrary function. Then the functions ¢; and 

$2 = (¢2—Tig:)/(N?—I1*)4, (19) 
with 


w= f extends, h= f tends (20) 


are normal and orthogonal. Consequently, by 
our lemma, 


Et Bx=minl f otHeds+ f o:*Hexds| (21) 


which leads to 





: Ae2—2hAwt 1A 
By=min| Hy~Ex+ | (22) 
N?-J;? 
where 
His= | oi*Heidx. (23) 


We note that (22) does not involve the correct 
¥; but only £; which can be easily and accurately 
estimated by the Ritz method. 

The ¢2’s for which the minimum is nearly 
reached are seen to be approximations to We. 

For practical purposes (22) is rather compli- 
cated. It simplifies considerably when g;=yi, 
becoming 


(E-E£,)I;? 


A 
, where E=——. (24) 
N?-—I,? N? 


E;= min| ~ 


Strictly speaking we are now obliged to use the 
exact wave function, ¥, for the computation of 
I,;. But fortunately, if ge is well chosen, the 
“un-orthogonality integral’ J; and consequently 
also the term following E in (24) are small. 
Therefore, to estimate E: great accuracy in the 
evaluation of J; is not necessary and a simple 
approximation of y¥ is quite sufficient. 

For the mth energy state one obtains, in an 
analogous notation, 


= (o- = 100s) / ue EI), (25) 
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a—l a—l 
E,=minl _ (Hu—B)+( Ha—2 > Hin 
1 i=1 


n—1 n- n—1 
+E dita) / (w-E 1)} (26) 

im1 jm i—1 
where ¢1, 92° * * ¢n—1 are any orthogonal functions 


and ¢, is arbitrary. In the special case when 
(i=, ros Pn—1 = Vn-1, we have 


n(n Et)/(v-E0) 
E e-son)/ 
("-¥u)} 0s 


Example: The Helium 2'S Term 


(27) 


E,= min| + 


To obtain an approximation to the wave 
function of the 2'S state it was assumed that 
the “inner’’ and ‘‘outer’’ electrons move in 
coulomb fields corresponding, respectively, to 
Z=2 and Z=1. This led to the function 


g2= ($r1—1)exp—[$ri] exp—[2re] 
+ ($r2—1)exp[$r2] exp—[2r,], 


where r; and rz denote distances from the nucleus. 
For the wave function of the ground state the 
approximation 


(27/16)* 


(29) 


(30) 





yi~ 


exp— [vt 


Tr 


was used. 

Substitution of (29) and (30) into (24) gave 
—4.24 Ry as approximation to EZ». This is 
still considerably higher than the observed —4.29 
Ry, but in view of the crudeness of our trial- 
function the improvement over the first-order 
perturbation value of —4.07 Ry may be con- 
sidered satisfactory. 
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UBLICATION of brief reports of important dis- 

coveries in physics may be secured by addressing them 
to this department. The closing date for this department is, 
for the issue of the Ist of the month, the 8th of the preceding 
month and for the issue of the 15th, the 23rd of the preceding 
month. No proof will be sent to the authors. The Board of 
Editors does not hold itself responsible for the opinions ex- 
pressed by the correspondents. Communications should not 
exceed 600 words in length. 





Corona from Fine Positive Points 


W. N. ENGLIsH* 
Department of Physics, University of California, Berkeley, California 
April 8, 1947 


|, = work of A. F. Kip’ has been extended to needle 
points of 0.002 mm radius of curvature. Three 
oscilloscopes were used to minimize the effect of amplifier 
pass band—a Navy type TS 28/UPN synchroscope and 
amplifier AM 73/UPA-1 (1000 cycles to 5 megacycles), a 
Dumont type 241 (20 cycles to 2 megacycles), and a 
Dumont type 171 (15 cycles to about 100 kilocycles) as 
used by Kip. Kip’s observation of the degeneration’ of 
streamer pulses with decreasing point size, until indis- 
tinguishable from burst pulses, was found to be instru- 
mental. Such pulses were clearly resolved on the synchro- 
scope for the finest points, and very similar in duration 
and shape to streamer pulses for large points, when reduced 
to the same amplitude. 

Dust-free points finer than about 0.01 mm radius of 
curvature gave no observable bursts or streamers. Corona 


int streamers. Right to left, point 


Fic. 1. Phetapaghe of pa 
radii 0.003, 0.003, 0.04, 0.1, 0.25 mm. Exposures 


of curvature: —0.00 
5 to 10 minutes at f:2. 
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Fic. 2. Synchroscope pulses at 6000 X amplification Corresponding 
photographs in Fig. 1. Sweep from right to left. ® 


onset was then evident on the oscilloscopes only by, 
slight and voltage-critical irregularity of the trace, Fig 
positive points dusted with MgO gave bursts and/e 
streamers at the onset voltage of the intermittent regine 
and above until steady corona occurred. The amplitude 
and frequency of the streamers varied somewhat from on 
dusting to another, and usually decreased with time, pr 
sumably because dust specks were removed from th 
critical region. This dust effect has been observed by 
Hudson? for smooth negative points as large as 0.25 mm 
radius. 

Photographs of the discharge at the point were taken to 
establish a connection between the synchroscope pulses 
and “visible” (i.e., photographable) streamers. Some of 
these are shown in Fig. 1 while the corresponding synchro 
scope pulses are drawn to scale in Fig. 2. It is noteworthy 
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Fic. 3. Streamer length on photograph plotted against amplitude d 
onoscope pulse at 6000 < none Observed variation ia 
pulse amplitude is indicated. 
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that all the streamer pulses are very closely of the same 
duration, 0.4 microsecond. A plot of pulse amplitude 
against “visible” streamer length, Fig. 3, shows that the 
charge per cm of streamer is essentially constant and 
independent of streamer length. Using an estimatéd pulse 
input impedance of 50 ohms, the slope of the curve gives 
a value of the charge per cm length of “visible” streamer 
of 10" ions. This agrees well enough with the value found 
by Kip of 5X 10° ions per cm. 

The radius of curvature of the steel needle points was 
measured from micro-photographs with up to 4500 times 
magnification. Values were checked among the points by 
plotting the onset voltages for intermittent corona. The 
experiments were done in room air at atmospheric pressure, 
with a 6-cm gap, and with 50 grams of thorianite ore at 
12 cm from the point. 

The results obtained indicate that streamers and bursts 
occur, and are clearly distinguishable for points down to 
0.002 mm radius of curvature, provided dust specks are 
present on the point. With larger points the streamer 
pulses have much greater amplitude than the bursts (a 
factor of about 10, for 0.25 mm radius). Both streamer 
and burst pulses decrease in amplitude with decreasing 
point size until they are barely detected, but the streamer 
pulse maintains its form and 0.4 microsecond duration, 
while the burst pulse, which is a little bit of stable corona, 
can have any duration, usually 100 to 300 microseconds 
or greater. 

The role of the MgO dust is not entirely clear. It may 
well be to provide triggering electrons, by the Paetow 
effect,? in the very small high field region. On the other 
hand, the variation of pulse size with different dustings 
suggests that the larger particles, which have a size 
comparable to the point size, may act to provide a less 
divergent field. This should cause a change in the inter- 
mittent regime onset, which has not been observed. 

The writer is indebted to Professor L. B. Loeb for 
advice and encouragement. 

* Eastman Kodak Fellow. 
14. F. Kip, Phys. Rev. 55, 549 (1939). 
F. Kip, 


?L. B. Loeb, A. G. G. Hudson, and W. H. Bennett, Phys. 


Rev. 60, 714 (1941). 
+L. B. Loeb, Fundamental Processes of Electrical Discharge in Gases 
(John Wiley and Sons, Inc., New York, 1939), p. 498 





Inversion Spectrum of Ammonia 


M. W. P. STRANDBERG, R. Kyui, T. WENTINK, JR., AND R. E. HILLGER 
Research Laboratory of Electronics, Massachusetts Institute of Technology, 
Cambridge, Massachusetts 
April 5, 1947 


OMPARISON of our’ results on the frequencies of 
the lines in the N"*H; inversion spectrum, with those 
of Good and Coles,? show a discrepancy of 4.56 Mc/sec. in 
the J=4, K =3 line, which is well outside the experimental 
error. A study of our results has disclosed that a computa- 
tional error was made in the reduction of the data, and our 
best value is 22,688.18 Mc/sec., which is in fair agreement 
with Good and Coles’ value. This would change the value 
of the anomalous frequency shift for this line in Table II 
of reference 1, to +1.76 Mc/sec. 
It is interesting to note that if the data for the N"H; 
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lines, given by Good and Coles, for the K =3 transitions, 
are analysed in a manner similar to that in reference 1 for 
the NH; lines, corresponding anomalous frequency shifts 
for the two isotopic molecules are obtained. — 








¥ measured Ypredicted 
J K NH; N“H, 
3 3 —0.32 Mc/sec. —0.30 Mc/sec. 
4 3 +1.76 +1.83 
5 3 — 6.96 —8. 
1 Strandberg, Kyhl, Wentink, and Hillger, Phys. Rev. 71, 326 (1947). 
2? W. E. Good and D. K. Coles, Phys. Rev. 71, 383 (1947). 





Nuclear Quadrupole Moment Effects in the 
Microwave Spectrum of Ammonia 


RICHARD J. WaTTs AND DuDLEY WILLIAMS 
The Ohio State University, Columbus, Ohio 
March 27, 1947 


HE hyperfine structure first observed in the inversion 

spectrum of ammonia by Good! has been the subject 
of much recent work. Coles and Good? have reported 
measuring the satellite separations for many of the lines 
and have published a formula which gives the: expected 
separations of the satellite lines from the main lines. 
Professor Van Vleck* has also developed a theory of the 
hyperfine structure. According to these theories, the 
satellite lines can be explained in terms of the interaction 
between the electrical quadrupole of the N™ nucleus and 
the electrical field of the NH; molecule. , 

The intensity of the satellite lines decreases rapidly as 
the rotational quantum number J increases. Dailey et al.* 
have reported measurements of the satellites of lines with 
J values as high as 4 and find these to be in excellent 
agreement with Van Vleck’s theory; but they have reported 
their inability to detect satellite lines for J values higher 
than 4. 

Recently, we have made a study of the ammonia 
spectrum and have been able to observe and measure the 
relative separations of satellites from main lines with J as 
high as 7. The results for lines 11, 22, 33, and 44 are in 
agreement with Dailey’s results and the relative separa- 
tions of all lines observed agree with Van Vleck’s theory. 
Our preliminary results are summarized in Table I. A 
more detailed account of this work will be published in 
the near future.* 


Taste I. Observed and calculated ratios of satellite displacements. 











} Line Observed ratio Calculated ratio 
JK Av’ /Ar * Av’ /Ar 
il 2.55 2.50 
22 1.56 1.56 
33 1.36 1.35 
1.28 1.26 
55 1.21 1.20 
66 1.18 1,17 


1.19 1.14 





. Phys. Rev. 70, 979 (1946), 
Van Vieck, and Wilson, Phys. Rev. 70, 


in connection with a research project 
by the Watson Laboratories of the Air Material Command, 
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Microwave Spectra: The Hyperfine Structure 
of Ammonia* 


WALTER GorRDY AND MYER KESSLER 
Depariment of Physics, Duke University, Durham, North Carolina 
April 2, 1947 


SING a 12-meter wave-guide cell, a sensitive detecting 
system, and a pressure of approximately 10-' mm Hg, 

we have re-examined! N“H; spectra in the region of 1.25 
cm, with the result that the small satellites of certain of 
the lines first detected by Good! are strikingly resolved. 
The upper curve in the accompanying photograph (Fig. 1) 
is a cathode-ray scope display of the 3,3 line*® with its 
satellite pairs on opposite sides. These satellites are be- 
lieved* to be a splitting of the center line by the nuclear 
quadrupole moment of N™. The wide separation of the 
satellite pair which are only about 0.6 Mc apart, illustrates 
the remarkable resolution now obtainable with electronic 
techniques in the microwave region. We have also made 
photographs showing the similar splitting in the 1,1; 2,2; 
4,4; 5,5; 6,6; and 8,8 lines.* The lower curve in the figure 
shows that obtained for the 6,6 line. The satellites are 
closest to the center line for the 1,1 case and move out 
from the center line as the quantum numbers of angular 
momenta increase. Accurate measurements of their posi- 
tions‘ are being made and will be reported later. The 
intensity of the satellites, as compared with that of the 
center line, decreases with increasing quantum numbers of 
angular momenta. We have not yet photographed the 
fine structure of the 7,7 and 9,9 lines because our oscillator 
would not operate in their regions. We observed several 





Fic. 1. Upper curve: 3,3 line of ammonia with satellite structure. 
lout curve: 6,6 line, with satellite structure. =10-§ mm 
g- 
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strong lines for which J and K are unequal but d 
detect any satellite structure.’ 

In this work we did not observe the broadening Caused 
by saturation which was reported by Townes.! The 
lines shown in the photographs were observed whe 
approximately 10 milliwatts of power were flowing through 
the cell. 

These results were obtained by a single crystal, yp. 
balanced detecting system which we have evolved, In 
brief, we use a slowly changing reflector voltage to modu. 
late the source and a low frequency filter in the amplifier 
so that the indicator does not respond to the gradual 
change of voltage over the envelope of the oscillator mode 
but is very sensitive to the abrupt changes in output 
caused by the sharp absorption lines. This system has 
the additional advantage that spurious signals caused by 
reflections are largely eliminated, since they are usually 
much broader than the absorption lines of a gas at low 
pressure. Our tests indicate that this system is superior 
to the balanced detectors which employ either two crystals 
and cancellation in the video stage or a single crystal and 
cancellation in the r-f stage. 


id not 


* This research was supported by a Frederick Gardner 
Special Grant from the Research Corporation and by a loan Pi 
ment from the Office of Research and Invention, United States Navy, 

1C, E, Cleeton and N. H. Williams, Phys. Rev. 45, 234 (1934); B, 
Bleamy and Penrose, Nature 157, 339 (1946); W. E. Good, Phys. Rey, 
70, 213 (1946); C. H. Townes, Phys. Rev. 70, 665 (1946); M. W. P. 
Strandberg, R. Kyhl, T. Wentink, Jr., and R. E. Hillger, Phys. Rey. 
71, 326 (1947). . 

2 Strandberg, Kyhl, Wentink, and Hillger (reference 1) and W. E. 
Good and D. K. Coles (Phys. Rev. 71, 383 (1947)) have provided 
accurately measured frequencies of the centers of the first six of these 


lines. 

3B. P. Dailey, R. L. Kyhl, M. W. P. Strandberg, J. H. Van Vleck, 
and E. B. Wilson, Jr., Phys. Rev. 70, 984 (1946); D. K. Coles and 
W. E. Good, Phys. Rev. 70, 979 (1946). 

4 Dailey, Kyhl, Strandberg, Van Vleck, and Wilson (reference 3) 
pave measured the satellite separations for the 1,1; 2,2; 3,3; and 44 
ines. . 





Bond Distances in OCS from Microwave 
Absorption Lines 


T. W. Dakin, W. E. Goop, AND D. K. CoLes 
Westinghouse Research Laboratories, East Pittsburgh, Pennsyloania 
April 8, 1947 


HE rotational absorption line for the J=1 to J=2 

rotational transition of the OCS linear molecule was 
reported by Dakin, Good, and Coles.' Absorption at higher 
pressures caused by the transition had been reported 
previously by Hershberger.? The absorption line reported 
above was caused by the molecule composed of the 
abundant isotopes: carbon 12, oxygen 16, and sulfur 32. 
Of these three elements, only sulfur has another isotope 
with an abundance of more than 1 percent of the total 
naturally occurring element. The sulfur isotope 34 has an 
abundance of about 3 percent. It was believed that the 
sensitivity of our microwave absorption apparatus de- 
scribed by Good,’ with a wave-guide absorption cell 36 
feet long, was just about enough to detect the absorption 
line caused by the OCS molecule with the sulfur 34 isotope, 
and after a preliminary calculation of the approximate 
frequency of this line, using the bond distances reported 
by Cross and Brockway,‘ it was searched for and eventually 
found. The resolution of this isotopic absorption line was 
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TABLE I. Data on bond distances in OCS. 











— Effect of uncertainty on 
Constant Value CO distance CS distance 
— 1.9823 +.0.00030 +0.00002 +0,00002 
i  semnaaieuns aie” oat” 
Atomic tC 12.003861 +0. ' ¥6; 
Atomic weight O 16,0000 
Atomic weight 1.9957 +0.002 +0.0047 + 0.0038 
ose 
- of  (1.66035-40.00031) +0,00011 +0,00015 
Mase af atom ( pear at 
naure Constant (6.6242 +0.0024) +0.00022 +40.00030 


X10-?? erg sec. 


Meise of inertia (137. 9540+0.0001)  £0.000009 + 0.000005 
+ be = x 
(* ) F000. 0.00075 +.0.00062 


(Se molecule). (3.4564 0:0003) 
Mesent X10-# g cm? 


+0.0058 +0.0049 











first reported by us at the National Research Council 
Conference on Insulation in November, 1946. 

The frequencies of these two absorption lines have since 
been measured more accurately by comparison with a 
harmonic of an oscillator which was heterodyned to zero 
beat with the 5-megacycle WWV National Bureau of 
Standards station as used by Good and Coles.*§ The fre- 
quencies of the absorption lines can thus be determined to 
about one part per million. They are f(S® line) = 24,325.92 
+0.02 Mc/sec., f(S* line) =23,731.332%0.03 Mc/sec. Un- 
fortunately, because the nature of the calculation and the 
uncertainty of other constants necessary, all of this accu- 
racy is not available in the two bond distances which may 
be calculated from these two frequencies. 

In Table I there is a list of all the constants used in the 
calculation of the bond distances from the frequencies. 
The values of these constants, together with their uncer- 
tainties are given, and the effect of these uncertainties on 
the calculated bond distances is given. It will be observed 
that in most cases the effect of the error is in the opposite 
direction on the two bond distances so that their sum can 
be calculated more accurately than either separately. It 
was hoped that this very accurate frequency measurement 
would give a much more accurate value of the bond 
distance. An examination of Table I shows that the 
principal uncertainty is the difference in atomic weight of 
the sulfur isotopes. This contributes by far the largest 
amount to the uncertainty of the bond distance. The 
possible error in the frequency difference between the 
absorption line of the two isotopic molecules is the next 
largest uncertainty. 

The bond distances C—O, (1.1612+0.0058) x 10-* cm, 
and C—S, (1.5604+0.0049) x 10-* cm calculated here are 
the average interatomic distances in the lowest vibrational 
state. The atoms are assumed to be point masses and the 
centrifugal stretching has been approximated from the 
expressions given by Nielsen® and the vibrational fre- 
quencies of Penney and Sutherland’ and found to be less 
than one part per million for the rotational levels involved 
here. The uncertainty in these bond distances can be 
reduced to about 0.1 percent when the difference in mass 
of the two sulfur isotopes is known to be better than 
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+0.0001 atomic weight units. The isotope atomic weight 
values used were obtained from the table of Evans.* The 
values of interatomic distances given here compare very 
well with the values obtained by Cross and Brockway 
from electron diffraction measurements: C—O = 1.16+0.02 
X10-* cm, C—S=1.56+0.03 x 10-* cm 


W. Dakin, W. E. Good, and D. K. Coles, Phys. Rev. 70, 560 


a) (ives). ershberger, R.C.A. Review 7, 422 (1946; J. App. Ph 
ev 8. 
17, ait (1946), an 


hem. Phys. 3, 821 (1935). 
, 383 (1947), 





. ). 
B. B. M. Sutherland, Proc. Roy. Soc. 156, 
R. D. Evans, Introduction to the Atomic Nucleus (M.1.T. 1944). 
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High Energy Ions in Crystal Lattices 
A. H. W. ATEN, Jr. 
Natuurkundig Laboratorium der N. V. Philips Gloeilampenfabricken, 
Eindhoven, Nederland 
April 7, 1947 


ANY complex anions are known to give rise to a 
Szilard-Chalmers effect when irradiated with slow 
neutrons.' In the case of a perchlorate, after neutron 
capture the activity is found entirely in the chloride 
fraction. This means more than the simple fact that the 
recoil separates the central chlorine from its surrounding 
oxygens. It is customary to consider the ClO,- ion as 
having a composition intermediate between the formulae 
Cl-O, and Cl"*(O*~),, but either all the activated chlorine 
is ejected as negative ions or, if part is ejected with a 
positive charge or as a neutral atom, it takes up electrons 
before it can react with any of the oxygen in the lattice. 
In other cases (e.g., arsenates and arsenites) the radio- 
activity is found entirely in ions of the same formula as 
those of which the crystal consists. Thus the capture of a 
neutron by a nucleus causes either a chemical reduction 
or it leaves the valency.of the activated ion unchanged. 
Comparison between different reactions of this type is 
complicated by the fact that both the chemical conditions 
of atoms ejected by recoil and their chances of subsequent 
reactions depend on the crystal lattice. 

In a note giving some preliminary figures,? it has been 
pointed out that the position is quite different in case the 
radioactive nucleus originates from an (n,p) or an (n,a) 
reaction. If fast neutrons of normal energy are used in 
such a process the kinetic energy of the recoil of the 
activated nucleus is very roughly 100 to 1000 times as 
large as the energy of a nucleus which has captured a 
neutron. Now this energy is not only sufficient to break 
the chemical bonds which connect the activated atom 
with the surrounding oxygens but even to strip it of many 
of its electrons during a series of collisions with neighboring 
particles. The highest positive charge which the activated 
atom attains will depend on the energy it has been given 


- by the nuclear reaction and will be influenced very little 


by the composition of the lattice. In general, this highly 
charged ion will be very unstable chemically, and after 




















losing its excess energy it will settle down in the Jowest 
energy state it can reach in the crystal. Which compound 
will be formed during this process will depend on the 
chemical composition of the lattice. This makes it possible 
to compare the fate of the same activated monatomic ion 
in different crystals. In general, we may expect to find the 
products of (,p) and of (,a) reactions in a highly oxidized 
state. Apart from this, a large concentration of oxygen in 
the lattice will favor the formation of oxygen-rich anions. 
On the whole, the results of the measurements in Table I 
are in agreement with these conclusions. 


TABLE I. Valency of radiophosphorus from sulfur and chlorine with 
fast neutrons. Fraction of total radiophosphorus (in phosphate +-phos- 
phite +hypophosphite) found in phosphite +hypophosphite in different 








F ~— of P® found 
phosphite 
Crystals irradiated +hypophosphite 





KsS0, 0.09; 0.10; 0.10 
KsSO,, dissolved in sulfite solution 0.09; 0.10 
NasSO,°10 aq. 0.02; 0.02 

NasSOs3°7 aq. 0.50; 0.56 

NasS hydrate 0.64; 0.78; 0.93; 0.93 
Ki ry 0.01; 0.01 

KC1O; —0.01; 0.01 

MgCl: hydrate 0.07 

MgCl: hydrate dissolved in sulfite solution 0.02 

NaCl 0.62; 0.68; 0.69 








A radium-beryllium neutron source was used, and the 
irradiated products were dissolved in water containing 
phosphate, phosphite, and hypophosphite ions, and, in 
the case of potassium chlorate, ammonia. Separations were 
carried out as described by Libby, hypophosphite following 
the phosphite fraction. Although it remains to be investi- 
gated in how far the results depend on the time elapsed 
between the irradiation and the separation and on the 
temperature of the experiment, the good agreement be- 
tween duplicates is sufficient proof of the reliability of the 
figures for comparative purpose (the value —0.01 found 
for KCIO; was included in the table to give an idea of the 
degree of experimental accuracy). It is at present an open 
question whether the phosphate and phosphite ions are 
formed inside the lattice or whether they adopt their 
final state at the moment the crystal is dissolved. In this 
connection it is of interest that the presence of reducing 
sulfite ions does not influence the result. 

In the NaCl crystal the primary product will consist of 
chlorine-phosphorus complexes, which hydrolyze when 
dissolved in water. Evidently the presence of oxygen in 
the MgCl; hydrate crystal favors the phosphate formation. 
In the NaS hydrate crystal many different types of 
thiophosphates will be formed which do not hydrolyze as 
easily as the chloro compounds. Their fate during the 
chemical separations will depend on absorption processes 
and will be quite irreproducible, as is seen clearly in the 
table. These figures for NasS hydrate are, of course, 
without meaning. 

1W. F. Libby, J. Am. Chem. Soc. 62, 1930 (1940). 
2A. H. W. Aten, Jr., Rec. Trav. Chim. Pays-Bas. 61, 467 (1942): 
considerations concerning - unds have been 


era com 
pu by Yankwitch, Rollefson, an orwick, Pe Chem. Phys. 14, 
31 (1946). 
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A Note Concerning the Motion of Arc Cathode 
Spots in Magnetic Fields* 


R. L. Loncini 
Westinghouse Research Laboratories, Pittsburgh, Pennsyleanig 
February 17, 1947 


A NUMBER of observers'* have reported that log 
pressure arc cathode spots, when in a magnetic 
move in a direction opposite from what one might expert 
by an application of Ampere’s law to the interaction of the 
magnetic field and arc current. The following explanatiog 
is consistent with the field emission theory of cq 
spots‘ which requires a strong electric field at the cathode 
This field is due to an excess of ions over electrons in the 
space between the plasma and the cathode. Electrons 
most of the current but, because of their relatively great 
velocity, have low density in this transition space an 
therefore make only a secondary contribution to th 
electric field at the cathode. As the field emission increaggs 
very sharply with field strength, the cathode-spot position 
will be that of the peak field. It will be shown that th 
electron space charge has a significant effect on the position 
of the field’s peak, and therefore influences the cathode 
spot position. 

Quantitative calculations of the effect of magnetic fields 
on cathode electric field intensity would require a thorough 
knowledge of space charge distribution. Since this js 
unknown, only a qualitative discussion will be given here 
For simplicity of description, we shall assume independent 
electric fields at the cathode, one caused by ion space 
charge (ion field) and one caused by electron space charge 
(electron field) (ignoring interactions between these). Figure 
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Fic. 1. Gain electric fields at the cathode surface. (M 
field into apex.) Figure A—Zero magnetic field. Figure 
magnetic field. Shows retrograde motion of the peak field or cathode 
aa. Mowe C—Large magnetic field. Shows the spot moving in th 

ion of the elementary particles. 







1 shows the two fields and their sum, the net field, th 
peak of which determines the cathode-spot position. 
Ions and electrons, having both opposite charges and 
directions of motion, will be deflected in the same direction 
by the magnetic field with the electrons being deflected 
hundreds of times farther than the ions because of their 
small mass. Thus when a small tangential magnetic field 
is applied, the average instantaneous position of the ions 
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and electrons (the space charge) will be displaced in the 
direction of their deflection. The electron field at the 
cathode will no longer partially neutralize the ion field 
symmetrically as in Fig. 1A but will have a greater neutral- 
izing effect in the direction of their displacement, be- 
cause of the larger electron displacement, as shown in Fig. 
1B. Thus, surprisingly enough, while each of the com- 
ponent fields is shifted to the right, the peak of the net field 
is shifted to the left, thus accounting for the observed 
motion of the spot. 

If the magnetic field intensity is increased greatly in 
magnitude, the electron and ion fields will be separated 
to the point of resolution, and then further displacement 
of the electron space charge would not influence the 
position of the net field maximum which would then be 
solely a function of the ion space charge. The displacement 
of the latter, small though it is, now solely determines the 
motion of the cathode spot (see Fig. 1C). This theory 
thus demands that the spot motion should reverse its 
direction as the magnetic field is greatly increased. This 
reversal has, indeed, been observed.* : 

Random motion will be superimposed on those described 
above because of irregularities in the cathode shape which 
will cause variations in the net field. 


* Taken in part from Master's thesis submitted to the University 


tsburgh. 
te N. Minorsky, J. de phys. et rad. 9, 127 (1928). 
?C. G. Smith, Phys. Rev. 62, 48 (1942); 64, 40 (1943); 69, 96 (1946). 
°C. 3. Lan ay and J. D. Cobine, Bull. Am. Phys. Soc. 22, 28 (1947). 
‘ , Gen. Elec. Rev. 26, 731 (1923); Zeits. f. Physik 46, 
( 





Activities Induced by Pile Neutron Bombard- 
ment of Samarium* 


Mark G. INGHRAM, RICHARD J. HAYDEN, AND Davin C, Hgss, Jr. 
Argonne National Laboratory, Chicago, Illinois 
March 27, 1947 


SAMPLE of Sm.0; was given a long neutron bom- 

bardment in the Hanford pile. By mass spectro- 
graphic analysis' active masses were observed at mass 
numbers 145, 151, 152, 154, 155, 156, 161, 167, 169, and 
176. The interpretation of these active masses is given in 
Table I. Parentheses are’ placed about those half-lives 
which are tentative. 

The assignment of samarium to certain of these active 
masses was made because samarium emits ions of form 
Sm* and SmO* from the heated oxide. The assignment 
of europium to other active masses was made because 
europium emits only as the metal. The 169 activity could 
conceivably have been caused by samarium, gadolinium, 
erbium, or ytterbium. Since no line was observed as mass 
153 the possibility of samarium was ruled out. Gadolinium, 
which emits mostly as GdO*, could have caused the line 
as could erbium or ytterbium emitting as Er*+ and Ybt. 
However, since gadolinium was a known impurity in the 
samarium sample and no erbium or ytterbium impurities 
could be detected, it was concluded that the 169 mass line 
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Mass No. Ion Active isotope Half-life 
145 Sm*#+ Sm >72 days 
151 Sm+ Sm! ~20 a 
152 Eu+ Eu ~ 5 years 
154 Eul+ Eu ~ 5 years) 
155 Eulss+ Eu —3 years 
156 Eulse+ Euls¢ 15.4 days 
161 Sm“4#O+ Sm“ >72 days 
167 Sm!#O+ Sm! (~20 years) 
169 Gd'8O+* Gda's >72 days 
176 Tbi@o+ Toe 72 days 








was caused by a gadolinium isotope of mass 153. The 
mass observed at 176 was caused by 72-day terbium? 
emitting at TbO*. 

The Sm™* was formed by (n,7) reaction on Sm™, It 
probably decays by X-capture or positron emission to 
615, Since the ratio of blackening at 145 and 161 positions 
is that characteristic of samarium, the half-life of the 61"* 
cannot be of the same order of magnitude as the half-life 
of Sm"*, The half-lives of Sm™* and Gd'™ were shown to 
be greater than 72 days by comparison with the decay of 
72-day Tb'®. The Sm‘ has previously been observed 
only in fission.* It decays with a half-life of roughly 20 
years to Eu'*!, The europium activities at masses 152 and 
154 were formed by (n,7) reactions on impurities of 
europium in the samarium sample. In addition an appreci- 
able part of the 154 observed was formed by the reaction 


(n,y) 8 (n,y) 
>Sm!8 —Eu!&——>Eu!, 
T; =47 hr. 
The 2-3 year Eu'** was formed by the reaction 


(n,y) B 
Sm}5——_Sm155 


Sm! 





Eu, 





T,=21 min. 


It has previously been observed only in fission.2 The 
15.4-day Eu'* had also been observed only in fission. 
Here it was formed by (n,7) reaction on the Eu"*, The 
cross section for this reaction was calculated to be roughly 
14,000 x 10-** cm*. This is the first example of the-deter- 
mination of an absorbing cross section for a radioactive 
isotope by means of the mass spectrometer. The half-life 
of the europium activity at the 156 position on the plate 
was shown to be 15.4 days by the following method. 
The active plate was placed successively against various 
parts of a larger photographic plate for times calculated 
to give equal intensity at the 156 position if the half-life 
were 15.4 days. Four exposures were taken for successive 
times of 5 days, 6.4 days, 9.2 days, and 15.8 days. The 
lines at masses 152, 154, and 155 became denser with 
successive exposure, but the intensity of the 156 lines 
remained constant thus verifying that the half-life of the 
Eu'** was actually 15.4 days. The half-lives of the Sm™*, 
Gd, and Tb'® were determined by a similar method. 

* This document is based on work performed under Contract No. 
W-31-109-eng-38 for the Manhattan Project at the Argonne National 
Rt Hayden and M. G. Inghram, Phys. Rev. 70, 89 (1946). 


? Bothe, Naturwiss. 31, 551 (1943). 
+ Rev. Mod. Phys. 18, 513 (1946). 
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The Quadrupole Moments and Spins of BR, 
Cl, and N Nuclei 
C. H. Townes, A. N. HoLpen, J. BARDEEN, AND F. R. MERRITT 


Bell Telephone Laboratories, Murray Hill, New Jersey 
April 4, 1947 


REQUENCIES and relative intensities of absorption 
lines near 24,000 megacycles in gaseous BrCN and 
CICN have already been reported These lines represent 
rotational transitions J=2 to J=3 for BrCN and J=1 
to J=2 for CICN, which are split into a hyperfine struc- 
ture. It has also been pointed out that the main features 
of this splitting appear to be caused by nuclear quadrupole 
moments.? 

More detailed analysis of the spectra yields information 
about nuclear moments of Br, Cl, and N are shown in 
Table I. 

The spins of the bromine nuclei are in agreement with 
values given by Tolansky,‘ but spins for the chlorine nuclei 
found by this work do not agree with the previously 
determined® value of 5/2. The spin of N™ is, of course, 
well known, and its quadrupole moment has been detected 
in the ammonia spectrum.® In the present case the coupling 
3.66 Mc is not far different from the value 4.1 Mc for the 
ammonia molecule, and here we may determine the 
quadrupole sign, which could not be obtained previously. 
To determine the quadrupole sign for each nucleus, the 
sign of &*V/d,? is assumed to be the same as that found in 
the hydrogen molecule.” 

About one-third of the BrCN molecules are excited to 
the lowest vibrational (bending) mode at room tempera- 
ture. Eight weak lines appear to be caused by these 
excited molecules, and in order to fit the observed positions 
of these eight lines, values must be assumed for the change 
in rotational constant with vibrational excitation and for 
l-type doubling. These values are consistent with what is 
known about such effects but cannot be exactly calculated. 
The rotational constant By for the ground state is taken 
smaller by Bp/340 than B,,,.1, the rotational constant for 
molecules excited to the first level of the bending vibra- 
tional mode. The /-type doubling, which is present in this 
excited state, gives a difference in the apparent value of 
Bv,<1) for the two states of By»,.1)/1000. Intensity ratios 
between two of the BrCN lines attributed to excited states 
and two lines attributed to the ground state were observed 
at +65°C and —40°C and found to change with tempera- 
ture as expected by approximately a factor of two. 

Agreement between experimental results previously 


TABLE I. Data on nuclear quadrupole moments. 
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3.66 +0.15 








he coup) is defined as (eQ)3*V /8.2 where Q is the nuclear 
q pole moment,’ is the electric potential produced at the 
nucleus by A rest of the molecule, and s is the coordinate along the 
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Fic. 2. Comparison of Cl*CN spectrum observed under low resolution 
with calculated pattern caused by Cl quadrupole moment. 
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published and theoretical expectations® for BrCN is shown 
in Fig. 1, and for CESCN in Fig. 2, assuming quadrupole 
splitting caused by Br and Cl nuclei only. One of the lines 
shown for CICN was not originally reported, since it had 
not at that time been resolved. 

In order to detect the smaller splitting caused by 
nitrogen which was superimposed on the spectra of Figs. 1 
and 2, it was necessary to increase the resolving power of 
the apparatus previously used. The oscilloscope trace of 
the two central lines of CICN shown in Fig. 2 is displayed 
under high resolution*® at two different pressures in Fig. 3. 
As a comparison, the theoretical splitting expected from 
N* is included in the figure. The weaker line of the central 
pair is not split at all by N™ effects and is the lowest 
frequency component of this group. Other CICN lines 
were observed under high resolution, and their splitting 
likewise agreed with expectations. The exact and detailed 
agreement between observations and theory leaves little 
doubt that the true value of the spins for the chlorine 
nuclei is 3/2 rather than 5/2. 

Figure 3 gives an interesting indication of the high 
resolving power which can be obtained: Here, in spite of 
considerable difference in intensity, lines are resolved 
which are separated by only 140 kilocycles or 5 X 10-* cm='. 


1 Townes, Holden, and Merritt, Phys. Rev. 71, 64 (1947). 
* Townes, Holden, and Merritt, New York Meeting Am. Phys. Soc. 


X38. 
Pape cliogg, Rabi, Ramsey, and Zacharias, Phys. Rev. 57, 677 (1940). 
4 Tolansky, Proc. Roy. Soc. 136, 585 (1936). 
* Elliott, Proc. Roy. Soc. 127, 638 (1930); Shrader, Phys. Rev. 64, 


1943). 
5 ey, Kyhl, Strandberg, Van Vleck, and Wilson, Phys. Rev. 70, 
984 (1946); Coles and Good, Phys. Rev. 70, 979 (1946). 
7 Nordsieck, Phys. Rev. 58, 310 (1940). 
8 Details of theory and of experimental methods for obtaining high 
resolution will be published later. 





Proposed Methods of Detecting the Neutrino 


K. C. WANG 
Department of Physics, National Chekiang University, Hangchow, China 
March 24, 1947 


T has been pointed out by. the writer! and others that 
the experimental investigation of K-electron capture 
of Be’ may give a relatively easy way of detecting the 
neutrino. Allen? has actually carried out such an experi- 
ment with a positive result. However, this cannot be 
regarded as a decisive proof of the existence of the neutrino 
because, as commented on by some authors,’ the method is 
still subject to relatively large corrections. We here propose 
some new methods of attacking the problem. 

Recent experiments of Baldwin and Klaiber‘ have shown 
that, with x-ray of nearly 100 Mev, nuclear reactions of 
the types Za(y,")Za-1 and Za(y,pn)(Z—1)a-2 can occur 
with rather large probabilities. It thus seems not too 
difficult to prepare the nucleus B® from C" or N™, If a 
cloud chamber is filled with gases containing C¥ or N“ 
and then exposed to a source of 100 Mev x-rays, B” will 
be produced in the chamber. Since B” has a very large 
beta-ray energy (12 Mev) and its mass is comparatively 
small, it may be much easier to detect and to measure 
the momentum of the recoil nucleus than in the case of 
Ci* (beta-ray energy 5 Mev) investigated by Crane and 
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Halpern.® If, because of its too short lifetime (0.022 sec.), 
the investigation of the B™ recoils as suggested above is 
found not to be practicable, we propose the use of N™ 
from the reaction O'*(n,p)N"*. Let the cloud chamber be 
filled with oxygen; under the exposure to neutrons from 
a Li+D source, some N"* nuclei may be formed in the 
chamber. The expansion of the chamber then follows 
immediately after the cutting off of the neutron source 
because the lifetime of N"* is only 8 sec. So far as the 
observation of the recoil nuclei is concerned, N" is also 
more favorable than Cl** although not so favorable as 
BY, . 
There is another method which may also prove valuable. 
Consider the nucleus Na™ which after §-radiation also 
emits a y-ray or the internally converted electron (e~- 
radiation). According to Konopinski,’ this y-ray may have 
a high multipole character and therefore may have a 
rather long life (about 10~? sec. for octupole radiation). 
Thus after the 8-emission, some of the recoil nuclei will, 
if they are allowed to move freely by placing the source 
and accessories in a high vacuum, travel a distance of the 
order of one centimeter before the emission of the y- 
radiation or the e~-radiation. This fact can be utilized to 
detect the neutrino by applying either the cloud chamber 
or “directional e~-coincidence counting” technique. By 
“directional coincidence" we mean that the relative 
positions of the counters, the source, and the slits are so 
arranged that coincidence counting can be produced only 
when the direction of the 8-emission and that of the recoil 
nucleus, which subsequently emits an e~-particle, make a 
certain definite angle @. If @ is always found to be 180°, it 
will be definite that no neutrino has been emitted during 
the 8-disintegration. Otherwise, conservation principles 
can be used to find the energy and momentum of the 
neutrino, which can then be statistically checked with the 
Fermi theory. Even if Na™ is not suitable for our purpose, 
it is very probable that, for certain By (not y$)-radiative 
substances, the y-radiation lifetimes lie between 107? and 
10~* sec. For these substances, there will be a high proba- 
bility of internal conversion, thus favoring this method. 
(It is of interest to note that by measuring the length of 
the recoil path, the lifetime of such y-radiation’can be 
determined.) Besides, certain fission fragments, which 
emit 6-rays quickly and successively, may also be useful. 
Finally, the natural Sa-radioactive substance ThCC’ may 
also be considered although the recoil mass is, in this case, 
too large, and the troublesome “aggregate recoil” may 
occur. After all, the possibility of carrying out experiments 
like these depends on the preparation of a “‘monolayer”’ 
source, since scattering process would probably confuse 
the results. 

A third method is to utilize the K-radioactive property 
of certain substances like Be’, Ca“, Cr®, etc. For some of 
these nuclei the K-capture may be followed by y-radiation 
of rather long lifetime. By placing the source in vacuum, 
the two isomeric forms can be separated and tested for 
recoil by K-capture. The success of this separation will 
unambiguously prove the existence of the neutrino. 

It may be mentioned that the investigation of the 
neutrino would also be possible from the meson disintegra- 

















tion. The general assumption, that positive mesons always 
decay into positrons and neutrinos while negative mesons 
are mostly absorbed by nuclei, cannot be regarded as well 
founded.* It seems to the writer that a coupling between 
a special counter arrangement similar to that designed by 
Rasetti’ and a cloud chamber will be suitable for this 














investigation. 
1K. C. W Phys. Rev. 61, 97 ie). 
Se oe ‘Allen, Phy hys. a. 6l, 692 (194: 2). 
Rev. Mod. Phys. 15, 209 (1943 


4G. bE ne Ree S. Klaiber, Phys. Rev. 70, Is (1946). 
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Y. Nishina, M. Takeuchi, and T. a Phys. Rev. 55, 585 (1939) 
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Velocity of Propagation of the Discharge in 
Geiger-Miiller Counters 


ERNEST WANTUCH 
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N important factor limiting the counting-rate readily 
achievable with counters is the time required for the 
avalanche to spread along the length of the central wire. 
When an ionizing event occurs in the volume of the 
counter, the electrons will be drawn toward the positive 
central wire. Near the central wire the electrons will 
acquire sufficient energy to ionize gas atoms by collision. 
During the chain of ionizations by collision, recombination, 
and excitation, photons, as well as new secondary electrons, 
are produced. The photons travel out in all directions, 
and eject electrons, thereby spreading the discharge down 
the length of the counter. The spread velocity of the 
discharge was measured by two experimental arrange- 
ments. 

1. The first arrangement consisted of a counter with a 
divided cylinder. The main section of the copper counter 
cylinder was 30 cm long and 1 cm in diameter. At each 
end was placed another section of the same diameter, 1 cm 
long. The separation between the sections was less than 
0.1 cm. One central wire served as anode for the three 
sections. 

The outputs from the end sections were applied to a 
timing circuit.! The resulting pulse which represented the 
time interval between the discharges of the two end 
sections was applied to the vertical plates of a DuMont 
208 oscilloscope. The background of the counter (approxi- 
mately 200 counts per minute) was used as our source of 
impulses. This resulted in a random distribution of time 
intervals which were observed on the oscilloscope. The 
maximum time-duration resulted from a discharge initiated 
at one end of the counter and travelling the length of the 

4 counter to be recorded at the other end. The timing circuit 
was Calibrated by a pulse generator giving a pair of pulses 
of known time separation. 

A critical analysis of this experimental arrangement 
showed that the measurements represented the time taken 
by the discharge to spread along the length of the counter 
plus any fluctuations in the time required to collect the 
positive ions at the cylinder. For this reason, an arrange- 
ment designed to minimize these fluctuations was devised. 
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TABLE I. Velocity of propagation of the discharge, 











— 
Gas (cm of Hg) _Max. time interval Spread velocity _ 
Method 1 aan 
0.6 cm ethyl alcohol 5.7 sec. 5.4 
9.4 cm argon aes cm/sec, 
1.0 cm ethyl alcohol 6.7 4.6 X106 
9.4 cm argon 
2.0 cm ethyl alcohol 8.0 3.9 X10 
8.0 cm argon 
y 4 cm ethyl alcohol 8.9 3.5 X10 
cm argon 
Method 2 
1.0 cm ethyl acetate 4.5 sec. 6.2 X10* cm /sec, 
8.0 cm argon ; 








a 

2. A second counter was constructed consisting of g 
cylinder 30 cm long and 1 cm diameter. Small circulg 
probes around the central wire were made of No, % 
platinum wire, and were of 0.14 cm inside diameter, The 
probes were located 28 cm apart. 

Pulses obtained from these probes were separately 
amplified by one stage of 6SJ7 tubes, then combined ing 
suitable differentiating network. This resulted in a double 
pulse with the separation between the peaks representing 
the travel time of the discharge. The double pulse was 
further amplified by a video amplifier and applied to the 
vertical plates of a TS-28 synchroscope. The maximum 
time difference between the peaks can be interpreted in 
terms of the time scale of the sweeps. The experimental 
results are summarized in Table I. These results support 
the theory that the positive-ion sheath spreads by photon 
emission and ionization. When the proportion of organic 
vapor is increased, the travel time should increase as the 
organic molecule should absorb these photons and then 
pre-dissociate, thus removing the photon from the chain 
of ionizations. The fluctuations in the time required to 
collect positive ions at the cylinder must be comparatively 
small. Our results show a slower propagation velocity than 
results previously reported,? Huber, e¢ al., having found 
velocities approximately 50 percent higher than the values 
reported here for corresponding mixtures. One serious 
objection to their first method is the finite width of the 
pulses applied to the electrical shutters inside their counter, 
In their second method the sharp dependence of their 
experimental values on applied voltage cannot be ac 
counted for. 

The author wishes to thank Professor S. A. Korff for 
helpful suggestions. 

1B. Rossi ong N. Nereson, Rew, Sci. Inst. 17, 65 (1 
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(19 P. a . Alden, and E. Baldinger, Helv. Phys. Acta 19, 24 
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Direction of Domain Magnetization in 
Powder Patterns 


H. J. WILLtiamMs 
Bell Telephone Laboratories, Murray Hill, New Jersey 
April 7, 1947 


N the course of an investigation of magnetic domain 
patterns, using colloidal particles of magnetite accord- 
ing to the well-known techniques,’ a number of new 
patterns have been observed, and the directions of magnet 
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pe. my Typical powder pattern on a (100) surface. Magnification 


Fic. 2, Pattern on a (100) surface with a set of parallel scratches 
made aru engine. 
Fic. 3. Pattern on a (100) surface with horizontal and vertical 
scratches made with a glass fibre brush. 

Arrows indicate directions of magnetization of the domains 
in 

nf 5. Modification of pattern in Fig. 1 by tension, within elastic 
limit, applied in a horizontal direction. 


a) Cavity with magnetic poles on the surface. (b) Domain 
1 sed by Néel. 


Fic. 7. Actual pattern around cavity. 
Fic. 8. Pattern around another cavity. 


ization in these patterns have been determined. Some of 
the results of the incomplete investigation are presented 
here. 

The (100) face of a single crystal specimen of iron 
containing 3.8 weight percent silicon was polished electro- 
lytically to give a strain-free surface. The pattern obtained 
is of a new and readily interpreted nature, and the inter- 
pretation has been checked by several new experiments. 
Figures 1, 2, and 3 are examples of the pattern for the 
demagnetized condition. The interpretation (Fig. 4 corre- 
. sponding to Fig. 3) shows the domains as magnetized in the 
directions of easy magnetization, < 100> for silicon iron, 
with the pattern arranged so as to minimize the poles 
produced on the boundaries in accordance with theory.* ® 

This interpretation is confirmed by two new experiments 
involving (1) scratching the surface, a technique suggested 
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by Dr. W. Shockley, and (2) applying mechanical stress. 
In Fig. 2 the horizontal dashes show portions of a series 
of parallel continuous scratches made by a conical sapphire 
point and a ruling engine, and reduced in intensity by an 
electrolytic polish. In the tapered regions the scratches 
show strongly; according to Fig. 4, in these regions the 
flux is perpendicular to the scratch; consequently, some 
flux leaves the metal and attracts the colloid. Where the 
scratch is parallel to the flux, no field emerges, the colloid 
is not collected, and the scratch is faint or imperceptible. 
Figure 3 shows the same effect with scratches running 
horizontally and vertically made by a brush of fine glass 
wool. 

Applying mechanical tension horizontally favors the 
horizontally magnetized domains from an energetic view- 
point since silicon iron has positive magnetostriction. 
When stress is applied, the tapered regions of Fig. 1 
(which are unfavorably magnetized) are seen to shrink, 
and the pattern of Fig. 5 is obtained. On releasing the 
stress, tapered regions reappear. All of this confirms the 
interpretation given in Fig. 4. 

The tapered regions of the pattern arise because . the 
face is not accurately parallel to (100). This phenomenon, 
together with the effects of applied magnetic fields, will be 
discussed at a later date. 

A further aid to the interpretation of Fig. 5 is furnished 
by the vertically elongated splotches. These are thought 
to be caused by roughness of the surface (prolonged 
electrolytic polishing produces a lemon-peel appearance) 
and act through the same mechanism as the scratches. 
Unlike the domain boundaries, the splotches appear 
always in the same place. They show up best when a 
concentrated colloid is used and serve, as the scratches ™ 
to interpret direction of magnetization. 

The effect of a hole in the silicon iron is shown in Figs. : 
and 8. Patterns of this kind, as shown in Fig. 6, have been 
predicted by Néel.* The magnetostatic energy caused by 
the concentration of poles on the surface of an inclusion 
(Fig. 6a) is reduced by the domain structure (Fig. 6b) 
which distributes the poles along the curved portions of 
the boundaries. As the domain elongates the magnetostatic 
energy decreases and the wall energy increases, and that 
structure occurs for which the sum of the two energies is 
a minimum. 

iF. Bitter, Phys. Rev. 38, 1903-1905 (1931); 

*L. bg H4mos and P. A es Zeits. 

*L. McKeehan Fy 'W. C. Elmore, 
529-531 1 (1988). 

«W. C. Elmore, Phys. Rev. 62, 486-493 (1942). 
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